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

Michael Schaferkotter (Contractor) Michael.Schaferkotter at nrlssc.navy.mil
Fri Mar 25 11:09:08 EST 2005


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



More information about the MITgcm-support mailing list