一:安装软件和准备数据

1.需要安装Vapor(注意安装路径不要存在空格)

注:版本2.4.2及以后

2.安装NCL,方法见http://www.cnblogs.com/striver-zhu/p/4574819.html

注:

  a.得包含这两个可执行文件:convert && psplit,一般安装时会自行安装进去

  b. NCL_5.1.0版本存在一个问题:wrf2geotiff.ncl中的492行的str_lower需要去掉,否则运行时报错

3.下载数据:jamingWRF2

4.3个ncl文件(NCAR\VAPOR\share\examples\NCL目录下):
  wrf_Height.ncl (this plots humidity, temperature, pressure and wind at a fixed elevation)
  wrf_Precip.ncl (this plots precipitation tendency with pressure isobars)
  wrf_crossSection2.ncl (this is a vertical plot of relative humidity)

二:数据转换
1.将WRF数据转换成VDC数据,从Vapor2.0开始,便不需要转换,但是大数据还是转换以下比较好

  通过cmd执行以下命令(进入命令所在目录下):

  wrfvdfcreate:扫描WRF数据,生成用于描述WRF数据的Vapor元数据
    (wrfvdfcreate wrfout_d02_2008-09-28.nc jangmi-09-28.vdf)
  wrf2vdf:将所有的WRF数据中的2D和3D变量转换成Vapor使用的数据格式,即VDC格式
    (wrf2vdf jangmi-09-28.vdf wrfout_d02_2008-09-28.nc   =》 生成的变量会存在jangmi-09-28_data目录下)

三:
1.将NCL绘制的数据图形转换成georeferenced格式(在工作目录下,需要添加.hluresfile文件,该文件只能通过cygwin的vi建立(:wq退出保存))

  wrf_Height.ncl:生成在一个特定平面上的图形()
  wrf_Precip.ncl:产生一个呈现降水趋势和海平面压力的图形
  wrf_CrossSection2.ncl:产生温度和相对湿度的垂直图形,沿着一个特定的平面上,平行于XZ轴
  

  我们将图形转换成.PS输出文件和将这些.PS文件转换成地理参考TIFF文件

 1).修改脚本,使其仅生成一个图像,在每个时间步长。 (如果想在相同的时间步长,使用多个图像,则需要为每个不同的图像的单独运行脚本)

 2).修改这个脚本,遍历WRF文件中的所有步长。这可能涉及外循环(循环多个文件)和内循环(循环每个文件中的时间步)

 3).修改脚本,输出为.PS文件

 4).将下面的新行到NCL脚本:

  a.在顶部,插入一行来加载wrf2geotiff.ncl:load “wrf2geotiff.ncl”(存在C:\NCAR\VAPOR\share\examples\NCL)
  (这个脚本和其他NCL的例子是默认安装子目录,unix下在share/vapor-x.x.x/examples/NCL,Windows下在$(VAPOR_HOME)/share/examples/NCL/)

  b.在NCL工作站创建之后,调用wrf2geotiff=wrf2geotiff_create(wks)

  c.如果这是一个垂直图形,禁止地理参考:
    wrf2geotiff_disableGeoTags(wrf2gtiff)

  d.设置图形属性res@gsnFrame=False

  e.图形每生成一次,插入两行:
    wrf2geotiff_write()
    frame(wks)

  f.在文件末尾调用wrf2geotiff_close(),将会合并图形生成tiff格式图形

四:获取地理文件

1.为模拟区域获得地形图片:需要经纬度范围,同时需要联网,进入getWMSImage.sh所在文件夹

cygwin下:getWMSImage.sh –o jangmiTerrain.tiff 115 15 130 30

05-08 08:44