标签:
利用比湿和温度计算相对湿度的函数是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)
MeteoInfoLab脚本示例:利用比湿、温度计算相对湿度
标签:
原文地址:http://www.cnblogs.com/yaqiang/p/4651487.html