[MITgcm-support] Diagnostics at end of run
Roland.Young at physics.ox.ac.uk
Mon May 11 06:11:16 EDT 2015
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:
dumpAtLast = .true.,
diag_mnc = .true.,
fields(1:4,1) = 'UVEL', 'VVEL', 'WVEL', 'THETA',
filename(1) = 'diagout4',
frequency(1) = -21600.,
timePhase(1) = 21600.,
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?
More information about the MITgcm-support