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

himansu pradhan oceancalling at gmail.com
Tue Jan 10 02:18:36 EST 2012


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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120110/afa0312c/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Current_dataexf.fig
Type: application/octet-stream
Size: 1182055 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120110/afa0312c/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Current_data.fig
Type: application/octet-stream
Size: 1162164 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120110/afa0312c/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp_dataexf.jpg
Type: image/jpeg
Size: 349773 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120110/afa0312c/attachment-0002.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp_data.jpg
Type: image/jpeg
Size: 100063 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120110/afa0312c/attachment-0003.jpg>


More information about the MITgcm-support mailing list