[MITgcm-support] Using pkg/calendar with diagnostics

[Dimitris Menemenlis]

Hi Samar, I am not very familiar with the timeInterval logic.
But I if you look at the time stamps of your output files,
if you start integration from January 1, they should be
31, 59, 90, 120, 151, etc. days for normal years and
31, 60, 91, 121, 152, etc. days for leap years.
Maybe the "timeInterval" computation is off?

To get 2 years of output you need to make sure
that endTime >= 31536000 for normal years
and endTime >= 31633400 for leap years.

On Mar 21, 2014, at 10:38 AM, Samar Khatiwala <spk at ldeo.columbia.edu> wrote:

> Thanks Dimitris!
> 
> I gave it a try but am puzzled by the output. I ran the model for 2 years expecting to find 24 output files 
> but instead got 23. Also, the timeInterval stamp in the .meta files looks odd. For example, compare the 
> second digit in the 2d line with the first digit in the 3d line. There's a gap which seems to explain why 
> none of the averaging intervals exceeds 29.5 days (except the first one). Any thoughts as to what I'm 
> doing wrong?
> 
> Thanks, Samar
> 
> timeInterval = [  3.600000000000E+03  2.635200000000E+06 ]; 30.4583  <-- number of days
> timeInterval = [  2.635200000000E+06  5.140800000000E+06 ]; 29
> timeInterval = [  5.270400000000E+06  7.819200000000E+06 ]; 29.5
> timeInterval = [  7.905600000000E+06  1.041120000000E+07 ]; 29
> timeInterval = [  1.054080000000E+07  1.308960000000E+07 ]; 29.5
> timeInterval = [  1.317600000000E+07  1.568160000000E+07 ]; 29
> timeInterval = [  1.581120000000E+07  1.836000000000E+07 ];