[MITgcm-support] Re : Query regarding "exf package".

Martin Losch Martin.Losch at awi.de
Tue Jan 10 09:13:06 EST 2012


Himansu,

1. please do not send large files to this list, rather post them somewhere and send the link.
2. I cannot read *.fig files
3. Your jpg-plots look like the surface forcing you specify does not make sense. Check (in both experiments) that the model sees the correct forcing, e.g. by writing the appropriate diagnostics with the diagnostics-package (oceTAUX and oceTAUY).

Please send your entire data and data.exf and data.cal files, as you use them in your experiments. From you previous emails it was not clear to me, what you are actually specifying. 

M.

On Jan 10, 2012, at 8:18 AM, himansu pradhan wrote:

> Hi Matrin Sir,
> 
> Thanks for the reply....... 
> your suggestation worked.... (both : setting the timeperiod and stress on C grid)
> 
> # ********data.exf*************
> # External Forcing Data
> # *********************
>  &EXF_NML_01
>  exf_iprec         = 64,
> #exf_verbose       = .TRUE.,
>  readstressOnCgrid=.TRUE.,
>  useExfCheckRange=.FALSE.,
>  &
>  &EXF_NML_02
>  ustressstartdate1=20110101,
>  ustressstartdate2=000000,
>  ustressperiod=0,
> #
>  vstressstartdate1=20110101,
>  vstressstartdate2=000000,
>  vstressperiod=0,
> #
>  ustressfile  = 'May23Ustr.bin',
>  vstressfile  = 'May23Vstr.bin',
> [.......]
> 
> ........but........ ?????
> the result (plot) for this run (i.e. for "current","temperature" , salinity ) are not fine...?????.... the figures attached can give a clear look at that.........
> 
> The results of the model run where u/v- windstress files are used in "data.exf" are in  Figure (a) and (b). The plot for vectors seems more distorted ( donot know why...???????........). 
> 
> (a) Current_dataexf.fig 
> (b) temp_dataexf.jpg   
> 
> The results of the model run where windstress files are used in "data" are Figure (c) and (d). This results are FINE.....(This is what i wish to get whn using u/v-windstress files in "data.exf")  
> (c)   Current_data.fig   
> (d) temp_data.jpg
> 
> 
> Any suggestations 
> 
> Thank you,
> 
> Himansu Kesari Pradhan
> India
> 
> 
> On Mon, Jan 9, 2012 at 2:27 PM, Martin Losch <Martin.Losch at awi.de> wrote:
> 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
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 
> <Current_dataexf.fig><Current_data.fig><temp_dataexf.jpg><temp_data.jpg>_______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list