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

[Patrick Heimbach]

Hong,

(at least) one more thing.
The following recommendations for new defaults
from tag-index are wrong:
   - SEAICE_advSnow = .true. is now the default
   - SEAICE_advSalt = .true. is now the default
   - SEAICE_advAge  = .true. is now the default
   - SEAICE_clipVelocities  = .false. is now the default
None of these flags exist in latest code
(hence my job on columbia just blew again).

Instead, they are called
   - SEAICEadvSnow = .true. is now the default
   - SEAICEadvSalt = .true. is now the default
   - SEAICEadvAge  = .true. is now the default
   - SEAICE_clipVelocities  = .false. is now the default

-Patrick



On Apr 14, 2009, at 12:22 PM, Patrick Heimbach wrote:

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