[MITgcm-devel] more seaice blues for the adjoint

[Jinlun Zhang]

Hi Martin,
Would this be due to some kind of finite differencing problem?
Jinlun

Martin Losch wrote:
> 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