[MITgcm-devel] MITgcm-devel Digest, Vol 89, Issue 2

[Ian Fenty]

> From: Martin Losch <Martin.Losch at awi.de>
> Hi there,
> the problem described below also exists in our seaice_solve4temp.F, the "legacy" > branch.
...
>Do you think that we should correct this "bug". Our code tends to have too little ice in >t he Arctic, ...

In case this issue is not yet cleared up, let me throw my hat in the ring.  Bounding TSURF to T_melt during the iterative loop or letting it exceed T_melt should make no difference to the final result, provided that TSURF = max(T_melt, TSURF) after the iterative loop *and* all the fluxes which depend on TSURF are then recomputed.  

My guess is that Frank's code didn't recompute the fluxes after bounding TSURF following the iterative loop, leaving an inconsistency between the newly bound TSURF and the ice-atmosphere fluxes calculated with the TSURF which exceeded T_melt.  

In my opinion, letting TSURF exceed T_melt in the search was useful because I could verify that the solver was functioning.

-Ian

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20101206/4cf6946f/attachment.htm>