加载中…grads计算水汽通量及单位
(2018-08-27 16:01:41)
标签：
水气通量、整层水汽通 |
分类： grads |
一、定义。水汽通量即表示水汽输送强度的物理量，定义为：单位时间内流经某一单位面积的水汽质量（用g/s表示，若单位面积取1cm^2,单位高度取hPa，一般单位为g/(s*hPa*cm)）。
二、单位。这里，u单位是m/s，q单位是g/kg，那么计算出来的单位就和天原书上g/(s*hPa*cm)相同，具体推到见天原书P636；若q单位为kg/kg,则水汽通量单位为kg/(s*hPa*cm)。经计算，夏季青藏高原地区500hPa平均水汽通量量级为个位数，量级为10g/(s*hPa*cm)。（Pa=N/m^2=kg*m/m^2/s^2=g/cm/s^2）
三、整层积分。
grads垂直积分时，默认已经*100（从hPa换成Pa，所以自己的地面气压数据代入公式前单位必须换成是hPa）并且默认除以9.8了。根据国家气候中心1000-300的整层积分的量级为800kg/m/s或1.0e5g/m/s。
（见气象家园帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=6383&extra=&page=3）
grads代码：
- 'sdfopen
E:\TDDOWNLOAD\nc_data\hgt.1998.nc'
- 'sdfopen
E:\TDDOWNLOAD\nc_data\shum.1998.nc'
- 'sdfopen
E:\TDDOWNLOAD\nc_data\uwnd.1998.nc'
-
'sdfopen E:\TDDOWNLOAD\nc_data\vwnd.1998.nc'
- 'sdfopen
E:\TDDOWNLOAD\nc_data\pres.sfc.1998.nc'
- 'set lat -10 80'
- 'set lon 60 160'
- 'set z 1 8'
- 'set time 12z01jun1998
12z31aug1998'
- ************vapar
flux***************
- 'define
uq=shum.2*uwnd.3'
- 'define
vq=shum.2*vwnd.4'
- 'define
uqs=vint(pres.5(lev=0),uq,300)'
- 'define vqs=vint(pres.5(lev=0),vq,300)‘
NCL积分时vint_qu=vibeta(lev,qu(latitude|:,longitude|:,level|:),linlog,spmon5,pbot,ptop)/g*100中除以g乘以100后才与上面grads中vint(pres.5(lev=0),vq,300)具有同样的量级。
四、水汽通量散度
整层水汽通量散度，单位是g/(m**2*s)，量级是0.1,若单位是kg/(m**2*s)，量级是0.0001.
五、边界积分
整层水汽通量对某一边界积分（求和）后单位为1.0e6kg/s
喜欢
0
赠金笔
前一篇：地形资料
后一篇：NCL特殊字符（上、下标等）