求助，初涉ncl绘图和气象数据

striver_zhu

以下信息是用ncl的printVarSummary打印出来的变量t的信息：其中最后两行提示错误，是什么错误啊？？？求解答，非气象专业的小白
谢谢了。。。

绘制等值线的程序，程序源代码：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
fi = addfile("./absolute.nc", "r")
t= fi->tem ; 温度数据
printVarSummary(t)
wks = gsn_open_wks("pdf", get_script_prefix_name())
gsn_define_colormap(wks,"WhBlGrYeRe")
res = True
plot = gsn_csm_contour_map(wks, t({5},:,:), res)
end

打印出来的变量t信息和提示的错误：
Variable: t
Type: short
Total Size: 62208 bytes
            31104 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 12] x [lat | 36] x [lon | 72]
Coordinates:
            time: [1..12]
            lat: [87.5..-87.5]
            lon: [-177.5..177.5]
Number Of Attributes: 5
  long_name :   CRU_Global_1961-1990_Mean_Monthly_Surface_Temperature_Climatology
  units :       celsius
  scale_factor :        0.01
  valid_min :   -9000
  missing_value :       -9999
fatal:Coordinate subscript type mismatch. Subscript (0) can not be coerced to type of coordinate variable
fatal:Execute: Error occurred at or near line 20 in file ztry.ncl

请问怎么解决啊？谢谢各位大神了。。。

longlivehj

这个不是printVarSummary的错误。问题应该处在plot那一行。可能的原因是，time坐标变量的类型不是整形。把该行大括号里面的5改成5.0试试看。

唐屠夫

学习了   

striver_zhu

longlivehj 发表于 2014-7-8 16:00
这个不是printVarSummary的错误。问题应该处在plot那一行。可能的原因是，time坐标变量的类型不是整形。把 ...

我也不知道，我把那个大括号去掉了，结果就好了

longlivehj

striver_zhu 发表于 2014-7-8 18:35
我也不知道，我把那个大括号去掉了，结果就好了

呵呵，好吧，加与不加大括号是有区别的，好好学习一下吧！

striver_zhu

longlivehj 发表于 2014-7-9 12:06
呵呵，好吧，加与不加大括号是有区别的，好好学习一下吧！

问一下子啊，那个加与不加的区别在哪里啊？小白完全不明白

longlivehj

striver_zhu 发表于 2014-7-9 12:28
问一下子啊，那个加与不加的区别在哪里啊？小白完全不明白

http://www.ncl.ucar.edu/Document ... es.shtml#Subscripts

striver_zhu

longlivehj 发表于 2014-7-9 14:00
http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#Subscripts

XIEXIENI....

BDRUSH

longlivehj 发表于 2014-7-9 14:00
http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#Subscripts

大神就是厉害