请大家进来看看帮帮忙，谢谢大家

longlivehj

Maziy 发表于 2014-8-13 16:31
我修改过后还是不行啊，:No coordinate variable exists for dimension (lev) in variable (out)

嗯，应该是有问题的。那个lev取出来是个4维的。
晚一点我用自己的wrfout调试一下，然后回复你。

Maziy

longlivehj 发表于 2014-8-13 16:40
嗯，应该是有问题的。那个lev取出来是个4维的。
晚一点我用自己的wrfout调试一下，然后回复你。

太谢谢你了，十分感谢

longlivehj

Maziy 发表于 2014-8-13 16:44
太谢谢你了，十分感谢

out&lev=lev(0, :, 0, 0)

longlivehj

你的脚本循环太慢了，用dim_max_n_Wrap替换了，你看看，运行效率快很多。

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
f      = addfile("./WRFV3/test/em_quarter_ss/wrfout_d01_0001-01-01_00:00:00.nc","r")
w      = wrf_user_getvar(f,"QCLOUD",-1)
lev    = wrf_user_getvar(f,"pressure",-1)

out = dim_max_n_Wrap(w, (/2, 3/))

out!1="lev"
out&lev=lev(0, :, 0, 0)

;================================================
; plot
;================================================
  wks   = gsn_open_wks ("pdf", "h_time")            ; open ps file
  gsn_define_colormap(wks,"BlWhRe")                ; choose colormap
  res                      = True                  ; plot mods desired  res@cnFillOn             = True                  ; turn on color
  res@cnLinesOn            = True                 ; no contour lines

;---This resource not needed in V6.1.0
  res@gsnSpreadColors      = True                  ; use full range of colors

;---This resource defaults to True in NCL V6.1.0
  res@lbLabelAutoStride    = True                  ; optimal labels

  res@cnLevelSelectionMode = "ManualLevels"       ; set manual contour levels
  res@cnMinLevelValF       = .001                  ; set min contour level
  res@cnMaxLevelValF       = .01                   ; set max contour level
  res@cnLevelSpacingF      = 0.0001                  ; set contour spacing

  plot = gsn_csm_pres_hgt (wks, out(lev|:, Time|:),res)
end

Maziy

longlivehj 发表于 2014-8-13 17:02
你的脚本循环太慢了，用dim_max_n_Wrap替换了，你看看，运行效率快很多。

load "$NCARG_ROOT/lib/ncarg/ ...

还是不行啊，
fatal:dim_max_n: Invalid dimension sizes to take maximum across, can't continue
fatal:["Execute.c":7556]:Execute: Error occurred at or near line 2652 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

fatal:["Execute.c":7556]:Execute: Error occurred at or near line 17 in file test2.ncl

Maziy

longlivehj 发表于 2014-8-13 17:02
你的脚本循环太慢了，用dim_max_n_Wrap替换了，你看看，运行效率快很多。

load "$NCARG_ROOT/lib/ncarg/ ...

可以了，我那个out = dim_max_n_Wrap(w, (/2, 3/))写错了，这个函数具体后面这个/2,3/是什么意思啊

Maziy

Maziy 发表于 2014-8-13 20:23
还是不行啊，
fatal:dim_max_n: Invalid dimension sizes to take maximum across, can't continue
fat ...

然后图上这样的，但是为什么中间数据会缺一块呢

longlivehj

Maziy 发表于 2014-8-13 20:31
然后图上这样的，但是为什么中间数据会缺一块呢

看了一下，我这边没有这样的情况。你根据白线出现的横坐标位置，检查一下数据吧！

longlivehj

Maziy 发表于 2014-8-13 20:27
可以了，我那个out = dim_max_n_Wrap(w, (/2, 3/))写错了，这个函数具体后面这个/2,3/是什么意思啊

(/2, 3/)的意思是对数组在第3、4维上求最大值。

Maziy

longlivehj 发表于 2014-8-13 21:34
(/2, 3/)的意思是对数组在第3、4维上求最大值。

麻烦你了，学到好多东西，十分感谢