[MITgcm-devel] multiple ptracers with divided adjoint

[Matthew Mazloff]

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