NCL模式数据处理目录

• wrfout文件处理和绘图
• • 文件处理基础
  • 文件处理进阶
  • • 插值
    • 平滑
• MODIS数据介绍
• • 文件名
• .shp数据介绍

wrfout文件处理和绘图

文件处理基础

• 在不打开NCL的情况下，用ncl_filedump或者ncdump查看文件信息
• 同时也可以用wrf_user_getvar查看变量信息，例如：hgt = wrf_user_getvar(f,“HGT”,0)
• 画图时，用到wrf_map_overlays（在WRF-ARW的地图背景上覆盖等高线图和矢量图，例如：plot = wrf_map_overlays(f, wks,(/plt_tem, plt_slp, vec/),
  pltres, mpres)
• 运用.shp文件方法

lnres = True
lnres@gslineColor = "Gray25"
lnres@gslineThicknessF = 3.
usa_id = gsn_add_shapefile_polyline(wks, plot, "xxx.shp", lnres)
draw(plot)
frame(wks) # 绘制地图绘制附加的线

• 其他需要注意
  • res = wrf_map_resources(f, res) ;检查文件，以确定附加到res上的地图绘制资源
  • res@tfDoNDCOverlay = True ； 叠加渲染层，以占据基础图的整个窗口
  • res@gsnAddCyclic = False ; 决定是否要把-180°与180°结合
  • 若想将数据放在不同的地图投影上，需要：tc@lat2d = wrf_user_getvar(f,“lat”,0); tc@lon2d = wrf_user_getvar(f,“lon”,0)
  • getvarfdims(f) – 获得所有的维度的名称
  • getfiledimsizes(f) – 获得对应的维度大小
  • getfilevarnames(f) – 获得文件f内的所有变量名

文件处理进阶

插值

• wrf_user_interp_3d
• wrf_user_interp_level

a = addfile("wrfout_1201.nc","r")
ht = wrf_user_getvar(a,"z",0)
p  = wrf_user_getvar(a,"pressure",0)
ht_500 = wrf_user_interp_level(ht, p, 500, False)

平滑

• wrf_smooth_2d(fld, iter)

MODIS数据介绍

文件名

• MOD–Terra, MYD–Aqua
• D3: daily, E3:every eight days, M3: monthly

.shp数据介绍

• point：城市或景点的位置，人口数据；
• polyline：河流、道路等非封闭的边界；
• polyon：封闭的地理边界，如国界、省界、地区界、湖泊等