[MITgcm-devel] pickups: mdsio and mnc with unorthodox fields

Dimitris Menemenlis menemenlis at sbcglobal.net
Mon Mar 20 09:17:05 EST 2006


Hi Martin, I have never used mnc pickups (see Jean Michel's e-mail for that) but
I'm pretty sure that current seaice_pickup.F is OK.  In terms of 2-D fields,
this is the way it goes:

#ifdef SEAICE_MULTILEVEL
TICES    field 1-7
YNEG     field 8
HSNOW    field 9
UICE     field 10-12
...
#else /* SEAICE_MULTILEVEL */
TICE     field 1
YNEG     field 2
HSNOW    field 3
UICE     field 4-6
...
#endif /* SEAICE_MULTILEVEL */

> My understanding is, that the 7th argument of MDSWRITEFIELD give the position
> the file, so 4 (as for UICE) means the 4th record, but where are record 2
> and 3?

Records 1-3 (in terms of nNz=3 records) for SEAICE_MULTILEVEL are filled in by
TICES, YNEG, and HSNOW.

> Doesn't that mean, that UICE (second record) overwrites YNEG (second record
> also), equivalently VICE and HSNOW?

TICE, YNEG, and HSNOW have nNz=1.  So first three records fill the first three 
2-D levels.  UICE has nNz=3, so second record fills in 2-D levels 4-6.  So 
nothing is overwritten.

>  Is this correct? I want to add more fields, so I need to know.

Yes.  I have tested that it is OK using a
2timesteps = 1timestep + 1timestep.

Dimitris



More information about the MITgcm-devel mailing list