- 积分
- 526
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-6
- 最后登录
- 1970-1-1
|
发表于 2020-10-7 00:12:23
|
显示全部楼层
老师,为什么我调用cbar_matlab.gs,总是显示不能画出color bar
'reinit'
'open L:\DJBY\ca_saf\201801\20180101.ctl'
i=1
while(i<=21)
'set t 'i''
'set gxout shaded'
'set mpdset cnworld'
'set map 9 1 1'
'set csmooth on'
'set grads off'
'define RH=rhprs'
'set lat 20 35'
'set lon 95 115'
'set lev 850'
'define g=9.8'
'define es=6.112*exp(17.67*(TMPprs-273.16)/(TMPprs-29.65))'
'define qs=0.62197*es/(lev-0.378*es)'
'define q=qs*RHprs/100'
'define qhdivg=hdivg(q*ugrdprs,q*vgrdprs)*1e8'
'define q1=UGRDprs*qs'
'define q2=VGRDprs*qs'
'set gxout shaded'
'set cint 5'
'set lat 21 30'
'set lon 96 107'
'L:\DJBY\ca_saf\201801\cbar_matlab.gs'
'd qhdivg'
'set gxout barb'
'set ccolor 2'
'd UGRDprs*2.5;VGRDprs*2.5'
'q dim'
it=sublin(result,5)
it=subwrd(it,6)
'draw title 'it'-850hPa-waterflux'
'printim L:\DJBY\ca_saf\201801\waterflux_850hPa1\'it'-850hPa-waterflux1.png white'
c
i=i+1
endwhile
'disable gxprint'
;
|
|