[gradsusr] Fortnightly anomaly of reanalysis 1 (R1) daily data

Geraldo Deniro geraldodeniro at hotmail.com
Sat Jan 24 18:41:30 EST 2015


Hi GrADS Users,I made a script to make fortnightly anomaly of the mean sea level pressure daily data of R1, when i play this scripts on GrADS a message with error appears. Does anyone knows how to solve it or a easy way to make this script?
The error:
ga-> slp.gsNo hardcopy metafile openAll files closed; all defined objects released;All GrADS attributes have been reinitializedToo many contour levels -- adjusting cintga->
The script:
------------------------------'reinit'
yr=1981while(yr<=2010) 'sdfopen slp.'yr'.nc'yr=yr+1endwhile
'sslp1=0.0''sslp2=0.0'
yr=1981ff=1while(yr<=2010) 'set dfile 'ff'set time 01nov'yr'slp1=ave(slp,time=01Sep'yr',time=15Sep'yr')''slp2=ave(slp,time=16Sep'yr',time=30Sep'yr')''sslp1=sslp1+(slp1/30)''sslp2=sslp2+(slp2/30)'ff=ff+1yr=yr+1endwhile
'sdfopen slp.2014.nc''set dfile 'ff'set time 01nov2014''slp141=ave(slp,time=01nov2014,time=15nov2014)' 'slp142=ave(slp,time=16nov2014,time=30nov2014)' 
'rgbset''set map 1 1 10''set mpdset brmap_hires''set lon 240 360''set lat -80 20'
'set vpage 0.0 11.0 0.0 8.5''set parea 0.5 10.2 1.0 7.5''set gxout shaded''set grads off''set xlopts 1 0.18 0.18''set ylopts 1 0.18 0.18''set xlint 20''set clevs -12 -9 -6 -3 3 6 9 12''set ccols 45 44 43 41 0 21 23 24 25''d (slp141-sslp1)/100''cbarm''set gxout contour''set grads off''set cthick 6''set clab forced''set cint 50''set digsize 0.18''d slp141''draw title NCEP/NCAR R1 \ PNMM & Anomalia de PNMM \ 01-15 novembro/2014''set display color white''printim slp_0115nov2014.jpg'
'q pos''c''set vpage 0.0 11.0 0.0 8.5''set parea 0.5 10.2 1.0 7.5''set gxout shaded''set grads off''set xlopts 1 0.18 0.18''set ylopts 1 0.18 0.18''set xlint 20''set clevs -12 -9 -6 -3 3 6 9 12''set ccols 45 44 43 41 0 21 23 24 25''d (slp142-sslp2)/100''cbarm''set gxout contour''set grads off''set cthick 6''set clab forced''set cint 50''set digsize 0.18''d slp142''draw title NCEP/NCAR R1 \  PNMM & Anomalia de PNMM \ 16-30 novembro/2014''set display color white''printim slp_1630nov2014.jpg'
return


Geraldo Deniro geraldodeniro at ufrj.br Graduate Student of Meteorology Department of Meteorology Institute of Geosciences Federal University of Rio de Janeiro 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150124/77601877/attachment-0001.html 


More information about the gradsusr mailing list