[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