[gradsusr] regarding VIMF

[Sanjukta Rani Padhi]

Dear GrADS members,

I have a doubt regarding the multiplication of the 1/g and vector field. I
wrote two types of script... I don't know which is correct.

my 1st script***************************************

'reinit'
'set display color white'
'clear'
'sdfopen shum.mon.mean.nc'
'sdfopen uwnd.mon.mean.nc'
'sdfopen vwnd.mon.mean.nc'

'set lat -30 30'
'set lon 30 120'

'set lev 1000 300'

'define psurf=uwnd.2*0.0+1000.0'


'define qu = shum.1*uwnd.2'
'define qv = shum.1*vwnd.3'
'set lev 1000'
'define uvt = vint(psurf,qu,300)'
'define vvt = vint(psurf,qv,300)'
'define ivt = sqrt(uvt*uvt+vvt*vvt)'

'set xlint 10'
'set ylint 10'
'set map 1 1 8'
'set grads off'

'set gxout shaded'
'd IVT'

'cbarn'
'set gxout vector'
*'d uvt;vvt'*

'gxprint IVT.pdf'

*****************************************************************

my 2nd script**************************

'reinit'
'set display color white'
'clear'
'sdfopen shum.mon.mean.nc'
'sdfopen uwnd.mon.mean.nc'
'sdfopen vwnd.mon.mean.nc'

'set lat -30 30'
'set lon 30 120'

'set lev 1000 300'

'define psurf=uwnd.2*0.0+1000.0'
'define g=9.80665'

'define qu = shum.1*uwnd.2'
'define qv = shum.1*vwnd.3'
'set lev 1000'
'define uvt = vint(psurf,qu,300)'
'define vvt = vint(psurf,qv,300)'

*'define uvt1=1/g*uvt''define vvt1=1/g*vvt'*
'define ivt = sqrt(uvt1*uvt1+vvt1*vvt1)'

'set xlint 10'
'set ylint 10'
'set map 1 1 8'
'set grads off'

'set gxout shaded'
'd IVT'

'cbarn'
'set gxout vector'
*'d u;v'*

'gxprint IVT.pdf'

*****************************************************************
Kindly help me with this topic.

Thank you and regards

Sanjukta Rani Padhi
Senior Research Fellow(DST INSPIRE)
Dept. of Marine Sciences
Berhampur University
Bhanja Bihar-760007
Berhampur, Odisha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20200902/31853f11/attachment.html>