【已解决】求助：算垂直水汽积分用vint报错

海soso

我要算的是21年的冬季平均的水汽垂直输送 贴上我用到的Ctl 和gs 求各位大神江湖救急 ！！！gs:
'reinit'
'open d:\wu\waterv\ndjf9415_pres.ctl'
'open d:\wu\waterv\ndjf9415_uq.ctl'
'open d:\wu\waterv\ndjf9415_vq.ctl'

'set fwrite d:\wu\waterv\vint_ndjf9415_uqs.grd'  
*'set fwrite d:\wu\waterv\vint_ndjf9415_vqs.grd'  
'set gxout fwrite'

i=1
while(i<=21)
'set lon 70 140'
'set lat 15 55'
'set t 'i''
'set z 1 8'
'define uqs=vint(pres.1/100,uq.2,300)/9.8'
*'define uqs=vint(uq.2,pres.1/100,300)/9.8'
*'define vqs=vint(pres.1/100,vq.3,300)/9.8'
'd uqs'
*'d vqs'
i=i+1
endwhile
'disable fwrite'
;

然后是ctl 一共三个
第一个：
DSET d:\wu\waterv\ndjf9415_pres.grd
UNDEF -999.0
TITLE NDJF mean
xdef   29 linear   70.00  2.500
ydef   17 linear   15.000  2.500
zdef   1 levels 1000
TDEF 21 linear 00z01jan1994 1yr
VARS 1
pres 0 99 Monthly Mean of pressure
ENDVARS

第二个：
DSET d:\wu\waterv\ndjf9415_uq.grd
UNDEF -999.0
TITLE NDJF mean
xdef   29 linear   70.00  2.500
ydef   17 linear   15.000  2.500
zdef   8 levels 1000 925 850 700 600 500 400 300
TDEF 21 linear 00z01jan1994 1yr
VARS 1
uq 8 99 Monthly Mean of uq
ENDVARS

第三个：
DSET d:\wu\waterv\ndjf9415_vq.grd
UNDEF -999.0
TITLE NDJF mean
xdef   29 linear   70.00  2.500
ydef   17 linear   15.000  2.500
zdef   8 levels 1000 925 850 700 600 500 400 300
TDEF 21 linear 00z01jan1994 1yr
VARS 1
vq 8 99 Monthly Mean of vq
ENDVARS


grads报的错误是：


谢谢各位！！！


【已解决】 其实不能算已解决 是我从别人那里拿了一个程序 是没问题的  贴出来吧：
'reinit'
'open d:\wu\waterv\ndjf9415_pres.ctl'
'open d:\wu\waterv\ndjf9415_u.ctl'
'open d:\wu\waterv\ndjf9415_v.ctl'
'open d:\wu\waterv\ndjf9415_shum.ctl'

'set fwrite d:\wu\waterv\vint_ndjf9415_qs.grd'
'set gxout fwrite'

i=1
while(i<=21)

'set dfile 1'
'set z 1'
'set t 'i''
'set lon 70 140'
'set lat 15 55'
'define p1=pres'
'set dfile 2 3 4'
'set z 1 8'
'set t 'i''
'set lon 70 140'
'set lat 15 55'
'define q1=shum.4';'define u1=uwnd.2';'define v1=vwnd.3'
'define qu1=q1*u1';'define qv1=q1*v1'
'set dfile 2 3 4'
'set z 1'
'set t 'i''
*'d vint(p1,qu1,300)'
*'d vint(p1,qv1,300)'
'd hdivg(vint(p1,qu1,300),vint(p1,qv1,300))'

i=i+1
endwhile

'disable fwrite'
*ret=close('e:\shuju\cyclone\850\2013-8.txt')
;

如果有人发现我原来的程序的问题还是欢迎指出:-D

男紫汗

pres.1/100这个为什么除以100？

石榴之爱

'define uqs=vint(lev(z=1),uq,300)/9.8'

海soso

石榴之爱 发表于 2017-4-1 07:30
'define uqs=vint(lev(z=1),uq,300)/9.8'

'define uqs=vint(pres.1(z=1),uq.2,300)/9.8'
我改成了这样 还是不行 报错是一样的T.T

男紫汗

'set z 1 8'放到循环外面，循环里面set z 1

海soso

男紫汗 发表于 2017-4-1 09:48
'set z 1 8'放到循环外面，循环里面set z 1

'reinit'
'open d:\wu\waterv\ndjf9415_pres.ctl'
'open d:\wu\waterv\ndjf9415_uq.ctl'
'open d:\wu\waterv\ndjf9415_vq.ctl'

'set fwrite d:\wu\waterv\vint_ndjf9415_uqs.grd'  
*'set fwrite d:\wu\waterv\vint_ndjf9415_vqs.grd'  
'set gxout fwrite'

'set z 1 8'
i=1
while(i<=21)
'set lon 70 140'
'set lat 15 55'
'set t 'i''
'set z 1'
'define uqs=vint(pres.1(z=1),uq.2,300)/9.8'
*'define vqs=vint(pres.1/100,vq.3,300)/9.8'
'd uqs'
*'d vqs'
i=i+1
endwhile
'disable fwrite'
;
改成了这样 还是一样的报错

石榴之爱

'set z 1 8'去掉试试

石榴之爱

'set x 1 144'
'set y 1 73'
i=1
while(i<=34)
'set t 'i
'set lev 1000'
'define qq1=vint(lev(z=1),q1,300)/9.8'
'define qq2=vint(lev(z=1),q2,300)/9.8'
'd qq1'
'd qq2'

i=i+1
endwhile
我是这样的

seafairy

我试过了，'set dfile 2 3 4'，这一个是关键