[MITgcm-devel] pkg/obcs land mask

[Dimitris Menemenlis]

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:

>
> Hi,
>
> printing out SST and surface SALINITY was a revelation to me...
>
> You can see the location in the full arctic domain in iglobal and  
> jglobal indices (thanks PAtrick again....)
>
> FLUX CO2 (mol kg-1 m s-1) =  NaN
> POTENTIAL SURF TEMP =   0.000000000000000E+000
> SURFACE SALINITY  =  0.000000000000000E+000
> SURF ALK  =   0.000000000000000E+000
> SURF DIC =   0.000000000000000E+000
> pCO2 Ocean =  NaN
> iglobal =          420
> jglobal =           86
>
> This case above should be the Pacific boundary, right ?
>
> Other cases only DIC goes to NaN although ocean physcal parameters
> are fine beacuse I guess the NaN get propagated in tiem stepping the  
> DIC tracer :
>
> FLUX CO2 (mol kg-1 m s-1) =  NaN
> POTENTIAL SURF TEMP =    7.33224004059763
> SURFACE SALINITY  =   32.9180385715927
> SURF ALK  =    2.10001697524426
> SURF DIC =  NaN
> pCO2 Ocean =  NaN
> iglobal =          418
> jglobal =           84
>
>
> Comments are welcome, as usual...
>
> Manfredi
>
>
> -- 
>
>
> Manfredi Manizza
> Room 54-1410
> EAPS-MIT
> 77 Massachusetts Avenue
> Cambridge, MA 02139-4307, USA
> phone: 617-253-3573
> fax  : 617-253-4464
> email: mmanizza at ocean.mit.edu
> web  : http://web.mit.edu/mmanizza/www/