[gradsusr] What's wrong in the script?

Jennifer Adams jma at cola.iges.org
Mon May 30 08:46:30 EDT 2011


Try changing statements like this:
    if (rainmax>10 & <=24)
to look like this:
    if (rainmax>10 & rainmax<=24)

--Jennifer

On May 30, 2011, at 1:30 AM, hassan quamrul wrote:

> Dear All,
>
> I am trying to plot rainfall from WRF output for different range  
> (say 10 mm to 250mm) for a perticular color cobination  
> colorbar .Please see below the script, but it not working (if  
> statement).
> What's wrong. Any can help? Please.
>
> 'reinit'
> 'open 2011052900.ctl'
> 'set display color white'
> 'c'
> 'set gxout shaded'
> 'set mpdset hires'
> 'set lat 18 28'
> 'set lon 84 97'
> 'set t 9'
> 'define rainmax=max(max(rainc,lon=84,lon=97),lat=18,lat=28)'
> 'd rainmax'
> rainmax=subwrd(result,10)
> if( rainmax=<10 )
> 'set clevs .1 1 2 4 6 8 10'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if(rainmax>10 & <=24)
> 'set clevs 1 4 8 12 16 20 24'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> *endif
> if(rainmax>24 & <=30)
> 'set clevs 1 5 10 15 20 25 30'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>30 & <=40 )
> 'set clevs 1 5 10 15 20 30 40'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>50 & <=60 )
> 'set clevs 1 5 10 20 30 40 60'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>60 & <=100 )
> 'set clevs 1 10 20 40 60 80 100'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>100 & <=140 )
> 'set clevs 1 10 20 40 70 100 140'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>140 & <=180 )
> 'set clevs 1 10 20 40 80 120 180'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> if( rainmax>180 & <=250 )
> 'set clevs 1 10 20 40 80 160 250'
> 'set ccols 0 3 10 7 12 8 2 6 6'
> endif
> 'd rainmax'
> 'cbar'
> 'draw title Rainfall(mm), 00Z,19.08.10 to 00Z,20.08.10'
> 'printim wrf_rain_2010081900_24.gif'
>
>
> regards.
>
> S.M.Quamrul Hassan
> Meteorologist
> Storm Warning Centre
> Bangladesh Meteorological Department
> Agargaon, Dhaka-1207
> Tel: +880 2 9135742
> Fax: +880 2 8118230
> Cell: +8801916255449
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110530/4197d1c1/attachment-0003.html 


More information about the gradsusr mailing list