关于整层水汽通量积分的错误

meiyao898

       各位大神，小弟今天遇到了一个问题，自己弄了一下午了还没解决，在此特向各位求助。

       我求整层积分不是向其他人一样，在NC资料中提取，而是自己算出每层水汽通量，而后把每层水汽通量一起写到一个GRD中，然后再对此grd中水汽通量进行积分，具体如下：

       第一步，写每层的GRD，GS如下


如上，计算过后写到每层的grd中，一共9层，分别为1000,925,850,775,700,600,500,400,300 共9个GRD文件

     第二步，把9个每层的grd写到一个GRD钟，GS如下
'reinit'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t1.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t2.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t3.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t4.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t5.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t6.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t7.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t8.ctl'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\t9.ctl'
'set fwrite C:\Users\a\Desktop\123\grads\shuiqi\3\shuiqi.grd'
'set gxout fwrite'
'set lon -19 19'
'set lat -9 14'
'set lev 1000'
'd uq.1'
'd vq.1'
'set lev 925'
'd uq.2'
'd vq.2'
'set lev 850'
'd uq.3'
'd vq.3'
'set lev 775'
'd uq.4'
'd vq.4'
'set lev 700'
'd uq.5'
'd vq.5'
'set lev 600'
'd uq.6'
'd vq.6'
'set lev 500'
'd uq.7'
'd vq.7'
'set lev 400'
'd uq.8'
'd vq.8'
'set lev 300'
'd uq.9'
'd vq.9'

'disable fwrite'

;


描述文件ctl，取其中一个，如下：
dset C:\Users\a\Desktop\123\grads\shuiqi\3\shuiqi1.grd
undef -999.0
title ****
xdef 39 linear -19 1
ydef 24 linear -9 1
zdef 1 linear 1000 1
tdef 1 linear jul1954 1mo
vars 2
uq 0 0 heights
vq 0 0 heights
endvars


第三步，对GRD文件整层积分，GS如下：
'reinit'
'open C:\Users\a\Desktop\123\grads\shuiqi\3\all1.ctl'
'enable print C:\Users\a\Desktop\123\grads\shuiqi\3\sq1.gmf'
'set lat -19 19'
'set lon -19 19'
'set z 1'

'define uqs=vint(1000,uq,300)'
'define vqs=vint(1000,vq,300)'
'define abq=mag(uqs,vqs)'
'd abq'
'print'
'disable print'
;


描述文件如下：
dset C:\Users\a\Desktop\123\grads\shuiqi\3\shuiqi.grd
undef -999.0
title ****
xdef 39 linear -19 1
ydef 24 linear -9 1
zdef 9 levels 1000 925 850 775 700 600 500 400 300
tdef 1 linear jul1954 1mo
vars 2
uq 9 0 heights
vq 9 0 heights
endvars



但是运行后，grads有如下错误





不知为何会这样，请各位高手帮帮忙啊，显示的好像是格点错误，我现在有两个可能感觉有问题的地方，第一，是我把每层的写到一个GRD中除了问题？但是感觉应该没有问题的。  第二，积分  vint(1000,vq,300)这样是不是有问题，难道只能从地面积分才可以么？

超台十八

set z 1 有需要吗？

做这个还是写个 fortran 好。

yeolsuu

我自己写的脚本，运行之后也是显示这个错误， 求解啊

taoqt

为何要每层分一个文件

congxiaochong

解决了么 我也遇到这个问题了

virl

同问{:eb302:}{:eb302:}{:eb302:}

virl

同求同求{:eb302:}

乖乖孙迪

同求，也是遇到这个问题

默戏

同问啊这个问题

DQ153LX

学习了，先看看