求整层水汽和逐层水汽

小甜甜

    我要求整层水汽，850-500,和500-200的水汽量
    这是求整层水汽的gs,
'reinit'
'sdfopen E:\data\specific-humidity\shum.day.1981-2010.ltm.nc'
'sdfopen E:\data\pres.sfc\pres.day.1981-2010.ltm.nc'
'set fwrite E:\bk2013y\w.1981-2010.ltm.nc'
'set gxout fwrite'
'set x 0 360'
'set y -90 90'
'set z 1 '
i=1
while(i<=365)
'set t 'i
'define wq=vint(pres.2(z=1),shum.1,300)'
*'define w=wq*1000'
'd wq'
i=i+1
endwhile
'disable fwrite'
;
‘shum.day.1981-2010.ltm.nc'：

pres.day.1981-2010.ltm.nc'

求高人指点，怎么做才对

曾记否

求850-500的，这个地方改了就好啊'define wq=vint(pres.2(z=1),shum.1,300)'
那句话是指从底层积分到300hpa

小甜甜

曾记否 发表于 2013-4-10 21:44

                               
登录/注册后可看大图

求850-500的，这个地方改了就好啊'define wq=vint(pres.2(z=1),shum.1,300)'
那句话是指从底层积分到300hp ...

我贴的程序画出来图本身就不对，求指导。而且你说求850-500的，改'define wq=vint(pres.2(z=1),shum.1,300)'
请问该怎么改，pres.2是地面气压，是变量，不能改成850常量。

笨北极熊笨

地面气压的单位是hPa还是Pa？若是Pa，这样写就不对。

小甜甜

笨北极熊笨 发表于 2013-4-12 16:55

                               
登录/注册后可看大图

地面气压的单位是hPa还是Pa？若是Pa，这样写就不对。

单位是pa,但我就是除以100也不对

笨北极熊笨

小甜甜 发表于 2013-4-13 10:26

                               
登录/注册后可看大图

单位是pa,但我就是除以100也不对

set z 1要加吗？下面是我写的积分gs，是对的
********从地面到100百帕对动能积分********
'reinit'
'open h:\data\udat\1979u.ctl'
'open h:\data\vdat\1979v.ctl'
'open h:\data\pdat\1979p.ctl'
'set gxout fwrite'
'set fwrite h:\data\vint\1979.dat'
'set lon 40 160'
'set lat 10 60'
nt=1;while(nt<=365)
'set t 'nt
'define k=vint(p.3/100,(u*u+v.2*v.2)/2,100)'
'd k'
nt=nt+1;endwhile
'disable fwrite'

7777

我也有相关的问题 不知道怎么处理

daodang2008

都在学习啊，这个好难搞的样子，能不能直接输出txt或者excel文件啊

15951676682

vint*9.8

quezhiping

楼主，问题解决了吗？好像vint函数只能计算从地面开始的，不能从中间层次开始积分，请问你850-500和500-200是怎么计算的呢