[gradsusr] What's wrong in the script?
fadzil mnor
fadzilmnor84 at gmail.com
Mon May 30 01:40:44 EDT 2011
Line 20.
Does removing asterisk will solve the prob?
regards,
Fadzil
On 5/30/11, hassan quamrul <smquamrul77 at yahoo.com> 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
More information about the gradsusr
mailing list