[MITgcm-support] obcs_prescribe and *.h conflict?
Michael Schaferkotter (Contractor)
Michael.Schaferkotter at nrlssc.navy.mil
Fri Mar 25 11:09:08 EST 2005
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
(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',
any pointers are appreciated.
More information about the MITgcm-support