- 积分
- 141
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
'reinit'
*==原始数据文件==*
'open i:/GFS/week/20120219/gfs_18z.ctl'
lon1='117 113 110 105.5 105.5 110 110 113 120 113 110 116'
lon2='122 117 113 110 110 113 113 120 125 120 116 120'
lat1='22 20 20 19 15 18 15 19 19 15 8 8'
lat2='26 22.3 22 21.7 19 20 18 22 22 19 15 15'
s=1
while(s<=12)
*==各海区起止经纬度==*
lon01=subwrd(lon1,s)
lon02=subwrd(lon2,s)
lat01=subwrd(lat1,s)
lat02=subwrd(lat2,s)
*f2='i:/GFS/Sea_ground.txt'
*bb=read(f2)
*bb1=sublin(bb,2)
*lon1=subwrd(bb1,1)
*lon2=subwrd(bb1,2)
*lat1=subwrd(bb1,3)
*lat2=subwrd(bb1,4)
'set mpdset hires'
'set lon 'lon01' 'lon02''
'set lat 'lat01' 'lat02''
*'set map 4 1 15'
t0=1
while(t0<=65)
f1='i:/gfs/week/20120219/datetime.txt'
aa=read(f1)
aa1=sublin(aa,2)
hour=subwrd(aa1,2)
day=subwrd(aa1,1)
'set t 't0
'set grads off'
*'enable print i:/GFS/week/picture/new/Windspeed.'hour'Z'day'Feb2012.gmf'
'set csmooth on'
******设置风向标大小set digsize******
'set digsize 0.07'
'set ccolor 7'
'set gxout barb'
'set cmark 0'
'set cthick 6'
'set rbcols auto'
*'d skip(UGRD10m,2);skip(VGRD10m,2)' 处理矢量的间隔
*'d UGRD10m;VGRD10m' 默认风矢量化成黑色
**********将风矢量画成彩色**************
'd UGRD10m*2.5;VGRD10m*2.5;mag(UGRD10m,VGRD10m)'
*'d mag(UGRD10m,VGRD10m)'
'draw title Windspeed at 10m 'hour'Z'day'Feb2012 (m/s)'
*'draw string 9.58 6.47 (`3.`1C)'
*'draw string 9.60 7.31 (m/s)'
'set line 2 1 8'
'run cbar.gs'
*'run d:/opengrads/contents/resources/scripts/cbar.gs'
'printim i:/gfs/result/'t0'.Windspeed_Sea's'.'hour'Z'day'Feb2012.gif white'
*'print'
'c'
*'disable print'
t0=t0+1
endwhile
s=s+1
endwhile
'reinit'
以上程序作了两重循环后,在printim的时候'hour'和'day'就循环不了了,
grads提示如下:
如果去掉最外面那一重循环,就米有问题,但这样只能画一个区域的,偶想实现画不同区域的,该怎么办呢?
|
|