[MITgcm-devel] partial cell under ice-shelf

[Jean-Michel Campin]

Hi Martin,

I move this to devel list since others might be interested.

I made some changes yesterday in calc_phi_hyd.F, to allow to account
for partical cell near the surface when using Z-coordinate (a case
relevant for pkg/shelfice), but I did not yet turn the switch on
(in set_params.F:
>       uniformFreeSurfLev = usingZCoords
> C- note: should be instead:
> c     uniformFreeSurfLev = usingZCoords .AND. .NOT.useShelfIce
)
I made a quick check, and I can set uniformFreeSurfLev=F
(to use the new calc_phi_hyd.F piece of code) and testreport results 
change only at the trunction level (which is expected), except 
the 3 isomip exp.

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.

I also found a couple of places in Non-Hydrostic code where I 
need to replace usingZCoords with uniformFreeSurfLev. Will do this
now before I forget.

Thanks,
Jean-Michel

On Tue, Apr 10, 2012 at 06:45:37PM -0400, Jean-Michel Campin wrote:
> Hi Martin,
> 
> I will try to make some changes in calc_phi_hyd.F and
> will see how it goes (might ask you to run a test at some point).
> Cheers,
> Jean-Michel