[MITgcm-support] obcs_prescribe and *.h conflict?
m schaferkotter
schaferk at bellsouth.net
Wed Mar 30 08:41:30 EST 2005
patrick:
i did find the culprit.
> 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?
in data.exf
# *********************
# External Forcing Data
# *********************
&EXF_NML
#
repeatperiod = 172800.,
obcsWstartdate1 = 19800101,
obcsWstartdate2 = 000000,
obcsWperiod = 600.0,
obcsEstartdate1 = 19800101,
obcsEstartdate2 = 000000,
obcsEperiod = 600.0,
exf_iprec = 64, <---- THIS
was '32'.
exf_yftype = 'RL',
&
the simulation now runs for a short while, and i have to make some
adustments.
michael
On Mar 25, 2005, at 11:39 AM, Patrick Heimbach wrote:
>
> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>
More information about the MITgcm-support
mailing list