[MITgcm-support] a problem with EXF

Matthew Mazloff mmazloff at MIT.EDU
Wed Sep 7 16:24:10 EDT 2005


Hi Ming,

I am a bit surprised; I just looked at the EXF & MNC package and I don't 
the model is able to read in NetCDF format input files.  Maybe someone 
in the know can tell us what the status is of using NetCDF input. 

Anyway, I know the model can read in binary input files and the 
precision is given by exf_iprec = 32 (or 64).  I actually sent an email 
under the thread name "Starting with MITgcm" that was about generating 
input files a few weeks ago which I'll append below.  You should easily 
be able to tell if the forcing fields are being read in correctly.   
Check the monitor statistics for your forcing fields and see if they 
seem reasonable.  Also, if the model crashes because it cannot find an 
input file, it should give an error message that includes the name of 
the file it last looked for.  Remember that if you give the input file 
name exf_input.data and the model can't find this file it may append a 
".data" on it and try again.  So if the error says it can't find 
exf_input.data.data that's why. 

And just in case you were referring to the previous issue of array 
precision, I'll add that exf_yftype refers to the precision with which 
the exf package should read arrays.  Your choices are RL for REAL*8 or 
RS for REAL*4. 

once again hoping I'm giving accurate info,
Matt


Here's the old email:

Hi Matej,

As Ed said, the verification experiments have matlab scripts that 
describe preperation of input files. For example, in MITgcm / 
verification / exp5 / input / gendata.m
a surface forcing file is made:

Q=Qo*(1+0.01*rand([nx,ny]));

it is size [nx,ny]
if it were 3-d (e.g. and initial theta file) it would be size [nx,ny,nz]

then it is written with

fid=fopen('Qnet.circle','w','b');
fwrite(fid,Q,'real*8');
fclose(fid);

The data is written by fwrite down the matrix columns.

it is written translated into real*8 because the mitgcm data file given 
in the verification experiment has
readBinaryPrec=64,

it is read in as a surface heating file by putting
# Input datasets
&PARM05
surfQfile='Qnet.circle',
&
in the data file
other possible initial files you can give the model are (as seen in  
MITgcm / model / src / ini_parms.F )
C--   Input files
     NAMELIST /PARM05/
    & bathyFile, topoFile,
    & hydrogThetaFile, hydrogSaltFile,
    & zonalWindFile, meridWindFile,
    & thetaClimFile, saltClimFile,
    & surfQfile, surfQnetFile, surfQswFile, EmPmRfile, saltFluxFile,
    & lambdaThetaFile, lambdaSaltFile,
    & uVelInitFile, vVelInitFile, pSurfInitFile,
    & dQdTFile, ploadFile,tCylIn,tCylOut,
    & eddyTauxFile, eddyTauyFile,
    & mdsioLocalDir,
    & the_run_name
CEOP



-matt



Ming Guo wrote:

> Hi, Matt,
>
> just wonder the format of the focing files, such as ustress and 
> vstress, should be NetCDF file, or something else?
>
> thanks,
> Ming
> **********************************
> Ming Guo,
> Physics and Physical Oceanography Department,
> Memorial University of Newfoundland,
> Office Phone:709-737-2407
> E-mail:gming at karluk.physics.mun.ca
> ***********************************
>
>
> On Mon, 5 Sep 2005, Matthew Mazloff wrote:
>
>> Hi Ming,
>>
>> The problem is
>>
>> yftype_exf='RL',
>>
>> in data.exf
>>
>> I looked in MITgcm/pkg /exf/exf_readparms.F
>> and found
>> exf_yftype but not yftype_exf.
>>
>> You should switch this.
>>
>> -Matt
>>
>>
>>
>>
>>
>> Ming Guo wrote:
>>
>>> Hi, Ed,
>>> I just set the files, but the errors still there,
>>>
>>> "
>>> namelist read: variable not in namelist
>>> apparent state: unit 11 named /tmp/tmp.FA6VrWo
>>> last formate: list io
>>> latey reading sequential formatted external IO
>>> Abort (core dumped)
>>> "
>>>
>>> The attachments are my "EXF" files.
>>> Can you help me to find out the problem?
>>> By the way, is there any spectial format needed for the force file?
>>>
>>> thanks,
>>> Ming
>>>
>>>
>>>
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20050907/4637fa00/attachment.htm>


More information about the MITgcm-support mailing list