[MITgcm-support] obcs_prescribe and *.h conflict?

Patrick Heimbach heimbach at MIT.EDU
Fri Mar 25 12:39:32 EST 2005


Michael,

there seem to have problems unrelated to open boundaries:

1) Your input theta, salt:
are you sure they're OK.
You have initial values of E+32.
That's typically the case if you have little endian files
and read them as big endian.
Could that be it?

2) Forcing:
what are you forcing with, or atmospheric state and bulk formulae?
Also, you seem to have set relaxation, i.e.
ALLOW_CLIMSST_RELAXATION
which is why the model wants a data.exf_clim.
Do you actually want relaxation?

-Patrick



Quoting Michael Schaferkotter <Michael.Schaferkotter at nrlssc.navy.mil>:

> ed:
> 
> following your lead, i was able to build.
> 
> i added exf,cal to code/packages.conf (as heimbach Sep22,2004 MITgcm-support
> Re: OBCS sponge and prescribe code ...
> i added useOBCprescribe=.TRUE. to data.obcs.
> 
> the first run died with:
> 
> PGFIO-F-209/namelist read/unit=9/'OLD' specified for file which does not
> exist.
>  In source file exf_clim_readparms.f, at line number 2279
> 
> the code indicated a needed namelist file data.exf_clim. found one and copied
> it to run directory.
> i set values in data.exf_clim for null effect (no climatology at this time).
> 
> with a time step of 1s and dumpFreq of 1s the code doesn/t even write the
> first output file:
> 
> MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF SOLUTION
> 
> STDERR.0000
> (PID.TID 0000.0001) SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=  -4.933E+32
>   5.849E+35
> (PID.TID 0000.0001) MON_SOLUTION: STOPPING CALCULATION
> 
> (PID.TID 0000.0001) %MON dynstat_theta_max            =  
> 5.8486031497624E+35
> (PID.TID 0000.0001) %MON dynstat_theta_min            = 
> -4.9326819537111E+32
> (PID.TID 0000.0001) %MON dynstat_theta_mean           =  
> 4.7703465478160E+32
> (PID.TID 0000.0001) %MON dynstat_theta_sd             =  
> 1.6703490233395E+34
> (PID.TID 0000.0001) %MON dynstat_theta_del2           =  
> 7.1675998670440E+32
> (PID.TID 0000.0001) %MON dynstat_salt_max             =  
> 3.5000000000000E+01
> (PID.TID 0000.0001) %MON dynstat_salt_min             =  
> 0.0000000000000E+00
> (PID.TID 0000.0001) %MON dynstat_salt_mean            =  
> 3.3660468750000E+01
> (PID.TID 0000.0001) %MON dynstat_salt_sd              =  
> 6.5692973000645E+00
> (PID.TID 0000.0001) %MON dynstat_salt_del2            =  
> 6.6976562500000E-01
> 
> i have a small 50x50x50 domain d[XYZ]=2.0. the salinity is constant 35.0 and
> the temperatue is
> a simple linear stratification 10.0 at top and 9.775 at bottom.
> (dTdz=0.005).
> 
> for the prescribed values i supply OBW[stuv].bin and OBE[stuv].bin.
> 
> u velocity is prescribed as a sine curve with period 2 days and amplitude
> 23cm/s.  there are 289 records at 10 minute intervals. (on East and West)
> u,v are zero elsewhere. 
> T is the same as (original distribution) and S is constant.
> 
> i/m not sure if i hosed something by commenting out most of data.exf_clim...
> only 
>  climtempfile       = '',
>  climsaltfile       = '',
>  climsstfile        = '',
>  climsssfile        = '',
> #
>  exf_clim_iprec     = 32,
>  exf_clim_yftype    = 'RL',
>  &
> 
> is uncommented.
> 
> 
> any pointers are appreciated.
> 
> michael
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> 


--------------------------------------------------------
Patrick Heimbach   Massachusetts Institute of Technology
FON: +1/617/253-5259                  EAPS, Room 54-1518
FAX: +1/617/253-4464             77 Massachusetts Avenue
mailto:heimbach at mit.edu               Cambridge MA 02139
http://www.mit.edu/~heimbach/                        USA




More information about the MITgcm-support mailing list