[MITgcm-support] a problem with EXF
Ming Guo
gming at karluk.physics.mun.ca
Thu Sep 15 12:33:04 EDT 2005
Hi,Matt,
I am still not sure what you mean ustre(x,y,t1),ustre(x,y,t2)....
For example, for a 2*2 grid field,
is the ustressfile like: ustre(1,1,t1),ustre(1,1,t2), ....,ustre(2,2,t1)
ustre(2,2,t2),....
or the ustre(x,y,t1),ustre(x,y,t2)... are just a series of stress
matrixes and each denots one time step?
I tried the latter method, but the value it reads seems not right.
Thanks,
Ming
On Wed, 14 Sep 2005, Matthew Mazloff wrote:
> Hi Ming,
>
> The model does need a ustress value at every grid point and at every
> timepoint. So the file must be like ustre(x,y,t1),ustre(x,y,t2)... There is
> an option to interpolate the exf files so you can give a coarse ustress field
> and have the model interpolate to your grid at each timestep, this would
> allow you to use vey small input files. Currently this is only setup to work
> on a lat, long grid.
> You can check out
> [MITgcm] / MITgcm_contrib / quarter_degree_global / input / data.exf
> maybe this link will work:
> http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/quarter_degree_global/input/data.exf?rev=HEAD&content-type=text/vnd.viewcvs-markup
> to see a sample of the data.exf for the interpolation routine.
>
> You'll also need to define USE_EXF_INTERPOLATION, in ECCO_CPPOPTIONS
>
> Matt
>
>
> Ming Guo wrote:
>
>> 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
>>>
>>>
>>>
>> _______________________________________________
>> 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