NCL 绘制FNL 位势高度图

hhumuer

NCL 绘制FNL （GRIB2格式）的位势高度图时，显示错误
warning:ContourPlotSetValues: Data values out of range of levels set by MANUALLEVELS mode
不知道脚本错在哪里？还请各位大神出手相助啊~~~~！
脚本如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
  fi  = addfile("./fnl_20140928_00_00.grib2","r")
  hgt=fi->HGT_P0_L100_GLL0(5,:,:)
  wks_type="pdf"
  wks_name="GHT500"
  wks = gsn_open_wks(wks_type,wks_name)
  gsn_define_colormap(wks,"gui_default") ;
  res = True
  res@mpMinLatF = 3
  res@mpMaxLatF = 53
  res@mpMinLonF = 73
  res@mpMaxLonF = 135
  res@cnFillOn = True  ; cn means ContourPlot
  res@cnLinesOn = True
  res@cnLineLabelsOn = True  
  res@cnInfoLabelOn = False ; set false will not set an information label
  res@cnLevelSelectionMode = "ManualLevels"
  res@cnMinLevelValF = 5700
  res@cnMaxLevelValF = 5900
  res@cnLevelSpacingF = 20
  res@gsnSpreadColors = True ; 6.1.0
  res@lbBoxLinesOn = False ; Label bar
  res@lbLabelAutoStride     = True            ; let NCL determine label spacing

  plot = gsn_csm_contour_map(wks, hgt, res)
end

MeteroCl

一楼正解，HGT_P0_L100_GLL0(lv_ISBL0, lat_0, lon_0)变量里面 lv_ISBL0 = 1000, 2000, 3000, 5000, 7000, 10000, 15000, 20000, 25000, 30000,
    35000, 40000, 45000, 50000, 55000, 60000, 65000, 70000, 75000, 80000,
    85000, 90000, 92500, 95000, 97500, 100000 ;高度层500hPa是第14个，5应改为13

尽头的尽头

高度层和你设置的等值线最大最小不对应，500hPa的高度应该不是第六层

海豚∩_∩浩

你用printMinMax（hgt,True）,看看hgt数值最大最小值。

hhumuer

MeteroCl 发表于 2015-10-22 16:20
一楼正解，HGT_P0_L100_GLL0(lv_ISBL0, lat_0, lon_0)变量里面 lv_ISBL0 = 1000, 2000, 3000, 5000, 7000,  ...

历害！！！我在官网上没有看到介绍lv_ISBL0或者其他维度的介绍信息，请问大侠，在哪里可以找到说明文档呢？

hhumuer

海豚∩_∩浩 发表于 2015-10-22 14:14
你用printMinMax（hgt,True）,看看hgt数值最大最小值。

嗯，这个很有用！谢谢啦！！！

hhumuer

尽头的尽头 发表于 2015-10-22 13:16
高度层和你设置的等值线最大最小不对应，500hPa的高度应该不是第六层

大神，正解！太谢谢了，请问哪里可以下载到FNL各变量不同层的说明文档之类的不？

MeteroCl

hhumuer 发表于 2015-10-22 16:33
历害！！！我在官网上没有看到介绍lv_ISBL0或者其他维度的介绍信息，请问大侠，在哪里可以找到说明文档呢 ...

我用ncl_filedump *.grib2 >& grb2.log 命令生成的说明文档

hhumuer

MeteroCl 发表于 2015-10-22 16:41
我用ncl_filedump *.grib2 >& grb2.log 命令生成的说明文档

这个输出得到的，只有关于lv_ISBL0的简单介绍，比如一共26层，单位是pa，但没有说具体每一层是多少pa呢

MeteroCl

hhumuer 发表于 2015-10-22 17:08
这个输出得到的，只有关于lv_ISBL0的简单介绍，比如一共26层，单位是pa，但没有说具体每一层是多少pa呢

还有一个办法，就是你把grib格式转成nc格式，然后ncdump来看就有详细信息啦