[MITgcm-devel] more seaice blues for the adjoint

[Martin Losch]

Hi Patrick,

I have found (or rather, was pointed to) a problem with the seaice  
solvers: The start to move spontaneously (in the absence of forcing),  
if the sea ice distribution is NOT uniform.
I have implemented a fix but this will cause problems with the  
adjoint: I need terms like
SQRT(deltaC), which used to be SQRT(MAX(deltaC,SEAICE_EPS_SQ)), so  
that the derivate code will be involve 1/sqrt(deltaC). Should I put  
this into #ifdefs?

Martin