[MITgcm-devel] pkg/obcs land mask

[Dimitris Menemenlis]

Martin, just to let you know that your modifs to pkg/obcs work!   
Cheers, D.

> I guess it works now thanks :D but I am also in the process of  
> updating the DIC(K) :D code that David has modified, bloody hell, it  
> never ends !!!!
>
> Infact after your modification it was not giving me NAN anymore but  
> then it was not outputting the dic* files with CO2 fluxes..
>
> On May 20, 2008, at 10:54 AM, Dimitris Menemenlis wrote:
>
>>
>> Manfredi, location (420,86), which is an open boundary, is wet but  
>> location (419,86) is dry land.  What happens is that pkg/obcs will  
>> mask all locations that are adjacent to land immediately next to  
>> the open boundary, so it will set to zero values at (420,86), which  
>> is what causes your DIC code to blow up.
>>
>> This is a known problem and it is supposed to have been fixed by  
>> changes to pkg/obcs checked in by Martin Losch on April 24-25.   
>> Please make sure that you are using code more recent than April 25  
>> and that OBCSfixTopo = .true. in data.obcs.  This is the default  
>> behavior.
>>
>> I am pretty sure that new code should fix your problem but let me  
>> know if it does not.  Dimitris
>>
>> Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
>> Jet Propulsion Lab, California Institute of Technology
>> MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
>> tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720
>>
>> On May 18, 2008, at 11:20 AM, Manfredi Manizza wrote:

Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720