- 积分
- 2276
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 278803532 于 2014-5-19 10:22 编辑
前几个星期在做水汽通量,看了很多论坛上贴着,也下载了很多GS,然后被清风版主封号了 ,后来根据别人的提示,修改了了程序,也看了http://www.iges.org/grads/gadoc/gadocindex.html这个网站一些函数的问题,搞定之后,发现以前做水汽通量散度的时候很多人做出来不是10^-7量级的,然后看了一些他们的程序,发现一点小问题,现在福利回送了。其实也就是一个单位的问题,具体如下,大家都忽略掉的一个就是算散度的时候,GRADS里面的函数用单位是米,但是水汽通量里面的用的是CM。
*---------------------------------------------------------------------------
*水汽压单位hPa,hPa=kg*cm^-1*s^-2,气温单位℃,比湿单位g/kg,由原来kg/kg乘以1000得到
*水汽通量单位g/(cm*hPa*s),uu,vv算整层水汽通量积分用到,单位(g*m)/(kg*s)
*水汽通量散度单位:g/(cm^2*hPa*s),除以hdivg(uq,vq)/100,除以100,考虑计算散度时,
*地球半径用的单位是m,变为cm,所以要除以100,
*----------------------------------------------------------------------------
程序大家都懂,就不多说了,描述文件我只了一个变量的,其他变量描述一样的,我就不重复了,整层积分那个函数用积分 'd vint(pres.2, qvv,500)',这里面变量QVV,对应水汽通量程序里面的'define vv=ugrdprs(t='it')*q',具体看上面网页那个描述。如果大家还有什么单位运算不明白的话,可以看《天气学原理》水汽那章,有单位的运算,或者《天气分析预报物理量计算基础》
|
评分
-
查看全部评分
|