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

MeteoInfo

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

脚本程序：

1. #Open and read data
   
2. fn = 'U:/data/climate/CFS/20150712/flxf2015081500.01.2015071200.grb2'
   
3. f = addfile(fn)
   
4. temp = f['Temperature_surface'][0,::-1,:]
   
5. sh = f['Specific_humidity_height_above_ground'][0,0,::-1,:]
   
6. #Calculate relative humidity
   
7. temp = temp - 273.15
   
8. rh = meteo.qair2rh(sh, temp)
   
9. #Plot
   
10. axesm()
    
11. geoshow('country', edgecolor=(0,0,255))
    
12. #layer = imshowm(rh, 20)
    
13. layer = contourfm(rh)
    
14. title('Relative humidity')
    
15. colorbar(layer, aspect=12)

#Open and read data <br /> fn = 'U:/data/climate/CFS/20150712/flxf2015081500.01.2015071200.grb2' <br /> f = addfile(fn) <br /> temp = f['Temperature_surface'][0,::-1,:] <br /> sh = f['Specific_humidity_height_above_ground'][0,0,::-1,:] <br /> #Calculate relative humidity <br /> temp = temp - 273.15 <br /> rh = meteo.qair2rh(sh, temp) <br /> #Plot <br /> axesm() <br /> geoshow('country', edgecolor=(0,0,255)) <br /> #layer = imshowm(rh, 20) <br /> layer = contourfm(rh) <br /> title('Relative humidity') <br /> colorbar(layer, aspect=12)

wsjqxjy

{:eb502:}{:eb502:}

Hhhhleo

老师 每次都出现name 'qair2rh' is not defined的错误，想问下怎么解决

MeteoInfo

Hhhhleo 发表于 2018-5-13 23:15
老师 每次都出现name 'qair2rh' is not defined的错误，想问下怎么解决

改为： meteo.qair2rh

Hhhhleo

老师，我看了很多脚本，但是还是不是特别懂黄色区域的内容应该怎么理解。

MeteoInfo

Hhhhleo 发表于 2018-5-14 22:45
老师，我看了很多脚本，但是还是不是特别懂黄色区域的内容应该怎么理解。

仔细看看这里：http://www.meteothinker.com/docs/meteoinfolab/user_guide.html

213de88

老师 假如我要是用这个几个物理量计算q比湿，应该如何写公式？最后如何存txt呢