NCL绘制色斑图-多要素叠加

原创

于 2020-03-25 21:06:16 发布 · 8.3k 阅读

33 ·

CC 4.0 BY-SA版权

本文介绍了如何下载GFS气象数据，并利用NCL（NCAR Command Language）来绘制不同高度层（500hPa至950hPa）的风、温度和湿度色斑图。通过示例代码展示了具体步骤，最终呈现了三个层次的图形结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

需求：下载GFS数据，绘制不同高度层（500hPa~950hPa）风温湿度；

GFS数据下载地址：https://www.cnblogs.com/tiandi/p/11045602.html

下载下来样例数据格式如：gfs.t06z.pgrb2.0p25.f000.grb2

代码如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
;---Open file and read data
  filename = "/cygdrive/d/learn/tonghang/gfs.t06z.pgrb2.0p25.f000.grb2"
  a    = addfile(filename,"r")
  lon = a->lon_0({100:125})
  lat = a->lat_0({15:40})
  lv=a->lv_ISBL0
  HGT= a->HGT_P0_L100_GLL0(:,{15:40},{100:125}) 
  u = a->UGRD_P0_L100_GLL0(:,{15:40},{100:125})
  v = a->VGRD_P0_L100_GLL0(:,{15:40},{100:125})
  rh = a->RH_P0_L100_GLL0(:,{15:40},{100:125})

;---Start the graphics
do i = 0,11
  wks = gsn_open_wks("png","/cygdrive/d/learn/tonghang/975-500wind-height-rh/wind&height&rh"+lv(i)/100)
  gsn_define_colormap(wks,"MPL_Blues")
;---Set some re