[MITgcm-devel] multiple ptracers with divided adjoint

Thu Mar 24 17:39:05 EDT 2011

Matt,

you actually have one more problem.
Your file 
obcs_ad_diff.list is out of date and missing the line
obcs_copy_tracer.f

which is the source of further errors you have.
This file got added to pkg/obcs in Nov. 2010.

-p.

On Mar 24, 2011, at 5:30 PM, Patrick Heimbach wrote:

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