[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