[MITgcm-devel] more seaice blues for the adjoint

[Martin Losch]

Patrick
I have checked in my solution, with #ifdefs. If you feel that the  
ifdef's are too much, then I can remove them and have the fortran IF  
for all case (I just thought it makes the code slower).

Martin

On 6 Jun 2007, at 16:13, Patrick Heimbach wrote:

>
> 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
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel