[MITgcm-support] diagnostics parameter

Jody Klymak jklymak at uvic.ca
Tue Sep 10 12:46:09 EDT 2013


Hi Martin,

Before Jean-Michel explains it properly:

On Sep 10, 2013, at  9:28 AM, Martin Losch <Martin.Losch at awi.de> 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    ', ….

I think you want frequency(1)=1200.  
I'm pretty sure that timePhase should be less than 1200. i.e. it says where in the 1200s the averaging occurs.  

C     frequency :: Frequency (in s) of Output (ouput every "frequency" second)
C     timePhase :: phase (in s) within the "frequency" period to write output

I am very confused by your math.  141192*1200sec = 169430400sec=1961 days, not 134, and none of your numbers below make sense ;-)

I don't know that there is a way to start writing diagnostics after a certain time w/o restarting to model.

Cheers,   Jody


> 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

--
Jody Klymak    
http://web.uvic.ca/~jklymak/







More information about the MITgcm-support mailing list