[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