[MITgcm-devel] partial cell under ice-shelf

[Jean-Michel Campin]

Martin,

For now, it would be with Linear Free-Surf only, to uncomment the line 68 
in set_params.F:
c     uniformFreeSurfLev = usingZCoords .AND. .NOT.useShelfIce
and (after changing the phi0surf file), to see if this works and
if results make sense. But there is no rush, when you find time.

Cheers,
Jean-Michel

On Wed, Apr 11, 2012 at 04:46:47PM +0200, Martin Losch wrote:
> Hi Jean-Michel,
> 
> not quite sure what I am supposed to test:
> On Apr 11, 2012, at 4:24 PM, Jean-Michel Campin wrote:
> 
> > Now regarding a real test with isomip set-up, I have the impression
> > that some of the input files would need to be changed (in particular
> > phi0surf.exp1.jmd95z ) to also account for partial cell 
> > in order to really make it consistent with the new (uniformFreeSurfLev=F)
> > partial cell code in calc_phi_hyd.F (this is what you described in the
> > JGR paper). I can try to test something, but I am not sure I will
> > do it right. So, when you have some time, if you can check and test this
> > piece of code, it would be great. In the mean time, will keep
> > set_params.F as it is now, so that nothing changed.
> 
> set uniformFreeSurfLev=F and reruns some of the isomip experiments with modified phi0surf file? For the case of a linear free surface? non-linear free surface? 
> 
> I'll look into your new code tomorrow (or the day after tomorrow, I have an all-day meeting tomorrow).
> 
> Martin
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel