各种湿度（水汽压，比湿，露点，相对湿度）的换算及位温，假相当位温Fortran源代码

wheelar

一般只要知道一个湿度参量，就可换算成其他的湿度，绝对量有水汽压，比湿，露点，相对量有相对湿度和温度露点差，另外程序里还提供了位温和假相当位温，其中，可以根据比湿或相对湿度或露点计算假相当位温
程序命名规则：FCN_SpecificHumidity_Rh(Tem,Lev,Rh)，FCN代表函数，根据相对湿度（rh）计算比湿SpecificHumidity
以上源代码均经过测试通过检验
例子：根据比湿计算假相当位温：
function FCN_ThetaSe_SpecificHumidity(Tem,Lev,q)
    implicit none
    real(kind=4),intent(in) :: Tem,Lev,q
    real(kind=4) :: FCN_ThetaSe_SpecificHumidity
    real(kind=4) :: Theta,vaporPressure
    real(kind=4) :: sglPow,sglDen,sglLog,sglCoff
    vaporPressure=FCN_VaporPressure(Lev,q)
    Theta=FCN_Theta_SpecificHumidity(Tem,Lev,q)
    sglLog=3.5*log(Tem)-log(vaporPressure)-4.805
    sglDen=55.0+2840/sglLog
    sglCoff=3376.0/sglDen-2.54
    sglPow=sglCoff*q*(1.0-0.81*q)
    FCN_ThetaSe_SpecificHumidity=Theta*exp(sglPow)
    return
end function FCN_ThetaSe_SpecificHumidity

libMoisture.f90 (8.77 KB, 下载次数: 464)

2015-5-23 17:46 上传

点击文件名下载附件

秋之期

wheelar 发表于 2015-5-25 20:16
所有的参数，vaporPressure为是水汽压，Es为饱和水气压，Tem：温度(K),dPoint：露点，Rh相对湿度，q比湿
...

那不知道各个参数的具体单位呢？比湿用的是g/kg，还是g/g？

xinyuan7068

谢谢楼主分享

daboy

能否给出程序中的参数说明？

wangyipu

下来先学习学习

wheelar

daboy 发表于 2015-5-25 09:44
能否给出程序中的参数说明？

所有的参数，vaporPressure为是水汽压，Es为饱和水气压，Tem：温度(K),dPoint：露点，Rh相对湿度，q比湿
Qs，饱和比湿，theta：位温，thetase：假相当位温

heqiurui

很有帮助！谢谢分享

heqiurui

第一个函数用温度代入求饱和水汽压的公式中，if else 中是否缺少一个判断条件，因为这个语句对abc的赋值是错误的

wheelar

heqiurui 发表于 2015-6-3 11:04
第一个函数用温度代入求饱和水汽压的公式中，if else 中是否缺少一个判断条件，因为这个语句对abc的赋值是 ...

谢谢提醒，只是这里关于饱和水汽压的计算，不同教材好像对系数确实不一样，其中有一种是根据以温度-12℃为界分为两种，还有直接就一种系数，我是两种都试过，好像差别并不大。可能是程序里面ifelse写的不明白，那块基本就没有用。不过即使就这么运行，虽然饱和水汽压无法验证，但通过饱和水汽压计算的比湿和原来模式输出的比湿比较，几乎是一样的，所以饱和水汽压应该也不会有太大问题

秋之期

挺好的代码