我用fnl1*1资料做1000——300整层水汽通量积分，哪位大神们帮我改改，运行不出来啊？

白考儿54

我用fnl1*1资料做1000——300整层水汽通量积分，哪位大神们帮我改改，运行不出来啊？或者推荐个更好用的整层水汽通量积分gs？谢谢啦！
'reinit'
'open fnl_20151211_12_00.ctl'
'set mpdset cnworld'
'set map scaled'
'draw map'
'set lon 70 150'
'set lat 15 70'
'set t 1'
s=24
while(s<=24)
'set z 's''
'set grads off'
'set grid off'
'define p=lev'
'define e=exp(17.67*(tmpprs-273.15)/(tmpprs-29.65))*(RHprs/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=ugrdprs*q/9.8'
'define q2=vgrdprs*q/9.8'
'set gxout barb'
'set cint 5'
'd q1;q2'
'run cbarn.gs'
s=s+1
endwhile
'gxprint guv_all_1112.png white x1000 y800'

D.B.L

zhangxg06 发表于 2016-1-24 20:42
xeixiefenxiang

他这是错的，整层水汽都没积分

liutong

把ctl贴出来看看呗

feixue_ren

还是不会用，在哪可以看到教这些个代码？

白考儿54

feixue_ren 发表于 2016-4-26 09:24
还是不会用，在哪可以看到教这些个代码？

有书啊  买本书或者网上有电子版书  可以看