[MITgcm-support] Re : Query regarding "exf package".
Martin Losch
Martin.Losch at awi.de
Mon Jan 9 03:57:28 EST 2012
Himansu,
for constant wind stress (one time slice) you need to have u/vtressperiod = 0. (which is the default if you do not specify anything) in data.exf, otherwise the model is trying to read the next field at t=t0+u/vstressperiod.
Further, readstressOnAgrid=.TRUE., makes the model believe that the stress components are defined on the center points of the C-grid (I think), so that they are interpolated to the U/V points. Specifying
zonalWindFile='pUstress.bin',
meridWindFile='pVstress.bin',
in data&PARM05 assumes stress at U/V points automatically. Set readStressOnCgrid = .true., if you want to do the same in with pkg/exf.
Martin
On Jan 9, 2012, at 7:57 AM, himansu pradhan wrote:
> hi MITgcm users,
>
> A doubt to clarify
>
> My u/v-windstress file is of one time slice (157 * 133) . This is what i was using in "data" file .
> *******data******
> &PARM03
> startTime=0,
> endTime=864000,
> deltaT=10.0,
> abEps=0.1,
> # momDissip_In_AB=.FALSE.,
> pChkptFreq=432000,
> # chkptFreq=432000,
> dumpFreq=21600,
> &
> # Input datasets
> &PARM05
> bathyFile='top_MOD.bin',
> hydrogThetaFile='15lvlT_may.bin',
> hydrogSaltFile='15lvlS_may.bin',
> # zonalWindFile='pUstress.bin',
> # meridWindFile='pVstress.bin',
> &
> but when i use the SAME u/v- windstress file (157*133) in "data.exf" . The model didnot run ..!!!!!!!!
> # *********data.exf************
> # External Forcing Data
> # *********************
> &EXF_NML_01
> exf_iprec = 64,
> readstressOnAgrid=.TRUE.,
> useExfCheckRange=.FALSE.,
> &
> &EXF_NML_02
> ustressstartdate1=20110101,
> ustressstartdate2=000000,
> ustressperiod=864000,
> #
> vstressstartdate1=20110101,
> vstressstartdate2=000000,
> vstressperiod=864000,
> #
> ustressfile = 'pUstress.bin',
> vstressfile = 'pVstress.bin',
> &
> [....]
>
> finally,
> i created u/v- windstress file having 2 time slabs (i.e. 157*133*2)............ The model "runs" SUCESSFULLY for the total period of intergration. .
>
> my doubt is :
>
> (a) why i need to give 2 time slabs (157*133*2) in "data.exf".....??? (whereas, In "data" i use only one time slab(157*133) for the complete period of integration)
>
> Am i doing right here ???
>
>
> Thank you ,
>
> Himansu Kesari Pradhan,
> India.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list