[MITgcm-devel] more seaice blues for the adjoint

[Patrick Heimbach]

Hi Martin,

not sure right now.
I'd rather avoid having them between #ifdef's,
but maybe it's unavoidable?
You could at least start by avoiding evaluation of sqrt()
for deltaC == 0.
Then worry about small denominators separately.

-Patrick



On Jun 6, 2007, at 8:53 AM, 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
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach