<br clear="all">Computation of stratospheric heat flux kms-1 using NCEP/NCAR reanalysis data for a single level<br><br>&#39;reinit&#39;<br>&#39;sdfopen F:\windows\ncep\ncep2\datanceplevel\<a href="http://air.2008.nc">air.2008.nc</a>&#39;<br>
&#39;sdfopen F:\windows\ncep\uwnd\<a href="http://vwnd.2008.nc">vwnd.2008.nc</a>&#39;<br>&#39;enable print 4575N_08_10h0360.gmf&#39;<br>&#39;set lat 45 75&#39;<br>&#39;set lon 0 360&#39;<br>&#39;set lev 10&#39;<br>&#39;set mproj scaled&#39;<br>
&#39;set parea 1.5 10 1.5 7&#39;<br>&#39;set grads off&#39;<br>&#39;set grid off&#39;<br>&#39;set ccolor 1&#39;<br>&#39;set cthick 8&#39;<br>&#39;set cstyle 1&#39;<br>&#39;set cmark 0&#39;<br>&#39;set strsiz 0.2&#39;<br>&#39;draw string 3.3 5.8 10hPa&#39;<br>
*&#39;set vrange -100 500 50&#39;<br>&#39;set xlopts 1 4 0.20&#39;<br>&#39;set ylopts 1 4 0.20&#39;<br>&#39;set time 3feb2008 14mar2008&#39;<br><br>* zonal average<br><br>&#39;define at=ave(air.1,lon=0,lon=360)&#39;<br>&#39;define av=ave(vwnd.2,lon=0,lon=360)&#39;<br>
*departure from the zonal average<br><br>&#39;define ate=air.1-at&#39;<br>&#39;define avw=vwnd.2-av&#39;<br>&#39;define az=ate*avw&#39;<br><br>*zonalmean overbar [v*T*]<br>&#39;set lon 60&#39;<br>&#39;define as=mean(az,lon=0,lon=360)&#39;<br>
<br><br>* latitude weighting<br>&#39;define a1=cos(lat*0.01744)&#39;<br>&#39;define a2=mean(a1,lat=45,lat=75)&#39;<br>&#39;define aa=as*a2&#39;<br>&#39;define ac=mean(aa,lat=45,lat=75)/mean(a2,lat=45,lat=75)&#39; <br>&#39;d ac&#39;<br>
&#39;draw ylab K m/s&#39;<br>&#39;print&#39;<br>&#39;disable print&#39;<br>&#39;!gxeps -c -i 4575N_08_10h0360.gmf -o 4575N_08_10h0360.eps&#39;<br><br>Regrads<br>-- <br>Resmi E A<br>Research Scholar<br>Cochin university of science and technology<br>
Cochin-16<br><br>