[MITgcm-devel] multiple ptracers with divided adjoint

[Patrick Heimbach]

Hi Matt,

yes, seems to be a problem with TAF,
and how it generates the files.
I'll get in touch with Fastopt.

Thanks for testing
-p.

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