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

[James T. Potemra]

Geraldo,

Do the plots look reasonable?  A message like this is sometimes caused 
by a corrupt data point (or points), e.g., a missing value is not set 
right and appears as some large number.

Jim

On 1/24/15 1:41 PM, Geraldo Deniro wrote:
> 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.gs
> No hardcopy metafile open
> All files closed; all defined objects released;
> All GrADS attributes have been reinitialized
> Too many contour levels -- adjusting cint
> ga->
>
> The script:
>
> ------------------------------
> 'reinit'
>
> yr=1981
> while(yr<=2010)
> 'sdfopen slp.'yr'.nc'
> yr=yr+1
> endwhile
>
> 'sslp1=0.0'
> 'sslp2=0.0'
>
> yr=1981
> ff=1
> while(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+1
> yr=yr+1
> endwhile
>
> '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
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150126/12967509/attachment.html