[MITgcm-devel] depressed Eta and sIceLoad

[Patrick Heimbach]

Hi Martin,

I implemented some initial diagnostics
(i.e. using pkg/diagnostics)
for the seaice package over the weekend.
It does include HSNOW.

So you could have a look at it.
As template for the use of diagnostics
check out labsea/input/
which has two seaice variables as part
of the diagnostics list
(you need to add 'SIhsnow ').

You need code updates in pkg/seaice, pkg/diagnostics
and model/src.

Hope it works.

-p.



On Mon, 2006-10-23 at 10:50, Jean-Michel Campin wrote:
> Hi Martin,
> 
> My first guess would be the snow-thickness (unfortunately, there was 
> not many diagnostic for this field in pkg/seaice, may be Patrick
> add one), since in the run I did, there was quiet a lot of snow.
> Cheers,
> Jean-Michel
> 
> 
> On Mon, Oct 23, 2006 at 10:26:09AM +0200, Martin Losch wrote:
> > 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
> 
> 
> > _______________________________________________
> > MITgcm-devel mailing list
> > MITgcm-devel at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
-- 
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach