整层水汽输送通量单位换算

自由行

    利用ncep资料2.5*2.5计算整层水汽输送后，计算最终换算为kg/m/s.所以由本来计算出的g/S/CM换算的话，只用除以10，然后由于积分时vint函数外还需除以9.8，所以最后结果是除以10*9.8（10*10），但是计算结果不对，反而应该总共除以1000.不知道为什么真是奇怪，急，请大神和版主帮忙。
这是最相关gs
reinit'
'open E:\jiduan\rain_reason\shuiqitongliang\strongyear\2.2.2.2\1quv.ctl'
'open E:\jiduan\rain_reason\shuiqitongliang\strongyear\2.2.2.2\1p.ctl'
'set gxout fwrite'
'set fwrite E:\jiduan\rain_reason\shuiqitongliang\strongyear\2.2.2.2\zheng-quv.grd'
'set lon 60 160'
'set lat 10 50'

it=1
while(it<=481)
  'set t 'it''
  'd vint(ps.2,fqu,300)/(9.8*10)'
  'd vint(ps.2,fqv,300)/(9.8*10)'
*vint表示积分，ps→300表示整层大气的积分
*(整层水汽通量的单位：kg/(m*s)--------*
   it=it+1
endwhile

'disable fwrite'
'reinit'

kongfeng0824

是2.5度的数据吗？那套数据呢？有网址吗？

羊羊羊

你一开始计算出来的单位应该是g*s/kg

15951676682

降水资料：https://www.esrl.noaa.gov/psd/data/gridded/data.cpc.globalprecip.html

其他资料：
https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html