- 积分
- 900
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-7-18
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
以下是我批量绘制500hpa高度场的gs脚本,一共14时次,
画出来是14张图,但是从第九张开始,高度场没有数字,都是0
我查过我的数据是没有问题的,请教各位是不是gs脚本的问题,这个问题出在哪里呢?
'reinit'
'open f:\grads\times\500_hgt_wind\hgt_wind_500_times.ctl'
'enable print f:\grads\times\500_hgt_wind\2004_hgt.gmf'
'set dfile 1'
ti=1
while(ti<=14)
'set lat 10 60'
'set lon 80 140'
'set lev 500'
'set t 'ti
'set gxout contour'
'set cint 60'
'set cthick 6'
' d HGTprs'
'set mpdset mres cnriver'
*'set mpdset cnriver'
'set map 4 1 6'
'draw map'
*string1=time1
* string2='500 hPa'
* rc=DrawStr(string1,string2)
ti=ti+1
'print'
'c'
endwhile
*'printim f:\grads\test\20070707_500dayhgtave.png'
'disable print'
*FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF
function DrawStr(string1,string2)
'q gxinfo'
xarea=sublin(result,3); yarea=sublin(result,4)
xstr =subwrd(xarea,4) ; ystr =subwrd(yarea,6)
'set strsiz 0.12'; 'set string 1 l 6 0'
'draw string 'xstr+0.1' 'ystr+0.30' 'string1
'draw string 'xstr+2.5' 'ystr+0.30' 'string2
* 'draw string 'xstr+1.2' 'ystr+0.25' 'string3
return
|
|