[MITgcm-devel] Re: [MITgcm-support] diagnostics at last timestep

[Martin Losch]

I could also threaten to change these lines in the repository to  
speed up the responses of the author of diagnostics (o:

But as a matter of a fact this "hack" does not work. Now I get out at  
the nIter0+nTimesteps-1st time step but not at the last one. I  
checked diagnostics_write.F again, and there a myItM1=myIter-1,
> C--     write snap-shot with suffix = myIter-1 to be consistent with
> C       state-variable time-step:
If I replace change that (wrITer=myIter and wrTime=myTime  also for  
freq<0) then I get output for all of my diagnostics at the last  
timestep (timePhase = 0,3600 and default). I guess that's what I  
want, but not for the case timePhase = -0.5*frequency (which is the  
default for snapshots).
Personally, I could live with this solution, but what does the  
diagnostics czar say?

Martin

On 1 Mar 2007, at 19:34, Patrick Heimbach wrote:

>
> 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
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support