- 积分
- 560
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
GRADS做出来的降水量图真的好难看,我不知道是个例转出来的数据有问题还是色标一类的设置有问题,麻烦各位大神帮我看看有哪里可以改进的吗?感激不尽!!!
以下是rain.gs
'reinit'
'open /home/mode/ARWpost/fh1.ctl'
j=0
k=23
i = 1
while(i<=6)
'set t 'i''
t=subwrd(result,4)
tc=substr(t,6,7)
'set gxout barb'
'set grid off'
'set grads off'
'set lat 10 30'
'set lon 110 135'
'set mpdset hires'
'set poli on'
'set display color white'
'set gxout shaded'
'set csmooth on'
'set clevs 9.9 24.9 49.9 99.9 249.9 299.9 349.9'
'set ccols 0 2 3 4 5 6 7 '
'define rain=rainnc+rainc'
'd rain'
'/home/mode/ARWpost/cbarn.gs'
*设置标题,标注等
'set string 1 l 2 0' ;*设置字符串颜色 位置 粗细 旋转角度
'set strsiz 0.11 0.15' ;*设置字符串的水平大小和高度大小
'draw string 1 0.3 Date: 'tc''
*'draw string 9 0.3 wt' ;*设置脚下标
*矢量图输出风场
*'set arrowhead 0.2'
*'set arrscl 0.5'
*'d skip(u,3,5);skip(v,3,5)'
*流线图输出风场
*'set gxout stream'
*'set line 0 1 6'
*'set lev 1000'
*'d u10;v10'
if(j=0)
'draw title 0807'k'00 Typhoon precipitation'
endif
if(j=6)
'draw title 0807'k'06 Typhoon precipitation'
endif
if(j=12)
'draw title 0807'k'12 Typhoon precipitation'
endif
if(j=18)
'draw title 0807'k'18 Typhoon precipitation'
endif
if(j=24)
k=k+1
j=0
'draw title 0807'k'00 Typhoon precipitation'
endif
j=j+6
'printim /home/mode/rain-'i'.png white'
'c'
i=i+1
Endwhile
End
|
|