[MITgcm-devel] depressed Eta and sIceLoad

[Patrick Heimbach]

Martin,

just checked in some modifs.
Essentially we only want to output 1st
level of the seaice fields.
So try now.

Re. KPP changes, I thought diagnostics would/should
replace all the other dump stuff,
but supposedly it shouldn't.
So I fixed some naming conflicts to enable
all variables to be dumped in every possible
way at the same time.

HSNOW should now also be part of monitor.

-p.



On Mon, 2006-10-23 at 11:38, Patrick Heimbach wrote:
> Martin,
> 
> hold off a bit with updating.
> We're changing diagnostics yet again
> (issue over whether we want the 3 levels
> in AREA, HEFF, etc. or not).
> 
> -p.
> 
> 
> 
> On Mon, 2006-10-23 at 11:05, Patrick Heimbach wrote:
> > 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