[MITgcm-support] shortwave, water type and calendar package

[Martin Losch]

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