[MITgcm-support] diagnostics parameter

Martin Losch Martin.Losch at awi.de
Tue Sep 10 16:18:20 EDT 2013


Hi Jean-Michel, thanks, I'll try that (frequency = 1200.,)
Martin

On Sep 10, 2013, at 9:52 PM, Jean-Michel Campin wrote:

> Hi Martin,
> 
> The pieces of code that decides when to write 
> 2D/3D diagnostics or stats-diags are the same, so it is strange that
> you get different behaviour between the 2. 
> Did you check (in STDOUT, just after reading data.diagnostics,
> the section starting with:
>> DIAGNOSTICS_READPARMS: active diagnostics summary:
> )
> that frequency and phase are really what you think they are ?
> 
> Apart from that, to get output only after 134 days, you need to 
> specify a frequency equal or larger than the time-step (when 
> frequency < timeStep, output is written every time-step 
> and timePhase is ignored).
> 
> Cheers,
> Jean-Michel
> 
> On Tue, Sep 10, 2013 at 06:28:23PM +0200, Martin Losch wrote:
>> Hi there,
>> 
>> this is a problem that I keep running into, and since the search function on mitgcm.org does not help very much (as usual), I have to ask this question again:
>> 
>> For debugging reasons, I want the model to start from a pickup at nIter0 = 131544. My timestep is deltaT = 1200. I want the model to write diagnostics at each timestep after timestep 141192, i.e. after the model has reached the time 141192*1200sec = 169430400sec=134days. I thought that I could do it like this in data.diagnostics (<http://mitgcm.org/pipermail/mitgcm-support/2010-August/006799.html>):
>> 
>>   filename(1) = 'myfile',
>>   frequency(1) = 1.,
>>   timePhase(1) = 169430400.,
>>   fields(1,1) = 'ETAN    ', ….
>> 
>> because it says in the standard data.diagnostics files in the verification experiments:
>> #  timePhase(n)     : write at time = timePhase + multiple of |frequency|
>> 
>> Obviously I got it wrong because now the model writes output at timesteps 132408 (after 12 days), 134568 (42days), 136800 (73days), etc.
>> Where did I go wrong, and how can I do it correctly?
>> 
>> Analogously in the diag_stats name list, I specified stats_freq=1 and stats_phase=169430400 and got output every time step from the very beginning (so already at time step 131544). Why are the two parameter sets behaving so differently?
>> 
>> Martin
>> 
>> 
>> 
>> 
>> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list