[MITgcm-devel] tices in global_ocean.cs32x15

[Martin Losch]

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