[MITgcm-support] Continuous decrease of sea surface height in model results

[Jonathan Whitefield]

I had the inverse problem (i.e. increasing ETAN) using a regional
configuration in the Arctic. If you are also using a regional setup,
the the chances are your EmPmR is imbalanced, and you are losing water
in your domain through evaporation. The workaround that I used was to
apply a land mask and remove a regional mean from the calculated ETAN
fields.

For example:

RAC=readbin('<directory>\RAC.data',[420 384]);
hFacC=readbin('<directory>\hFacC.data',[420 384]);
RAC=RAC.*hFacC;

tmp(hFacC==0)=NaN; %where tmp is the loaded ETAN field
tmp=tmp-nansum(tmp(:).*RAC(:))/sum(RAC(:));

%Then write out the tmp field and save it.

Hope this makes sense