Bug in quadratic ice shelf cavity stress

[bgalton]

Hi all,

There is a bug in the calculation of the u-component quadratic ice shelf cavity stress (sustr).

Line 353 in set_vbc.F:

cff2=SQRT(u(i,j,N(ng),nrhs)*u(i,j,N(ng),nrhs)+cff2*cff2)

should read:

cff2=SQRT(u(i,j,N(ng),nrhs)*u(i,j,N(ng),nrhs)+cff1*cff1)

Cheers,

Ben Galton-Fenzi

[arango]

Yes, thank you. I corrected this bug in the adjoint codes but I forgot to correct it in the nonlinear model. I updated the repository, check the following track ticket.

Thank you for reporting this bug