利用比湿和温度计算相对湿度的函数是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)

转载于:https://www.cnblogs.com/yaqiang/p/4651487.html

MeteoInfoLab脚本示例:利用比湿、温度计算相对湿度相关推荐

  1. MeteoInfoLab脚本示例:多Y轴图

    数据范围相差比较大的数据序列进行对比的时候多Y轴图就很重要了.MeteoInfoLab中提供了一个twinx函数来根据已有的坐标系(Axes)生成一个新的Axes,这个命令会使得已有的Axes不绘制右 ...

  2. MeteoInfoLab脚本示例:AIRS Swath HDF数据

    例子中的AIRS Swath HDF数据在Polar Stereographic(南极)投影中接近矩形,需要先从数据中读出经纬度及相关数据数组,利用surfacem函数绘制Swath数据(散点),在s ...

  3. MeteoInfoLab脚本示例:Streamline流线图

    绘制Stramline流线图的函数是streamline,需要两个变量(U/V分量或者风向/风速). 脚本程序: f = addfile('D:/Temp/GrADS/model.ctl') u = ...

  4. MeteoInfoLab脚本示例:TRMM 3B43 HDF数据

    TRMM 3B43是卫星观测月平均降水量产品,是HDF的格点数据.需要注意的是数据中降水变量维的顺序里经度维在前纬度维在后,这与通常的设置(纬度维在前经度维在后)相反,需要对获取的二维数组进行转置,使 ...

  5. MeteoInfoLab脚本示例:OMI Grid HDF数据

    OMI卫星格点数据的例子,全球臭氧柱总量分布. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_L3-OMTO3e_2005m ...

  6. MeteoInfoLab脚本示例:MODIS Sinusoidal投影HDF数据

    MODIS卫星很多陆面数据都是Sinusoidal投影,数据被分为一个个10*10度(赤道地区)的瓦片(http://modis-land.gsfc.nasa.gov/MODLAND_grid.htm ...

  7. MeteoInfoLab脚本示例:站点填图

    打开包含站点填图的站点数据文件(比如micaps 1)之后,用文件对象的smodeldata函数获取StationModel数据对象,然后用stationmodel函数绘制站点填图图层. 脚本程序: ...

  8. MeteoInfoLab脚本示例:获取气团轨迹每个节点的气象数据

    读取HYSPLIT输出的轨迹数据文件和相应时间的气象数据文件,生成轨迹图层,循环每条轨迹的节点,读出该节点的经度.纬度.气压.时间,通过对气象数据插值获得该节点的气象数据. 脚本程序: #------ ...

  9. MeteoInfoLab脚本示例:OMI Swath HDF数据

    这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_ ...

  10. MeteoInfoLab脚本示例:CloudSAT Swath HDF数据

    读取CloudSAT HDF Swath数据,绘图分上下两部分,上面是时间和高度维的Radar Reflectivity Factor二维图,下面是卫星轨迹图. 示例程序: # Add file f ...

最新文章

  1. R语言使用DALEX包的model_performance函数对caret包生成的多个算法模型进行残差分布分析并可视化每个模型的残差反向累积分布图
  2. 【计算机网络】关于数据链路层的讨论(看不懂你来打我!)
  3. linux命令 su和sudo,Linux中sudo和su的区别
  4. 如何分析堆外内存使用情况_堆上与堆外的内存使用情况
  5. python多线程logger_python多线程logger_Dblogger:一个Python语言实现支持多进程和多线程日志打印高效Python日志库......
  6. 是无数像老钟叔的p8u8
  7. 基于JAVA+SpringMVC+Mybatis+MYSQL的线上电器商城
  8. 计算机硬件如何分类,计算机硬件分类教学.doc
  9. CSS综合征病例,ChurgStrauss 综合征 (css) 变应性嗜酸性肉芽肿课件
  10. python出行轨迹记录_Python中的自动路径绘制,python,轨迹
  11. C++语言特性的性能分析
  12. MacBook鼠标指针乱窜/不受控制问题的解决方法
  13. iOS开发之定位神器-超简单方式解决iOS后台定时定位
  14. python提取html中的href标签,如何使用Python从HTML获取href链接?
  15. MySQL Overview
  16. img ,background-img,background-size
  17. h5网页服务器选择,h5制作选择虚拟主机还是服务器
  18. Linux下NFS服务器的配置
  19. 【线性代数】矩阵与线性方程组的几何意义
  20. MS【1】:Metric

热门文章

  1. 给Edge添加chrome主题
  2. Electron构建魔兽世界插件更新安装
  3. MIPS汇编程序设计——Cousera
  4. 赏析角度有哪些_从哪些角度赏析句子
  5. java进行多个照片合并操作
  6. 游戏感:虚拟感觉的游戏设计师指南——第七章 响应的测量方法
  7. SuperMap iClient3D for WebGL制作立体地图
  8. Qt使用libvips处理超大图片(GB以上图片)
  9. Masimo与Penington研究所合作,提升公众对处方阿片类药物过量危害的认识
  10. Google Play的崩溃与ANR