求助：用最基础的gsn_csm_contour_map画图，仍然画不出来

ireneshixiao

总报错：contour = gsn_csm_contour(wks,sm1,res)fatal:ContourPlotDraw: CPCLAM/AREDAM - UNABLE TO INSERT NODE INTO MAP
但是我的目标数据就是二维的，该定义的空值也定义了，求帮忙！！
printVarSummary(sm1)

Variable: sm1
Type: double
Total Size: 6600 bytes
            825 values
Number of Dimensions: 2
Dimensions and sizes:   [latitude | 25] x [longitude | 33]
Coordinates:
            latitude: [  36..  30]
            longitude: [ 115.. 123]
Number Of Attributes: 4
  long_name :   Variable
  units :       m^3.m^-3*100
  missing_value :       -9999
  _FillValue :  -9999



具体代码就这几行：
fn = "try.nc"
fi = addfile(fn, "r")
lon=fi->longitude
lat=fi->latitude

sm1=fi->sm
lat@units ="degrees_north"
lon@units ="degrees_east"

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

wks = gsn_open_wks("png","copy")
res                      = True               ; plot mods desired
res@tiMainString         = "Default Color"    ; main title
res@cnFillOn             = True               ; turn on color fill
res@gsnAddCyclic =True
res@cnFillPalette       = "gui_default"

plot = gsn_csm_contour_map(wks,sm1,res)


求帮忙啊！要画的图就是sm的空间分布，sm我之前已经在matlab处理成只有Longitude和latitude二维的数据了，时间维度已去除掉。但就是画不出

rolemodel

把sm1=fi->sm换成sm=fi->sm1，把最后一行改成plot = gsn_csm_contour_map(wks,sm,res)，试试

ireneshixiao

rolemodel 发表于 2017-10-4 07:47
把sm1=fi->sm换成sm=fi->sm1，把最后一行改成plot = gsn_csm_contour_map(wks,sm,res)，试试

已经找出原因啦 是原数据中的NAN在Ncl中无法识别 要用Ncl的函数把NaN替换一下呢 多谢你啦