gs脚本编辑

独行客

vel.gs内容
'reinit'
'set gxout fwrite'
'set fwrite d:\FNL\vel.dat'
'open D:\FNL\fnl201405.ctl'
i=28
while(i<=34)
'define vel=aave(vvelprs,lon=112,lon=115,lat=22,lat=23)'
'set t 'i''
'set lev 1000'
'd vel'
'set lev 900'
'd vel'
'set lev 800'
'd vel'
'set lev 700'
'd vel'
'set lev 600'
'd vel'
'set lev 500'
'd vel'
'set lev 400'
'd vel'
'set lev 300'
'd vel'
'set lev 200'
'd vel'
'set lev 100'
'd vel'
i=i+1
endwhile
'disable fwrite'
'reinit'
然后运行文件 产生vel.dat文件
接着编写
vel.ctl内容
DSET d:\FNL\vel.dat
TITLE Upper Data
undef -9.99E33
xdef 4 linear 112 1.0
ydef 2 linear 22 1.0
zdef 10 levels 1000 900 800 700 600 500 400 300 200 100
tdef 7 linear 0z21May2014 6hr
vars 1
vel 10 99 Pressure vertical velocity
Endvars
最后
omega.gs文件
'reinit'
'open d:\FNL\vel.ctl'
'set x 1'
'set y 1'
'set lev 1000 100'
'set t 1 7'
'define_colors'
'set cterp on'
'set csmooth on'
'set grads off'
'set grid off'
'd vel'
'enable print d:\fnl\omega0.gmf'
'print'
'disable print'
'reinit'

独行客

运行第一个文件时，得到vel.dat数据  可以有定义语法错误

独行客

独行客 发表于 2014-9-27 13:50
运行第一个文件时，得到vel.dat数据  可以有定义语法错误

估计是定义变量出错啦 请高手能帮我看看吗

独行客

{:5_275:}{:5_275:}{:5_275:}{:5_275:}

传说中的谁

你做区域平均后不是一个数据么，ctl里怎么x方向有4个点？

传说中的谁

你做区域平均后不是一个数据么，ctl里怎么x方向有4个点？

智文

楼主辛苦了，资料很全面，非常好！

千里草

Grads中好像不可以连续使用d命令的吧，你的结果正确吗？

liu541226

aave的语法错误了吧？怎么选的是区域平均

对雪吟歌

既然已经做了区域平均，怎么经纬度还是不固定呢？你不是区域平均了吗