MeteoInfoLab脚本示例：利用比湿、温度计算相对湿度

利用比湿和温度计算相对湿度的函数是qair2rh(qair, temp, press=1013.25)，三个参数分别是比湿、温度和气压，气压有一个缺省值1013.25，因此计算地面相对湿度的时候也可以不给气压参数。下面的例子打开CFS模式的地面结果，并从比湿和温度数组计算出相对湿度数组并绘图。

脚本程序：

#Open and read data
fn = ‘U:/data/climate/CFS/20150712/flxf2015081500.01.2015071200.grb2‘
f = addfile(fn)
temp = f[‘Temperature_surface‘][0,::-1,:]
sh = f[‘Specific_humidity_height_above_ground‘][0,0,::-1,:]
#Calculate relative humidity
temp = temp - 273.15
rh = qair2rh(sh, temp)
#Plot
axesm()
mlayer = shaperead(‘D:/Temp/map/country1.shp‘)
geoshow(mlayer, edgecolor=(0,0,255))
#layer = imshowm(rh, 20)
layer = contourfm(rh)
title(‘Relative humidity‘)
colorbar(layer, aspect=12)