[gradsusr] Grads pgm for the computation of stratospheric total heat flux

[resmi ea]

Computation of stratospheric heat flux kms-1 using NCEP/NCAR reanalysis
data for a single level

'reinit'
'sdfopen F:\windows\ncep\ncep2\datanceplevel\air.2008.nc'
'sdfopen F:\windows\ncep\uwnd\vwnd.2008.nc'
'enable print 4575N_08_10h0360.gmf'
'set lat 45 75'
'set lon 0 360'
'set lev 10'
'set mproj scaled'
'set parea 1.5 10 1.5 7'
'set grads off'
'set grid off'
'set ccolor 1'
'set cthick 8'
'set cstyle 1'
'set cmark 0'
'set strsiz 0.2'
'draw string 3.3 5.8 10hPa'
*'set vrange -100 500 50'
'set xlopts 1 4 0.20'
'set ylopts 1 4 0.20'
'set time 3feb2008 14mar2008'

* zonal average

'define at=ave(air.1,lon=0,lon=360)'
'define av=ave(vwnd.2,lon=0,lon=360)'
*departure from the zonal average

'define ate=air.1-at'
'define avw=vwnd.2-av'
'define az=ate*avw'

*zonalmean overbar [v*T*]
'set lon 60'
'define as=mean(az,lon=0,lon=360)'


* latitude weighting
'define a1=cos(lat*0.01744)'
'define a2=mean(a1,lat=45,lat=75)'
'define aa=as*a2'
'define ac=mean(aa,lat=45,lat=75)/mean(a2,lat=45,lat=75)'
'd ac'
'draw ylab K m/s'
'print'
'disable print'
'!gxeps -c -i 4575N_08_10h0360.gmf -o 4575N_08_10h0360.eps'

Regrads
-- 
Resmi E A
Research Scholar
Cochin university of science and technology
Cochin-16
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120508/221d311f/attachment-0003.html