[MITgcm-devel] tices in global_ocean.cs32x15

[Martin Losch]

On second thought, I do think we should do something about this and change SEAICE_multDim to 1 (default), this changes the results.

Martin

On Dec 10, 2014, at 5:40 PM, Martin Losch <Martin.Losch at awi.de> wrote:

> Hi Jean-Michel,
> 
> while looking for the problem that I introduced I stumble over this:
> 
> global_ocean.cs32x15 uses SEAICEmultDim = 7 and pickupStrictlyMatch = .FALSE., so the model assumes that the pickup files are OK and just reads them. nbFields = -1, and doMapTices = .FALSE. in seaice_read_pickup.F, but the pickup_seaice does not contain 7 seaice TICES fields as far as I can see. So that TICES(i,j,2,bi,bj) already only contains zeros and you get a division by zero in seaice_solve4temp (when computing mm_log10pi).
> 
> I guess it does not matter, because the 1/0 that are produced are overwritten very soon, because it’s an iteration, but I would how this can go undetected and not influence the results. Should we do anything about this?
> 
> Martin
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel