[MITgcm-devel] initialization of MNC -- where and when?

[Alistair Adcroft]

Currently in write_state:

  call MNC_CW_RL_W_D('state',0,0,'iter',-1,iter,myThid)
  call MNC_CW_RL_W_D('state',0,0,'model_time',0,iter,myThid)
  call MNC_CW_RL_W_D('state',0,0,'U',0,iter,myThid)
  call MNC_CW_RL_W_D('state',0,0,'V',0,iter,myThid)
  call MNC_CW_RL_W_D('state',0,0,'T',0,iter,myThid)
  ...

Wew suggested:

  call MNC_CW_SET_RECDIM('state',record_number,myThid)
or
  call MNC_CW_INCR_RECDIM('state',myThid)

  call MNC_CW_RL_W_D('state',0,0,'iter',myThid)
  call MNC_CW_RL_W_D('state',0,0,'model_time',myThid)
  call MNC_CW_RL_W_D('state',0,0,'U',myThid)
  call MNC_CW_RL_W_D('state',0,0,'V',myThid)
  call MNC_CW_RL_W_D('state',0,0,'T',myThid)

MNC_CW_SET_RECDIM specifies which time-level (record number) to write in.
MNC_CW_INCR_RECDIM places the record at the end of the record dimension
All subsequent calls to MNC for files "state.*" use that record number until
another call to change the record number is made.

We don't have to delete the iter argument (but it isn't necessary).

A.
--
Dr Alistair Adcroft            http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative        tel: (617) 253-5938
EAPS 54-1523,  77 Massachusetts Ave,  Cambridge,  MA,  USA