求助，，NCL 画WRF结果中的地面长波

无敌泡

这是我的代码，总是第27行有错误，请各位大神指点
;   Example script to plot a field from a single metgrid file
  2 ;   November 2008
  3
  4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  5 load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
  6
  7 begin
  8 ;
  9   a = addfile("wrfout.d01.2011_03_19_00:00:00.nc","r")  ; Open a file
10
11
12 ; We generate plots, but what kind do we prefer?
13   type = "x11"
14 ; type = "pdf"
15 ; type = "ps"
16 ; type = "ncgm"
17   wks = gsn_open_wks(type,"plt_landflux_1")
18   
19
20   res = True                           ; Set up some basic plot resources
21   res@MainTitle = "GROUND LONG WAVELENTH"
22   res@Footer = False
pltres = True
25   mpres = True
26
27   glw    = wrf_user_getvar(a,"GLW",0)
28
29
30
31     opts = res                                ; Set some plotting resources
32     opts@cnFillOn = True
33     opts@ContourParameters = (/ 10. /)
34
35     contour = wrf_contour(a,wks,glw,opts)
36     delete(opts)
37     plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
38
39 end

longlivehj

我执行你的代码，27行没有问题啊！
提示的什么错误？

无敌泡

longlivehj 发表于 2014-4-4 22:09
我执行你的代码，27行没有问题啊！
提示的什么错误？

fatal:Number of dimensions in parameter (2) of (wrf_contour) is (3), (2) dimensions were expected
我这次用的是a->LH读取数据，结果到了 contour = wrf_contour(a,wks,v,opts)，就提示错误阿

longlivehj

无敌泡 发表于 2014-4-5 09:25
fatal:Number of dimensions in parameter (2) of (wrf_contour) is (3), (2) dimensions were expected
...

直接a->LH得到的数组是3维的，当然会出现错误啊！
wrf_user_getvar第三个参数为0，是获取第一个时次的数据，所以得到的数组是2维，是可以出图的。

无敌泡

longlivehj 发表于 2014-4-5 19:09
直接a->LH得到的数组是3维的，当然会出现错误啊！
wrf_user_getvar第三个参数为0，是获取第一个时次的数 ...

谢谢您，现在出来新的错误了，warning:ContourPlotInitialize: scalar field is constant; ContourPlot not possible:[errno=1102]
是说我的数值设置的不对，我看他们全球例子中也是（-300，300，10），这个数值我该怎么确定呢？
谢谢

无敌泡

longlivehj 发表于 2014-4-5 19:09
直接a->LH得到的数组是3维的，当然会出现错误啊！
wrf_user_getvar第三个参数为0，是获取第一个时次的数 ...

就是出来的图，只有底图，，

longlivehj

无敌泡 发表于 2014-4-6 09:07
就是出来的图，只有底图，，

取其它时次看看，就是wrf_user_getvar的第三个参数不要取0，取1，2……看看呢！

无敌泡

longlivehj 发表于 2014-4-6 10:21
取其它时次看看，就是wrf_user_getvar的第三个参数不要取0，取1，2……看看呢！

我太爱你了！！！谢谢！！！

rjys1219

好东西必须顶