[MITgcm-devel] seaice ridging: adjoint

[Martin Losch]

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