[MITgcm-devel] problem with TAF tang-Lin exp. hs94.1x64x5

[Martin Losch]

Hi Jean-Michel,

I tried reverting the code to the previous version, but that does not help. Going back to version 5.8.0 does not either, not the combination of both.
Going back to pre PR 517 works. 
Using “staf -keep” works, but we don’t want that.
That’s weird, so TAF is confused.

This works:

        ilgen=ilnblnk( xx_gentim2d_file(iarr) )
        temp_genarr_fnamA = xx_gentim2d_file(iarr)

        call ctrl_init_rec ( temp_genarr_fnamA(1:ilgen),
[…]

but does it work with OpenAD/Tapenade?

Looking at the code the character array "fnamebase(1:maxCtrlTim2D)” is not necessary at all, definitely it does not need to be an array. Maybe that’s what upsets TAF?

Martin

BTW, the second assignment of fnamebase(parr) can be dropped, too, it’s not used after that. I say we drop the entire variable.

> On 21. Sep 2023, at 17:37, Jean-Michel Campin <jmc at mit.edu> wrote:
> 
> Hi Martin and Dan,
> 
> With the changes in ctrl_map_ini_gentim2d.F from Dan's PR #751 (got merged yesterday)
> the TAF tang-lin experiment "hs94.1x64x5" does not run anymore with gfortran and -devel
> see, without MPI:
>  http://mitgcm.org/testing/results/2023_09/tr_villon-b_20230921_1/summary.txt
> and with MPI:
>  http://mitgcm.org/testing/results/2023_09/tr_engaging-gfoTlm_20230921_0/summary.txt
> 
> It looks like a TAF problem, current Version 5.8.8, which remove lines 122 & 124:
> 122         ilgen=ilnblnk( xx_gentim2d_file(iarr) )
> 123         temp_genarr_fnamA = xx_gentim2d_file(iarr)
> 124         fnamebase(iarr) = temp_genarr_fnamA(1:ilgen)
> in the Tang-Lin version of this routine (i.e., in ctrl_map_ini_gentim2d_tl.f),
> just before the call to "ctrl_init_rec".
> And since this experiment is using "debugLevel = 2" (the only TLM test with debugLevel > 1), 
> it triggers an error in ctrl_init_rec (with -devel, local var fnamebase(iarr) is not set at all).
> 
> With the older version of ctrl_map_ini_gentim2d.F (before PR #751), TAF was keeping
> these 2 lines:
>>          ilgen = ilnblnk(xx_gentim2d_file(iarr))
>>          fnamebase(iarr) = xx_gentim2d_file(iarr)(1:ilgen)
> 
> If you know a trick to convince TAF to keep these 2 lines, let me know.
> Otherwise, will likely submit the problem to Ralf Giering.
> 
> Cheers,
> Jean-Michel
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel