[MITgcm-devel] 'stop in exf_interp.F: value of filePrec not allowed'

[Patrick Heimbach]

Martin,

On Dec 15, 2006, at 9:24 AM, Martin Losch wrote:

> Patrick, sorry, me again:
> if the initialization of AREA(:,:,3,:,:) cause (small) problems I  
> would be happy to change this to a local array (makes more sense to  
> me anyway):
>
> instead of
> AREA(I,J,3,bi,bj) = MAX(A22,AREA(I,J,2,bi,jb)
> do
> areaLoc(I,J) = MAX(A22,AREA(I,J,2,bi,jb)
> and leave the field AREA alone at this point (it's not really clean  
> to abuse AREA(3) anyway).
>
> What do you think?

Yep, would be cleaner.

Re. heffm, it's a good question.
It either points to a problem, or,
since heffm is passed through several levels of S/R calls,
TAF doesn't keep track, and conservatively assumes activity
(I've seen such behaviour in the past).
Not sure...

Looking in the adjoint code, it looks benign, i.e.
it's really only read-fro/written-to common once
(no other occurences such as generation of adheffm).

-p.



> Martin
>
> On 15 Dec 2006, at 15:16, Martin Losch wrote:
>
>> Patrick,
>>
>> why do you have to store HEFFM in seaice_model.F (this is the land  
>> mask for the seaice model and should really be replaced by maskC 
>> (i,i,1,bi,bj)). Doesn't this point to something strange?
>>
>> I just hope that is wasn't too difficult the fix the adjoint ...
>>
>> Martin
>>
>>
>>
>>
>>
>> _______________________________________________
>> MITgcm-devel mailing list
>> MITgcm-devel at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-devel
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach