[MITgcm-support] diagnostics at last timestep

[Patrick Heimbach]

Hi Martin,

I looked at the code diagnostics_write.F
and it seems clear that dumpAtLast (which is supposed to do the trick)
is only implemented in conjunction with positive frequency,
i.e. time-averages.
Don't know why that is, and don't see a good reason for it.

If you want to be courageous and test, then change line

           IF ( dumpAtLast .AND. myTime.EQ.endTime
      &                    .AND. freqSec.GE.0. ) dump2fileNow = .TRUE.

(remove ".AND. freqSec.GE.0."),
or wait for the author of diagnostics to reply.

-p.



On Mar 1, 2007, at 12:20 PM, Martin Losch wrote:

> Hi there,
>
> I can't make the diagnostics package dump a snapshot of a  
> diagnostics variable, say UVEL, at the last time step. I have tried  
> all sorts of things:
> frequency(1) = -3600.,
> timePhase(1)= 0.,3600., nothing
> dumpAtLast = .true.,
> but nothing works and I don't get output for the last timestep  
> (which is at a multiple of 3600.sec of course). For time averaged  
> quantities (frequency>0) it seems to work.
> What have I overlooked?
>
> Martin
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

---
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach