[MITgcm-support] a problem with EXF
Ming Guo
gming at karluk.physics.mun.ca
Wed Sep 14 15:29:54 EDT 2005
Hi, Matt,
As I am using uniform time-dependent wind stress, how to generate the
stress file? just create ustressfile like, ustre(t1), ustre(t2)... and the
model can give ustress value to all the grids, or I also need to specify
the value of every grids as gendata.m in exp0?
thanks
Ming
On Wed, 7 Sep 2005, Matthew Mazloff wrote:
> 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
>>
>>
>> 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
>
>
More information about the MITgcm-support
mailing list