- 积分
- 787
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请各位大侠指点:原想把色标竖放的,确成这样了?!(见图)
以下是我的脚本:
'reinit'
'set rgb 29 205 0 255'
'set rgb 28 255 0 200'
'set rgb 27 255 0 0'
'set rgb 26 255 100 0'
'set rgb 25 255 160 0'
'set rgb 24 255 200 0'
'set rgb 23 255 255 0'
'set rgb 22 200 255 0'
'set rgb 21 140 255 0'
'set rgb 20 180 250 170'
*****************************blue
'set rgb 41 30 255 255'
'set rgb 42 0 200 255'
'set rgb 43 0 145 255'
'set rgb 44 0 80 255'
'set rgb 45 0 0 255'
'set rgb 46 85 0 255'
'set rgb 47 130 0 255'
'set rgb 48 180 0 255'
*********************************
'open h:\fnl\11\2017062008.ctl'
'enable print e:\tu\20170620-26_700rf.gmf '
i=1
while(i<=28)
'set t 'i
'set grads off'
'set grid off'
'set mpdset hires'
'set mpdset cnworld'
'set lon 10 150'
'set lat 5 85'
'set lev 700'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set gxout shaded'
'set cmin 1'
*'set white -0.2976 0.2976'
'set cint 2'
'set clskip 2'
'set cthick 6'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(TMPprs-273.15))/(TMPprs-29.65)))'
'define qs=(0.62197*es/(prs-0.278*es))'
'define q=qs*rhprs /100'
'define q1=q*ugrdprs/g'
'define q2=q*vgrdprs/g'
'define sqtl=mag(q1,q2)*1000'
'set lon 30 120'
'set lat 10 60'
'd sqtl '
'cbarn 1 1 10.4 4.3'
*'cbarn 1 1'
'set gxout vector'
*'set strmden 2'
'd ugrdprs/g;vgrdprs/g'
'q w2xy 86.06 47.28 '
x1=subwrd(result,3)
x2=subwrd(result,6)
'draw mark 3 'x1' 'x2' 0.10'
'run d:\opengrads\contents\resources\scripts\cbarn.gs'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
'draw title 700hPa rf'recs''
'enable print e:\tu\'recs'700hPa-rf.gmf'
'print'
'c'
i=i+1
endwhile
'disable print'
;
|
-
|