MeteoInfoLab脚本示例：露点和温度计算相对湿度

MeteoInfo

MeteoInfoLab的meteo包包含一些气象常用函数，比如利用露点和温度计算相对湿度的函数：dewpoint2rh，参数是露点和温度。MICAPS第1类数据中包含了露点和温度，下面的例子给出了计算相对湿度的代码。

1. fn = 'D:/Temp/micaps/10101420.000'
   
2. f = addfile_micaps(fn)
   
3. lon = f['Longitude'][:]
   
4. lat = f['Latitude'][:]
   
5. temp = f['Temperature'][:]
   
6. dp = f['DewPoint'][:]
   
7. rh = meteo.dewpoint2rh(dp, temp)
   
8. 
   
9. #Plot
   
10. axesm()
    
11. geoshow('country', edgecolor='k')
    
12. layer = scatterm(lon, lat, rh, edge=False, zorder=0)
    
13. colorbar(layer)
    
14. title(u'相对湿度', fontname=u'黑体')

fn = 'D:/Temp/micaps/10101420.000'<br /> f = addfile_micaps(fn)<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> temp = f['Temperature'][:]<br /> dp = f['DewPoint'][:]<br /> rh = meteo.dewpoint2rh(dp, temp)<br /> <br /> #Plot<br /> axesm()<br /> geoshow('country', edgecolor='k')<br /> layer = scatterm(lon, lat, rh, edge=False, zorder=0)<br /> colorbar(layer)<br /> title(u'相对湿度', fontname=u'黑体')

guo066600

王老师，我用CRU中的vapour pressure数据计算露点温度的代码如下：
fn = "F:/CRU/vap/cru_ts4.06.2001.2010.vap.dat.nc"
f = addfile(fn)
vap = f['vap']
x = vap[:,:,:]
T_dew = meteo.dewpoint(x)
运行结果中有很多NaN，查看原始的vap也有很多NaN

请问王老师，这是数据异常吗？

MeteoInfo

guo066600 发表于 2022-10-4 09:26
王老师，我用CRU中的vapour pressure数据计算露点温度的代码如下：
fn = "F:/CRU/vap/cru_ts4.06.2001.201 ...

meteo.dewpoint(x)