[MITgcm-devel] [MITgcm-cvs] MITgcm/pkg/seaice CVS Commit

[Martin Losch]

Thank you, Jean-Michel,

to be honest, I didn’t check global_ocean.cs32x15 because I thought that lab_sea was the only experiment with adjoint dynamics. The changed results are expected, so it wouldn’t have made a difference. The rms of the gradient deviations changed from 3.98e-6 to 4.00e-6, I would call that “noise” in terms of agreement between AD and FD gradients.

The cost function in both lab_sea and global_ocean.cs32x15.seaice is sum(rA*HEFF) which only very indirectly tests the sensitivities of the drift velocities; they are still totally wrong, but the degree of wrongness is no reduced and the recomputations are all gone (for SEAICE_VECTORIZE_LSR/SEAICE_LSR_ADJOINT_ITER defined and also SEAICE_VECTORIZE_LSR_ZEBRA, we should probably think of new defaults). At some point I’ll have to talk to Patrick about this.

Martin

PS. I ***was*** going to generate and check-in updated results today after some more tests, but good that you’ve done it! Thanks again.

On Jul 5, 2014, at 4:54 AM, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:

> Hi Martin,
> I did not see any changes for TLM output.
> Going to check-in the 2 new output_adm.txt (before the set of
> daiy testreport starts).
> Cheers,
> Jean-Michel
> 
> On Fri, Jul 04, 2014 at 03:21:58PM -0400, Jean-Michel Campin wrote:
>> Hi Martin,
>> 
>> Now that the game is over (and the team that plays better won),
>> seems that not only lab_sea (standard test) but also 
>> global_ocean.cs32x15.seaice AD output have changed.
>> I assume it was expected.
>> 
>> I did run an AD testreport on baudelaire so I can check-in
>> the 2 new output_adm.txt files if you prefer.
>> Going to check what the TLM is doing.
>> 
>> Cheers,
>> Jean-Michel
>> 
>> On Fri, Jul 04, 2014 at 07:56:31AM -0400, Martin Losch wrote:
>>> Update of /u/gcmpack/MITgcm/pkg/seaice
>>> In directory forge:/tmp/cvs-serv1200/pkg/seaice
>>> 
>>> Modified Files:
>>> 	seaice_lsr.F 
>>> Log Message:
>>> my usual pre-weekend checkins that will destroy all relevant
>>> verification experiments:
>>>  - allays initialize local arrays in seaice_lsr_tridiagu/v
>>>  - add more local tapes to seaice_lsr_tridiagu/v
>>>  - move storing u/vice(:,:,bi,bj) out of seaice_lsr_tridiagu/v, but 
>>>    keep it in bi/bj loop; this seems to be essential
>>>  - gradients are still wrong but only by 4-7 orders of magnitude
>>>    instead of 15 orders of magnitude, lab_sea improves a little
>>> 
>>> update of lab_sea is comming up
>>> 
>>> 
>>> _______________________________________________
>>> MITgcm-cvs mailing list
>>> MITgcm-cvs at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-cvs
>> 
>> _______________________________________________
>> 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