[MITgcm-devel] seaice ridging: adjoint

[Martin Losch]

OK, done, thanks,

The adjoint works but gives slightly different numbers. for SEAICE_ITD, there is no adjoint, as far as I can see.

Martin

On Apr 23, 2014, at 12:38 PM, Patrick Heimbach <heimbach at MIT.EDU> wrote:

> Hi Martin,
> go for it. I’m around, so can take a look.
> p.
> 
> On Apr 23, 2014, at 5:04 AM, Martin Losch <Martin.Losch at awi.de> wrote:
> 
>> Hi Patrick,
>> 
>> in preparation for real seaice ridging code, I am moving the regularization part (most of “PART I”) out of seaice_growth into a new routine.
>> 
>> Unfortunately this changes the adjoint gradient values (making the agreement with the FD gradient a little worse). All forward runs are OK:
>> G D M    C  A  F
>> e p a R  o  d  D
>> n n k u  s  G  G
>> 2 d e n  t  r  r
>> 
>> Y Y Y Y 16> 8<16 FAIL  global_ocean.cs32x15.seaice
>> Y Y Y Y 16> 7<16 FAIL  global_ocean.cs32x15.seaice_dynmix
>> Y Y Y Y 16> 5<16 FAIL  lab_sea  (e=0, w=18)
>> Y Y Y Y 16> 6<16 FAIL  lab_sea.noseaicedyn
>> 
>> The taf_ad.log is the same with and without my modifications (let alone line numbers and the extra routine), so no extra. I have no idea where to look. Can I check in my changes anyway and we have a look together?
>> 
>> Martin
>> _______________________________________________
>> 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-1420 | 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