[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