加入地形的剖面图循环出图只有第一张能加入地形

Crystal_qiong

加入地形的剖面图想循环出图，但只能第一张图上能加入地形，其他的都没有，哪位大神能帮下我，万分感谢
gs文件如下：'reinit'
'open E:\lt\meri\orog.ctl'
'sdfopen E:\lt\data\uwnd.1993.nc'
'sdfopen E:\lt\data\omega.1993.nc'
'set strsiz 0.2'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set grid off'
'set xlint 5'
'set ylint 100'
'set ylopts 1 6 0.18'
'set zlog on'
'set lon 75 105'
'set lat 38'
'set t 1'
'set lev 1000 100'
'set gxout shaded'
'set rbcols 1'
'set cmin 0'
'd smth9(orog.1)'

'set dfile 2'
'set lon 75 105'
'set lat 38'
i=1169
while (i<=1179)
'set t 'i''
'set lev 1000 100'
'define u=uwnd.2'
'define w=(omega.3)*(-100)'
'set gxout stream'
'set arrowhead 0.4'
'set strmden 10'
'd u;w'
'draw title Zonal circulation along longitude 38N'
'draw xlab Longitude';'draw ylab Pressure (hpa)'
'printim E:\lt\zonal\'i'(38).png white '
'c'
i=i+2
endwhile
'reinit'
;

四叶草

脚本显示'd smth9(orog.1)'命令只运行过一次，肯定只有第一张图会有地形信息。

Crystal_qiong

四叶草 发表于 2018-11-16 09:43
脚本显示'd smth9(orog.1)'命令只运行过一次，肯定只有第一张图会有地形信息。

但是地形文件只有一个时次，把它放在循环里面也会提示错误，应该怎样解决这个问题呢？

四叶草

Crystal_qiong 发表于 2018-11-16 09:47
但是地形文件只有一个时次，把它放在循环里面也会提示错误，应该怎样解决这个问题呢？

变量带着时间参数，时间参数永远显示t=1'd smth9(orog.1,(t=1))'

Crystal_qiong

四叶草 发表于 2018-11-16 09:58
变量带着时间参数，时间参数永远显示t=1'd smth9(orog.1,(t=1))'

还是不行啊

四叶草

Crystal_qiong 发表于 2018-11-16 10:11
还是不行啊

应该是变量名有问题，你把'set dfile 2'这个命令去掉看看。理论上是行得通的，我自己也没实践过。

Crystal_qiong

四叶草 发表于 2018-11-16 10:16
应该是变量名有问题，你把'set dfile 2'这个命令去掉看看。理论上是行得通的，我自己也没实践过。

行不通啊，错误更多了

四叶草

Crystal_qiong 发表于 2018-11-16 10:23
行不通啊，错误更多了

那就单时次绘制得了

付亚男

'reinit'
'sdfopen E:\DATA\NC_2P5\uwnd.2017.nc'
'sdfopen E:\DATA\NC_2P5\omega.2017.nc'
'open E:\DATA\OROG\orog.ctl'


'set lon 70 110'
'set lat 35 40'
'set lev 1000 100'
'set t 1 20'

'define u=uwnd'
'define w=(omega.2)*(-100)'



i=2
while (i<=9)


'set lon 75 105'
'set lat 38'
'set lev 1000 100'
'set t 'i

'set gxout stream'
'd u;w'


'set gxout shaded'
'set rbcols 1'
'set cmin 0'
'set csmooth on'
'd orog.3(t=1)'


'printim E:\'i'.png white '

'c'
i=i+1
endwhile

;

付亚男

试试像我在楼上这么写，应该是可以的，我试过了但是你的omega应该是不可以直接和uwnd画流线的，因为前者的单位是Pa/s，最好换算成m/s才有物理意义