求教高手！！！~NCL画taylor泰勒图运行不灵

leizhang_zl

NCL官网有直接的例子，http://www.ncl.ucar.edu/Applications/taylor.shtml除了taylor_8.ncl 其他的都能画出图，而这个taylor_8.ncl在执行任务后，一直没有动静，也不提示出错，也不会运行结束。向高手求教！

leizhang_zl

;**********************************
; taylor_8.ncl
;**********************************
;
; These files are loaded by default in NCL V6.2.0 and newer
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
;
; This file still has to be loaded manually
load "./taylor_diagram_cam.ncl"
;**********************************
begin
;**********************************
; Assume the following have already been computed:
;   _ratio are the ratio:  Case_Variance/Reference_Variance
;   _cc    are the cross correlation coef of Case to Reference
;   _bias  are the bias estimates
;**********************************

; Cases [Model]
  case      = (/ "Data A", "Data B" /)
  nCase     = dimsizes(case )                 ; # of Cases [Cases]

; variables compared
  var       = (/ "Globe","20S-20N","Land", "Ocean", "N. America", "Africa" /)
  nVar      = dimsizes(var)                   ; # of Variables

; "Case A"                        
  CA_ratio   = (/1.230, 0.988, 1.092, 1.172, 1.064, 0.990 /)
  CA_cc      = (/0.958, 0.973,-0.740, 0.743, 0.922, 0.950 /)
  BA         = (/2.7  ,-1.5  , 17.31,-20.11, 12.5 , 8.341 /)

; "Case B"
  CB_ratio   = (/1.129, 0.996, 1.016, 1.134, 1.023, 0.962 /)
  CB_cc      = (/0.963, 0.975, 0.801, 0.814,-0.946, 0.984 /)
  BB         = (/1.7  , 2.5  ,-17.31, 20.11, 19.5 , 7.341 /)

; arrays to be passed to taylor plot
  ratio      = new ((/nCase, nVar/),typeof(CA_cc) )  
  cc         = new ((/nCase, nVar/),typeof(CA_cc) )
  bias       = new ((/nCase, nVar/),typeof(CA_cc) )

  ratio(0,:) = CA_ratio
  ratio(1,:) = CB_ratio

  cc(0,:)    = CA_cc
  cc(1,:)    = CB_cc

  bias(0,:)  = BA
  bias(1,:)  = BB

;**********************************
; create plot
;**********************************

  res   = True                           ; default taylor diagram
        
  res@Markers      = (/16, 16/)               ; make all solid fill
  res@Colors       = (/"red", "blue" /)         
  res@varLabels    = var
  res@caseLabels   = case

  wks   = gsn_open_wks("png","taylor")        ; send graphics to PNG file
  plot  = taylor_diagram_cam(wks,ratio,cc,bias,res)
end

njyjy

请问楼主现在知道原因了吗我跟你是一样的问题