我使用这个广泛接受方程计算RH (%):
$ RH = e \ * \压裂{100}{es (T)} $
在哪里e美元蒸汽压力和吗es (T)美元在温度的饱和蒸汽压力吗元新台币。
我已经观察到e美元(从0$ - $30 hPa)。我计算es美元在R使用下列方程(取决于湿球温度Tw美元相对于零(占蒸汽压力超过液体或固体水)):
f.es1 < -函数(T) * exp (17.38 * 6.107 T / (239。+ T) # Tw > = 0 f.es2 < -函数(T) 6.107 * exp (22.44 * T / (272.4 + T)) # Tw < 0
我注意到,RH美元(使用上面的第一个方程)非常大在某些情况下(例如,4352.567)。这是发生在e美元是0.6,es美元是0.01378497。
我知道我可以规模RH 1:10 0之间的数据,但我想知道是否有更好的办法处理这个吗?它发生时e美元和es美元是极其微小的。我猜这也许真的是在干旱地区吗?有更好的方法计算吗RH美元这些地方吗?
这是一些关于数据集的信息:
1变量(不包括维度变量):短vap(经度、纬度、时间)(分块:[1]2160年,30日)standard_name: vapor_pressure long_name: vapor_pressure单位:kPa add_offset: 0 scale_factor: 0.01 _FillValue: -32768 missing_value: -32768描述:蒸汽压维度:经度纬度时间coordinate_system: WGS84, EPSG: 4326三维:尺寸:12 * * *是无限的* * * standard_name:时间long_name:时间单位:天1900-01-01就是日历:公历轴:T朗大小:2160 standard_name:经度long_name:经度单位:degrees_east轴:X lat大小:1080 standard_name:纬度long_name:纬度单位:degrees_north轴:Y 5全球属性:CDI:气候数据接口版本1.9.9rc1 (https://mpimet.mpg.de/cdi)约定:cf - 1.6历史:2021年9月15日14:21:54结婚:cdo remapcon, r2160x1080 TerraClimate19812010_vap。数控WCvap_terraclimate_1981_2010。数控方法:这些层从TerraClimate创建使用气候上辅助插值的每月异常CRU Ts4.0和日本55年再分析与WorldClim v2.0气候学(JRA-55)数据集。CDO:气候数据运营商版本1.9.9rc1 (https://mpimet.mpg.de/cdo)
标题下面的图片:
ind_e = e
ind_T_1500 =温度(下午3点)
ind_es =饱和蒸汽压力
ind_Rh1500 =相对湿度计算出3点。
ind_Tx =最高温度
ind_Tn = Tmin
这是更新后的dataframe(下图)。e是通常略低,但它仍然是不正确的。
T_1500 = Tx - 0.116 * DTR
,DTR是每日温度范围)提供一个温度为3点。e来自Terraclimate平均温度的结果,但这不能让我能想象一个巨大的区别吗? \ endgroup美元