[MITgcm-devel] depressed Eta and sIceLoad

[Martin Losch]

Hi Jean-Michel,

unfortunately I cannot say when, but at sometime the sIceLoad fields  
started to cause problems in the height field Eta. With useSeaice  
= .true., useRealFreshWaterFlux=.true., exactConserv=.true., and  
LINEAR free surface (#undef NONLIN_FRSURF) I get local depressions of  
ETA underneath seaice  along the coast of -40m and more in 100years.  
They go away, when I set useRealFreshWaterFlux = .false. I need to  
mention, that I have a global EmPmR imbalance of approximately 10cm  
per year, so that my mean Eta is approximately -11m after 100years,  
but that wasn't a problem earlier. I attach a figure which shows Eta  
and surface velocities averaged over the years 91 to 100 of a 100  
year integration. What I find strange is that the strong depression  
in Eta does not have any corresponding signal in the surface  
velocities. T/S do not have any striking anomalies in these regions  
around Antartica, either, that could compensate Eta-gradients of more  
10m per 100km. The model does not explode, the mean effective ice  
thickness is way below 2m in the areas of Eta-depression.

Do you have any idea what's going on?

Martin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sload.png
Type: image/png
Size: 29311 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20061023/de4b56de/attachment.png>