[MITgcm-support] caledar package
Jean-Michel Campin
jmc at ocean.mit.edu
Wed Jun 6 17:03:59 EDT 2012
Hi Martha,
David is right:
1) pkg/diagnostics shapshots ignores calendarDumps. I will try
to change this soon.
2) for pkg/diagnostics shapshots, the default phase (timePhase in
data.diagnostics) is half of the corresponding frequency (whereas
for time-averaged diagnostics, the default phase is zero), so that
if you want shapshots output with no phase lag, you will need to
specify explicitly in data.diagnostics: timePhase(n)=0.
But regarding this half day shift:
> >For the time average diagnostics, the real calender months works fine (except
> >for the first output is at 30.5 days rather than as it should be at 31
I don't know; might be useful to send your parameter files (data, data.cal,
data.diagnostics).
Cheers,
Jean-Michel
On Fri, May 18, 2012 at 01:51:58PM -0400, Daiwei (David) Wang wrote:
> Hi Martha,
>
> I can confirm that diagnostics snapshots do not conform to calendar
> months. I believe they are implemented evenly spaced in time.
>
> To save snapshots at the end (rather than the middle) of each
> interval, you need to set timePhase to be zero.
>
> To calculate (tracer) tendencies, I'd use TOTTTEND/TOTSTEND averages
> instead THETA/SALT snapshots. They can be true calendar monthly
> diagnostics.
>
> HTH,
> David
>
> On 5/18/12 12:34 PM, Martha W Buckley wrote:
> >I've found some strange behavior of the calendar when real months are output.
> >I want to output both time average diagnostics and snapshots monthly.
> >
> >For the time average diagnostics, the real calender months works fine (except
> >for the first output is at 30.5 days rather than as it should be at 31 so all
> >the outputs are shifted from what they should be by 12 hours: ie output #2 is
> >at (30.5+29)*24=1428 (since my first year happens to be a leap year) while it
> >should be at (31+29)=1440).
> >
> >However, when snapshots are output, the real calender months do not seem to be
> >used, which means that the snapshots aren't the right ones for calculating
> >tendencies. Even more strangely, the snapshots are centered at the middle of
> >the month-- first one is at 366=15.25*24, second is at 1098=(15.25+30.5)*24,
> >etc. I'd have expected them to be output at the beginning and end of the
> >month.
> >
> >thanks,
> >Martha
> >
> >_______________________________________________
> >MITgcm-support mailing list
> >MITgcm-support at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> --
> Daiwei (David) Wang, Postdoctoral Associate
> 54-1517, Program in Atmospheres, Oceans, and Climate
> Department of Earth, Atmospheric and Planetary Sciences
> Massachusetts Institute of Technology
> Tel: 617-253-5458
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list