[MITgcm-devel] Problem in yersterday changes

[Patrick Heimbach]

Fixed (I think).
Was related to adding snowprecip for thsice adjoint
and not updating seaice adjoint.
-p.

On Oct 18, 2010, at 10:22 AM, Jean-Michel Campin wrote:

> Hi Martin,
>
> I've checked that the problem comes from changes in pkg/autodiff
> made on Sat. Oct 16 (all tests with TAF version 2.0.0).
> I am not going to do more tests now, because I don't know much about
> thoses changes.
> Can someone check those changes ? Or is it clear that they are OK
> in which case, might be a more tricky Pb.
>
> Cheers,
> Jean-Michel
>
> On Mon, Oct 18, 2010 at 12:51:56PM +0200, Martin Losch wrote:
>> Jean-Michel
>>
>> timeint = deltaTclock = modelstep. Modelstep is the "cal" variable  
>> name. it is set once in cal_set and then kept in a common block  
>> and it should 3600 in this experiment, right? cal_getmonthsrec  
>> calls cal_timeinterval(-modelstep,...) so I assume that there is  
>> some memory problem. TAF has been changed a lot recently, so it's  
>> very well possible that suddenly the order or arguments etc. might  
>> be changed in the adjoint routines, so that memory gets  
>> overwritten? Can you figure out what happens with older versions  
>> of TAF?
>>
>> Martin
>>
>> On Oct 17, 2010, at 5:28 PM, Jean-Michel Campin wrote:
>>
>>> Some update:
>>> I can reproduce the problem on my laptop with ifort and without MPI,
>>> and the problem is still there with seaice_growth.F v.1.89 from
>>> Oct.15 (looks like Gael's changes are not the problem).
>>>
>>> Here is what I get when I run the noseaice test with idb:
>>>> Program received signal SIGFPE
>>>> cal_timeinterval (timeint=-9.2559631349317831e+61, timeunit= 
>>>> (...), date=(...), mythid=1, .tmp.TIMEUNIT.len_V$10=4) at /home/ 
>>>> jmc/mitgcm/gcm_ifc/verification/lab_sea/bld_ad1p/ 
>>>> cal_timeinterval.f:628
>>>> 628         date(1) = int(timeint/float(secondsperday))
>>> I guess timeint/secondsperday is too big to fit into an integer.
>>>
>>> Jean-Michel
>>>
>>> On Sun, Oct 17, 2010 at 10:51:41AM -0400, Jean-Michel Campin wrote:
>>>> Hi,
>>>>
>>>> Looks like there is a problem with all lab_sea AD tests
>>>> as seen on aces with ifort+mpi:
>>>> http://mitgcm.org/testing/results/2010_10/tr_aces-adm_20101017_0/ 
>>>> summary.txt
>>>>
>>>> My impression is that it's not directly related to mpi (to be  
>>>> confirmed).
>>>>
>>>> This was before Martin's lates changes in seaice-solve4temp.F
>>>> and caused by one of yesterday modifs.
>>>>
>>>> I propose to hold on changes until this is fixed.
>>>>
>>>> Cheers,
>>>> Jean-Michel
>>>>
>>>> _______________________________________________
>>>> 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
>>
>>
>> _______________________________________________
>> 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

---
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