关于可降水量单位问题

hustlh

最近在做可降水量的计算，用的是ncep/ncar再分析数据的shum.mon.mean,pres.mon.mean，用的公式是：

                                                                       
其中，w为整层大气可降水量（也称整层水汽含量），q为比湿，ps为地表气压，p取300hpa，g为重力加速度，ρ为液态水密度。

计算GS如下：

'reinit'
'set grid off'
'set grads off'
'sdfopen g:\sqtl\shum.mon.mean.nc'
'sdfopen g:\sqtl\pres.mon.mean.nc'
'set lon 70 140'
'set lat 15 55'
'set z 1 8'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld'
'set map 1 1 1'
'draw map'
i=625
while(i<=636)
'set t 'i''
'set z 1'
'define w=vint(pres.2(z=1),shum.1,300)/1000'
'set gxout shaded'
'd w'
'cbarn'
'set gxout contour'
'd w'
'printim g:\kjs\'i'.tif x1000 y800 white'
'print'
'c'
i=i+1
endwhile
'disable print'
;


不知道怎么搞的，最后求得2000年1-12月各月可降水量（单位是mm）怎么比我已有的降水量数据小，也不知道是不是在求w的时候单位数量级搞错了。

求得的可降水量w大约是在5-55mm范围内，但我已有的降水量数据在6,7,8,9,10月的值都在100mm以上，7月的值超过了200mm。

不知道怎么搞的，降雨数据应该没有问题，按理说可降水量应该是大于降水量的

望大家帮我看看！！！


file:///C:\Documents and Settings\Administrator\Application Data\Tencent\Users\648465979\QQ\WinTemp\RichOle\HZUV%2PADJ4WZ]{B69~VWGB.jpg
file:///C:\DOCUME~1\ADMINI~1\LOCALS~1\Temp\msohtml1\01\clip_image002.giffile:///C:\DOCUME~1\ADMINI~1\LOCALS~1\Temp\msohtml1\01\clip_image002.gif

木林森(lynn)

我觉得是不是资料的不一致性导致出现你这样的问题啊，你算可降水量用的是再分析资料，而降水量用的是观测资料，所处才出现这样的问题。具体我也不知道 ，瞎猜的

暮雪千山

大气可降水量数值范围是没有问题的，哪怕是热带地区，大气可降水量可能也就是60-70mm，但是降雨量大多了。7楼说的比较对，降水的水气不一定是本地的大气可降水量，有可能是其他地方输送过来的。。所以，降雨量大于大气可降水量是没问题的啊

hustlh

大家快来啊！

hustlh

，大家都不来啊！

liuzf406

也遇到了这样的问题， 咋没人来指点一下呢？

不相信

计算整层可降水量的这个问题已经解决好了吗？我也需要求助了

sunlight

是不是 这个仅能算出来的是空气中水汽凝结成的水
对于已经凝结的液态固态水 没有纳入计算？

yuejun0026

可降水量只是这个区域上方空气块中的水汽含量，降水量不仅仅只是可降水量，在降水过程中还有不断的水汽输送过来降水的，所以可降水量小于降水量是很正常的

无敌泡

楼主，我按照你的算法，出来的才是零点几，怎么办啊

datan

降水量就应该大于可降水量，这正常啊