[MITgcm-devel] multiple ptracers with divided adjoint

Thu Mar 24 17:30:23 EDT 2011

Hey,
sorry, I had changed something else in my testing, and forgot to undo.
Nothing to do with DIVA.
Will update repository in a minute.
-p.

On Mar 24, 2011, at 4:17 PM, Matthew Mazloff wrote:

> Hi Patrick,
> 
> I just tried compiling and received the error:
> 
> TAF ERROR   ad_input_code.f:292085 : tape not defined
> 
> Though there are several other errors I received...see attached taf_ad.log
> 
> The only difference I can think of between our setups is that I am compiling the DIVA.  Do you have any ideas regarding these errors?
> 
> Thanks
> -Matt
> 
> 
> <taf_ad.log>
> 
> 
> 
> 
> On Mar 24, 2011, at 5:23 AM, Patrick Heimbach wrote:
> 
>> 
>> Hi Matt,
>> 
>> I've updated the ptracers_ad_check_lev?_dir.h
>> with working syntax (suggested by Ralf).
>> Please update and set
>> #define AUTODIFF_PTRACERS_SPLIT_FILES
>> in ECCO_CPPOPTIONS.h
>> I've tested it in tutorial_tracer_adjsens
>> 
>> I'm also appending what Ralf wrote back.
>> 
>> Cheers
>> -Patrick
>> 
>> ---
>> 
>> Hi Patrick,
>> 
>> The assumption is that a store-directive for a variable provides all
>> parts of that variable which are required.
>> So several store-directives for different parts won't work.
>> A solution is to store different parts of that array in one directive:
>> 
>> CADJ STORE ptracer(:,:,:,:,:,1),ptracer(:,:,:,:,:,2),ptracer(:,:,:,:,:,3)
>> CADJ &     = tapelev2_ptr, key = ilev_2
>> 
>> or
>> 
>> CADJ STORE ptracer(:,:,:,:,:,1:5),ptracer(:,:,:,:,:,6:nTracer)
>> CADJ &     = tapelev2_ptr, key = ilev_2
>> 
>> Another solution would be to modify ADOPEN/ADREAD/ADWRITE/ADCLOSE to write several
>> files in case of to large files.
>> 
>> Cheers,
>> Ralf
>> 
>> 
>> 
>> On Mar 18, 2011, at 9:00 PM, Matthew Mazloff wrote:
>> 
>>> Hi Patrick,
>>> 
>>> I tried this (and several other similar things).  Taf just keeps dropping the store commands...they will be in ad_taf_input.f but gone from d_taf_input_ad.f and ad_taf_output.f
>>> 
>>> Not sure why it won't store it....if you have any ideas please let me know
>>> 
>>> Thanks
>>> Matt
>>> 
>>> 
>>> 
>>> On Mar 14, 2011, at 10:19 PM, Patrick Heimbach wrote:
>>> 
>>>> Hi Matt,
>>>> 
>>>> ok.
>>>> It ought to be possible to remedy this
>>>> by having a loop over all tracers in
>>>> ptracers_ad_check_levN_dir.h   ( N> 1 )
>>>> something like
>>>> 
>>>>   DO iTracer = 1, PTRACERS_num
>>>> CADJ STORE ptracer(:,:,:,:,:,iTracer)
>>>> CADJ &     = tapelev2, key = ilev_2
>>>> etc...
>>>>   ENDDO
>>>> 
>>>> That should generate a separate file for each tracer.
>>>> 
>>>> -p.
>>>> 
>>>> On Mar 14, 2011, at 10:57 PM, Matthew Mazloff wrote:
>>>> 
>>>>> Hi Patrick,
>>>>> 
>>>>> Just to follow up on this issue -- the problem has to do with the size of the tapelevel files -- so its machine limitations, not the code.
>>>>> 
>>>>> For my setup I can run with up to  PTRACERS_num =  6 -- then at  PTRACERS_num =  7 the outer tape level files get too large and the model poofs.
>>>>> 
>>>>> -Matt
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> On Mar 7, 2011, at 7:41 PM, Matthew Mazloff wrote:
>>>>> 
>>>>>> Hi Patrick,
>>>>>> 
>>>>>> Well it poofed in the outer-most tapelevels for PTRACERS_num =  11.
>>>>>> But actually, it poofed somewhere else for PTRACERS_num =  2.
>>>>>> 
>>>>>> I am going to re-run the PTRACERS_num =  2 case to make sure I know
>>>>>> exactly where it poofed
>>>>>> 
>>>>>> I haven't tried it non-DIVA
>>>>>> 
>>>>>> Thanks!
>>>>>> -Matt
>>>>>> 
>>>>>> 
>>>>>> 
>>>>>> 
>>>>>> On Mar 7, 2011, at 5:17 PM, Patrick Heimbach wrote:
>>>>>> 
>>>>>>> Hi Matt,
>>>>>>> does it also poof in same setup, but non-DIVA?
>>>>>>> Also, does it crash writing tapes or the DIVA snapshots?
>>>>>>> Also, does it crash writing outer-most tapelevels or intermediate?
>>>>>>> -p.
>>>>>>> 
>>>>>>> On Mar 7, 2011, at 7:29 PM, Matthew Mazloff wrote:
>>>>>>> 
>>>>>>>> Hello,
>>>>>>>> 
>>>>>>>> I tried to run the DIVA with multiple tracers but it crashes while
>>>>>>>> writing the tapelevels.  With PTRACERS_num =  1 everything is fine
>>>>>>>> --- I switch to PTRACERS_num =  2 and poof.
>>>>>>>> 
>>>>>>>> I'll look into whats going on later -- but wanted to check if
>>>>>>>> anyone was already aware of this problem.
>>>>>>>> 
>>>>>>>> Thanks
>>>>>>>> Matt
>>>>>>>> 
>>>>>>>> 
>>>>>>>> 
>>>>>>>> _______________________________________________
>>>>>>>> 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
>>>>>>> 
>>>>>>> 
>>>>>>> 
>>>>>>> _______________________________________________
>>>>>>> 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
>>>> 
>>>> 
>>>> 
>>>> _______________________________________________
>>>> 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
>> 
>> 
>> 
>> _______________________________________________
>> 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