[MITgcm-support] shortwave, water type and calendar package
Patrick Heimbach
heimbach at MIT.EDU
Sun Apr 18 14:17:04 EDT 2004
Approved.
And no comments, please! ;o)
-p.
Quoting Martin Losch <mlosch at awi-bremerhaven.de>:
> Hi,
> it's me again. As usual I sent my email before I figured the problem
> out myself:
> myIter is not passed to swfrac but it is set to 0. If all dates are
> set correctly in the calendar package
> (theCalendar='model',startDate_1=00000101, in my case), this value of
> myIter makes cal_getMonthsRec return count0=12 and count1=1, so that
> jwtype(jerl(count0))=2, and everything is fine. But my question
> remains: It's not very transparent that the date is computed to be the
> same at all times, but at the same time the cal-package is tricked into
> always return the same values. I would suggest to set jwtype=2
> explicitly, comment out the call to cal_getMonthsRec and write a proper
> annotation. I can volunteer to do that, but I cannot promise that my
> comments won't be sarcastic (o:
>
> Martin
>
> On Sunday, April 18, 2004, at 01:42 PM, Martin Losch wrote:
>
> > Hi,
> > I am trying to run the MITgcm with dynamic seaice and climatological
> > forcing (12 monthly mean records). All sorts of problem pop up but
> > most prominently the model stalls with an address error in swfrac.F.
> > The reason is that with seaice you need the exf-package and
> > consequently the cal-package; then swfrac calls cal_getMonthsRec,
> > which for my parameter settings (theCalendar='model') screws up and
> > returns count0=0, so that the statement jwtype=jerl(count0) is not
> > defined (jerl is defined as "integer jerl(12)", so count0 should range
> > from 1 to 12).
> > Now, I could go ahead and check, why cal_getMonthsRec gives the wrong
> > answer (and I will try to do that, but mainly because this routine
> > also screws up in different place). But before I do that I would like
> > someone to explain to my why jerl (which obviously has to do with 12
> > months) is a globally uniform constant. Does the water type (and
> > consequently the absorption properties of sea water) change uniformly
> > over the globe over the year? Hard to believe, but I have to admit
> > that I don't know much about this water type business.
> > I would think that I would make more sense to remove the call to
> > cal_getMonthsRec from swfrac by default and have the water type be a
> > global constant. If someone wants to include some dependence on the
> > season (I guess biological production is a candidate that could change
> > the water type), it can be done for the specific problem.
> > Any thoughts on this? Where am I wrong?
> >
> > Martin
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>
_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/
Patrick Heimbach Massachusetts Institute of Technology
FON: +1/617/253-5259 EAPS, Room 54-1518
FAX: +1/617/253-4464 77 Massachusetts Avenue
mailto:heimbach at mit.edu Cambridge MA 02139
http://www.mit.edu/~heimbach/ U.S.A.
More information about the MITgcm-support
mailing list