[MITgcm-support] Diagnostics at end of run

Roland Young Roland.Young at physics.ox.ac.uk
Mon May 11 06:11:16 EDT 2015


Hi,

I’m having trouble outputting diagnostics at the end of a run. I’m currently running a test job with time step 600s for 12 hours (72 timesteps). I want to get diagnostic output of some instantaneous fields after 6 hours and 12 hours, but I can’t get it to output diagnostics in the final timestep.

Here are the relevant lines in my data.diagnostics:

 &DIAGNOSTICS_LIST

  dumpAtLast = .true.,
  diag_mnc   = .true.,

  fields(1:4,1) = 'UVEL', 'VVEL', 'WVEL', 'THETA',
  filename(1)   = 'diagout4',
  frequency(1)  = -21600.,
  timePhase(1)  =  21600.,

 &
 
&DIAG_STATIS_PARMS
 &

This gives me output at timesteps 0 and 36 (model start and 6h).

timePhase(1)  =  0. also gives me output at timesteps 0 and 36.
timePhase(1)  =  21000., gives me output at timesteps 35 and 71 (i.e. 1 time step before the last one).
timePhase(1)  =  600., gives me output at timesteps 1 and 37.
timePhase(1)  =  10800., gives me output at timesteps 18 and 54 (3h and 9h).
Leaving timePhase(1) unset also gives me output at timesteps 18 and 54.

None of which is what I want. What am I doing wrong here?

Thanks,

Roland






More information about the MITgcm-support mailing list