Dear All,<br><br>I like to plot the olr and win vectors in the same plot. i am trying to plot, but only olr has been plotted. i want to show both in the same plot. please give me solution for them.<br><br>script:1<br><br>'reinit'<br>
'sdfopen h:\ncep\ncep\<a href="http://olr1980.nc">olr1980.nc</a>'<br>'sdfopen h:\ncep\ncep-winds\10m\<a href="http://uwnd.10m.gauss.2010.nc">uwnd.10m.gauss.2010.nc</a>'<br>'sdfopen h:\ncep\ncep-winds\10m\<a href="http://vwnd.10m.gauss.2010.nc">vwnd.10m.gauss.2010.nc</a>'<br>
'enable print h:\ncep\ncep\olrwin-80.gmf'<br>'set lat -35 25'<br>'set lon 265 345'<br>tt=1<br>while(tt<=12)<br>'set t 'tt<br>'define ol=ulwrf.1'<br>'set gxout contour'<br>
*'set gxout shaded'<br>'set grads off'<br>'set grid on'<br>'set ccolor -1'<br>'set clevs 180 190 200 210 220 230 240 250 260 270 280 290 300'<br>*'set ccols 180 190 200 210 220 230 240 250 260 270 280 290 300'<br>
'd ulwrf.1'<br>*'cbarn'<br>'q time'<br>a=sublin(result,1)<br>b=subwrd(a,3)<br>'draw title OLR 'b<br><br>'define u1=uwnd.2'<br>'define v1=vwnd.3'<br>'set gxout vector'<br>
'set grads off'<br>'set grid on'<br>'d u1;v1'<br>'print'<br>'c'<br>tt=tt+1<br>endwhile<br>'disable print'<br>'reinit'<br><br>result: getting only contours, no vectors.<br>
<br><br>another, i want to calculate the anomalies. i have both year file (with monthly mean) at the same long term monthly mean, i would like to calculate the anomaly and plot the same. monthly mean data was taken but the long term mean file showing error, file number out of range. <br>
<br>Script: 2<br><br>'reinit'<br>'sdfopen h:\ncep\ncep-olr-mon\<a href="http://1975olr.nc">1975olr.nc</a>'<br>'sdfopen h:\ncep\ncep-olr-mon\<a href="http://olr.mon.ltm.nc">olr.mon.ltm.nc</a>'<br>'enable print h:\ncep\ncep-olr-mon\1975olranom.gmf'<br>
*'set lat -35 25'<br>*'set lon 165 345'<br>tt=1<br>while(tt<=12)<br>'set t 'tt<br>'define ol=olr.1'<br>'define lt=olr.2'<br>'define anom=ol-lt'<br>'set lat -35 25'<br>
'set lon 165 345'<br>'set gxout contour'<br>*'set gxout shaded'<br>'set grads off'<br>'set grid on'<br>'set ccolor -1'<br>*'set clevs 180 190 200 210 220 230 240 250 260 270 280'<br>
'd anom'<br>*'cbarn'<br>'q time'<br>a=sublin(result,1)<br>b=subwrd(a,3)<br>'draw title OLR 'b<br>'print'<br>'c'<br>tt=tt+1<br>endwhile<br>'disable print'<br>'reinit'<br>
<br>result: error.<br><br>please give me solution for those two. i am using windows based grads (Grads 2.0.1.oga.1 in windows vista) <br clear="all"><br><br>thanks in advance.<br>-- <br>M.V.Subrahmanyam<br><br>