移动坐标下的aave

Jepeng

各位大神，最近小弟挣扎于做一个以台风路径为中心的整层水汽通量，用的是aave函数，但是它积分的区域是在变化的，下面是我的gs文件（其中get_vflux.ctl是已经算好的整层水汽通量，我想把算出来的结果存放到output_vflux.dat文件并写上相应的ctl文件用来做时间序列图，track.txt是台风的经纬度信息）
'reinit'
'open D:\get_vflux.ctl'
'set fwrite D:\output_vflux.dat'
'set gxout fwrite'
i=1
while(i<=25)
read_file=read('D:\track.txt')
read_line=sublin(read_file,2)
lon0=subwrd(read_line,1)
lat0=subwrd(read_line,2)
*say lon0 ' ' lat0
minlon=lon0-8
maxlon=lon0+8
minlat=lat0-8
maxlat=lat0+8
say minlon ' ' maxlon ' ' minlat ' ' maxlat
'set t 'i
'set z 1'
'set x 1'
'set y 1'
'define invap=aave(vflux,lon='minlon',lon='maxlon',lat='minlat',lat='maxlat')'
'd invap'
i=i+1
endwhile
'disable fwrite'
'reinit'

下面是自己写的output_vflux.ctl
dset D:\output_vflux.dat
options byteswapped
undef 1.e30
title INTEG WATER VAPOR INFLOW
xdef  1 linear 1 1
ydef  1 linear 1 1
zdef  1 linear 1 1
tdef 25 linear 12Z24AUG2012 180MN      
VARS 1
invap 1 0 NET WATER VAPOR INFLOW
ENDVARS

最后想出的图类似于这篇文献

但是我做出来的图是这样

这个量级明显有问题，但是我直接生成的整层水汽通量图是正常的，如下

希望有大神可以给予我帮助，谢谢~