[MITgcm-support] Using pkg/calendar with diagnostics

[Samar Khatiwala]

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 ];
timeInterval = [  1.844640000000E+07  2.103840000000E+07 ];
timeInterval = [  2.108160000000E+07  2.363040000000E+07 ];
timeInterval = [  2.371680000000E+07  2.630880000000E+07 ];
timeInterval = [  2.635200000000E+07  2.890080000000E+07 ];
timeInterval = [  2.898720000000E+07  3.157920000000E+07 ];
timeInterval = [  3.162240000000E+07  3.425760000000E+07 ];
timeInterval = [  3.425760000000E+07  3.667680000000E+07 ];
timeInterval = [  3.689280000000E+07  3.935520000000E+07 ];
timeInterval = [  3.952800000000E+07  4.194720000000E+07 ];
timeInterval = [  4.216320000000E+07  4.462560000000E+07 ];
timeInterval = [  4.479840000000E+07  4.721760000000E+07 ];
timeInterval = [  4.743360000000E+07  4.989600000000E+07 ];
timeInterval = [  5.006880000000E+07  5.257440000000E+07 ];
timeInterval = [  5.270400000000E+07  5.516640000000E+07 ];
timeInterval = [  5.533920000000E+07  5.784480000000E+07 ];
timeInterval = [  5.797440000000E+07  6.043680000000E+07 ]; 28.5


On Mar 21, 2014, at 1:06 PM, Dimitris Menemenlis <dmenemenlis at gmail.com> wrote:

> Yes.  That's all you have to do if you want calendar-month output.  It works.
> 
> On Mar 21, 2014, at 9:59 AM, Samar Khatiwala <spk at ldeo.columbia.edu> wrote:
> 
>> I'd like to write out monthly mean diagnostics and I was wondering how I can accomplish that with the calendar 
>> package. Looking at diagnostics_write.F it would seem all I have to do is specify a nominal averaging frequency 
>> of, say, 86400*30.5 = 2635200.0 and set calendarDumps=.TRUE. in data.cal. But that seems suspiciously easy! 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support