[MITgcm-devel] seaice code beyond checkoint61j blows for ECCO-GODAE

[Patrick Heimbach]

Hong,

Jean-Michel gave me a few hints on how to restore
previous defaults to remain compatible with checkpoint61j

As far as I can see most of them are mentioned in Martin's email.
But I think one CPP option is missing:
I now put in SEAICE_OPTIONS.h
#define SEAICE_OLD_AND_BAD_DISCRETIZATION
Not sure which of the changes is the relevant one.
Currently have a job in queue on columbia to find out.

In the meantime, if you only revert pkg/seaice/
to checkpoint61j then you should be ok.

The changes to mom_calc_visc.F should be independent and
you can keep those (but also need consistent change to the_main_loop.F)

Good luck!
-Patrick



On Apr 14, 2009, at 12:04 PM, Hong Zhang wrote:

> On Tue, 2009-04-14 at 03:49 -0700, Martin Losch wrote:
>>> I do not know if Dimitris has gotten around to test this on CS-or  
>>> curvilinear grids.
> Hi
> We also got the problem on CS:
> got all seaice fields NaN both on levl 28 and level 50 config.
> We updated the code on April 7
> and want to fix the potential blowup (sensitivity) on CS-corner
> with Patrick's check-in (see :
>> Modified Files:
>>     mom_calc_visc.F
>> Log Message:
>> Overlaps had been forgotten in calculating ijk keys
>> (spotted by jmc using gfortran with check-bounds)
> Have to revert?
>
> cheers,
> hong
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
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