vint函数的用法

凉茶

本人编了一个求整层水汽积分的gs，用到vint（slp，shum，top）函数，经过调试，也出了数据，从原理上讲所出图应该和对流层下层的水汽含量差不多吧 可是我出的数据却大的离奇，检查了好多遍都没看出毛病来，我在想是不是vint函数的用法我用的不对，求大师指导下！以下是我的gs文件！
'reinit'
'open D:\lunwen\zcjf\shum.ctl'
'open D:\lunwen\zcjf\slp.ctl'
'set gxout fwrite'
'set grads off'
'set fwrite d:\lunwen\zcjf\vint.grd'
'set lon 70 140'
'set lat 0 60'
i=1
while(i<=54)
'set t 'i''
'define zcjf=vint(slp.2,shum.1,300)'
'd zcjf'
i=i+1
endwhile
'disable fwrite'
;

做个霸气的木头

你slp.2的气压单位是hpa吗？
话说top和slp的单位必须为hpa的

易小凯

你的量级大概是多少？你算的是整个气柱中的水汽量，量级大一些应该也是正常的

凉茶

易小凯 发表于 2012-5-11 13:15

                               
登录/注册后可看大图

你的量级大概是多少？你算的是整个气柱中的水汽量，量级大一些应该也是正常的

大一些正常没错 但能达到3个数量级么 大气中水汽不是大部分集中在对流层下层么 按道理来说应该和850 的水汽含量差不多啊  

凉茶

凉茶 发表于 2012-5-11 13:49

                               
登录/注册后可看大图

大一些正常没错 但能达到3个数量级么 大气中水汽不是大部分集中在对流层下层么 按道理来说应该和850 的水 ...

为了验证所作的正确性 我把vint（slp，shum，top）中top的值设置为900 得出的结果也是原来湿度大小的1000倍左右 是不是我单位弄错了

凉茶

申明一下 我所用的资料是necp的海平面气压月平均资料 单位应该就是hap吧

凉茶

做个霸气的木头 发表于 2012-5-11 12:33

                               
登录/注册后可看大图

你slp.2的气压单位是hpa吗？
话说top和slp的单位必须为hpa的

我所用的资料是necp的海平面气压月平均资料 单位应该就是hap吧 而我最高层 即top也是直接写的常量300 也是指300hap高度场啊

做个霸气的木头

凉茶 发表于 2012-5-11 13:55

                               
登录/注册后可看大图

我所用的资料是necp的海平面气压月平均资料 单位应该就是hap吧 而我最高层 即top也是直接写的常量300 也是 ...

为保险起见，你还是输出slip，看看数值是不是hpa单位的，如果输出等值线显示为5980等，明显不是hpa，需要要除以十。

陌小默的

NCEP的资料气压单位不是pha

凉茶

做个霸气的木头 发表于 2012-5-11 18:24

                               
登录/注册后可看大图

为保险起见，你还是输出slip，看看数值是不是hpa单位的，如果输出等值线显示为5980等，明显不是hpa，需要 ...

我把海平面气压文件转成txt 数据很正常啊 但程序运行结束后 出来的数据很大