[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