怎样画任意两点的时间剖面

冰枫

路边摊1990 发表于 2012-5-28 13:03

                               
登录/注册后可看大图

我已经画出来了，现在我的问题是想和我发的图一样在x轴上即写经度又写纬度

这个你就可以查看论坛里的剖面图的GS了，它上面有写坐标的语句

fish

路边摊1990 发表于 2012-5-28 13:03

                               
登录/注册后可看大图

我已经画出来了，现在我的问题是想和我发的图一样在x轴上即写经度又写纬度

那种坐标貌似grads画不出来，我是用一种比较笨的方法给写上去的~

流浪的心

利用这脚本函数anghov.gs

anyback

需要需要~

river

楼主共享一下方法吧，大家都想知道呢

river

楼主分享下方法吧，让大家学习一下

路边摊1990

river 发表于 2012-11-25 22:05

                               
登录/注册后可看大图

楼主分享下方法吧，让大家学习一下

'open f:\qcat\spd2009.ctl'
'set lon 120 140'
'set lat 20 40'
'set t 60'
*'define windspd=sqrt(u*u+v*v)'
'd windspd'
say 'click a'
'q bpos'
x0=subwrd(result,3)
y0=subwrd(result,4)
say 'click b'
'q bpos'
x1=subwrd(result,3)
y1=subwrd(result,4)
'draw line 'x0' 'y0' 'x1' 'y1''
'q xy2w 'x0' 'y0''
lon0=subwrd(result,3)
lat0=subwrd(result,6)
'q xy2w 'x1' 'y1''
lon1=subwrd(result,3)
lat1=subwrd(result,6)
lx=(lon1-lon0)*112
ly=(lat1-lat0)*223
'd sqrt('lx'*'lx'+'ly'*'ly')'
l=subwrd(result,4)
*decimal=l/55.8
decimal=l/50
decimal_part=substr(decimal,3,20)
n=decimal-decimal_part
inclat=(lat1-lat0)/n
inclon=(lon1-lon0)/n
i=1
while(i<=n)
lon.i=lon0+(i-1)*inclon
lat.i=lat0+(i-1)*inclat
i=i+1
endwhile
'set fwrite f:\line\hov.dat'
nt=1
while(nt<=91)
i=1
'set t 'nt''
while(i<=n)
'set lon 'lon.i''
'set lat 'lat.i''
'set gxout fwrite'
'd windspd'
*dummy=subwrd(result,4)
*data=subwrd(dummy,4)
*rc=write('d:\gad.dat',dummy,append)
i=i+1
endwhile
nt=nt+1
endwhile
'disable fwrite'
fname = 'f:\line\hovtmp.ctl'
rc = write(fname,'dset d:\line\hov.dat')
rc = write(fname,'undef -999.000')
rc = write(fname,'title Time Mask')
rc = write(fname,'xdef 'n' linear 'lon1' 'lon.1' 'lon.2' 0.25')
rc = write(fname,'ydef 1 linear 'lat0' 'lat.1' 'lat.2' 0.25')
rc = write(fname,'zdef 1 linear 0 1')
rc = write(fname,'tdef 18  linear 00Z003NOV2007 7dy')
rc = write(fname,'vars 1')
rc = write(fname,'tmask 0 99 Time Mask')
rc = write(fname,'endvars')
rc = close(fname)
*close('d:\gad.dat')
*return
;

river

路边摊1990 发表于 2012-11-27 17:39

                               
登录/注册后可看大图

'open f:\qcat\spd2009.ctl'
'set lon 120 140'
'set lat 20 40'

很强大，谢谢分享

river

路边摊1990 发表于 2012-5-28 13:03

                               
登录/注册后可看大图

我已经画出来了，现在我的问题是想和我发的图一样在x轴上即写经度又写纬度

使用如下语句
'set xlabs 35N,70E|37N,82E|39N,94E|41N,106E|43N,118E|45N,130E'

007wojiushi

做个霸气的木头 发表于 2012-5-26 19:29

                               
登录/注册后可看大图

Grads里面设定经纬度没有间隔设定的功能，有一种笨的方法可以把这图画出来：
你把需要的直线上的每个点每个 ...

嗯嗯 这个方法理论上是可以的！