[MITgcm-support] Sensitivities, OpenAD, tiles
Martin Losch
Martin.Losch at awi.de
Fri Sep 4 03:13:54 EDT 2015
Hi Heriberto,
I haven’t worked with OpenAD yet, but if your model simulations does not crash, then I assume this this is a benign problem associated only with the output. I had a similar problem with adjoint seaice variables (and TAF code), that weren’t properly exchanged before writing to disk. Please have a look at pkg/autodiff/addummy_in_stepping.F which is the routine that organizes the output of the adjoint variables. You’ll see that there are plenty of “call adxexch”, but unfortunately only if ALLOW_OPENAD is not defined. I guess you can probably try and fix the code for the OpenAD case youself (and contribute your fixes).
Martin
> On 04 Sep 2015, at 02:48, Heriberto Vazquez <heriberto1mx at gmail.com> wrote:
>
> Hello everyone
>
> I am doing some research in the Gulf of Mexico, so far I have an implementation of MITgcm and its adjoint built with OpenAD. Everything looks just fine except because in the sensitivity plots I can see the tiles' edges with relative high values in sensitivities and in the corner of the tiles with the highest values.
>
> Some peculiarities are the fact that the sensitivity in adjetan works just fine without problem with the edges but in sensitivity of adjuvel and adjvvel I can see the lines; for adjuvel are more visible the vertical lines (north-south) and for adjvvel are more visible the horizontal lines (east-west), but for salinity and temperature both lines (horizontal and vertical) are visible and in the corners of the tiles appear the biggest sensitivities.
>
> I tried changing advection schemes and the size in the overlapping but the lines are still there.
>
> It looks like is an exchange problem, but I don't know where to look at.
>
> Do you have any recommendations or suggestions?
>
> Thanks in advance for your help
>
> Cheers,
>
> Heriberto
> Posdoctoral Researcher
> Scripps Institution of Oceanography
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list