ncl画合成差值场 差值大的地方没有通过检验 反而小的地方通过检验？

jojingwj

请教大牛！
ncl画合成差值场，差值大的地方没有通过检验，反而小的地方通过检验？也用其他数据画过，只有这一个是这样的，程序如下，请问是什么问题呢？
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"

begin
   type = "eps"
  wks = gsn_open_wks(type,"4a-s-c-chazhi-cha")
gsn_define_colormap(wks,"BlueWhiteOrangeRed")
;读取topo
  b        = addfile ("H:\wrf-out\CTRL\c-526", "r")
  gtopo    = wrf_user_getvar(b,"HGT",0)
  lon2d    = wrf_user_getvar(b,"XLONG",0)
  lat2d    = wrf_user_getvar(b,"XLAT",0)
    olon = fspan(64.5,135.5,72)
  olat = fspan(14.5,55.5,42)
  ;data1 = new((/248,152/),"float")
  olon!0          = "lon"
  olon@long_name  = "lon"
  olon@units      = "degrees-east"
  olon&lon        = olon
  olat!0          = "lat"
  olat@long_name  = "lat"
  olat@units      = "degrees_north"
  olat&lat        = olat
   grid_gtopo = rcm2rgrid_Wrap(lat2d,lon2d,gtopo,olat,olon,0)
;创建新数组
varx = new((/5,42,72/),float)
vary = new((/5,42,72/),float)

;读取五组试验  
do i=0,4
  a     = addfile ("H:\wrf-out\CTRL\c-"+(i + 526), "r")
  c     = addfile ("H:\wrf-out\all-tp-big\"+(i + 526), "r")
  times  = wrf_user_getvar(a,"times",-1)  ; get all times in the file
  ntimes  = dimsizes(times )         ; number of times in the file
print(i)   
  rainc0    = wrf_user_getvar(a,"RAINC",-1)
  rainnc0   = wrf_user_getvar(a,"RAINNC",-1)
  rainTot0  = rainc0 + rainnc0         
  x0     =  (rainTot0(ntimes-1,:,:)-rainTot0(ntimes-94,:,:))/92

  grid_x0 = rcm2rgrid_Wrap(lat2d,lon2d,x0,olat,olon,0)
  varx(i,:,:) =grid_x0
   rainc1    = wrf_user_getvar(c,"RAINC",-1)
  rainnc1   = wrf_user_getvar(c,"RAINNC",-1)
  rainTot1  = rainc1 + rainnc1         
  x1     =  (rainTot1(ntimes-1,:,:)-rainTot1(ntimes-94,:,:))/92
  grid_x1 = rcm2rgrid_Wrap(lat2d,lon2d,x1,olat,olon,0)
  vary(i,:,:) =grid_x1
delete(grid_x0)
delete(grid_x1)
delete(times)
delete(ntimes)
delete(rainc0)
delete(rainnc0)
delete(rainTot0)
delete(x0)
delete(rainc1)
delete(rainnc1)
delete(rainTot1)
delete(x1)
end do
xAve = dim_avg_n_Wrap(varx,0)
  yAve = dim_avg_n_Wrap(vary,0)
  print("ok1---------------------------------------")
; Set some Basic Plot options
cha = new((/42,72/),"float")
cha=yAve-xAve
copy_VarCoords(vary(0,:,:),cha(:,:))

xVar = dim_variance_n_Wrap(varx,0)   
yVar = dim_variance_n_Wrap(vary,0)
  print(yVar)
  
n1=5
n2=5
  siglvl = 80
  iflag=True             ; population variance not similar
    prob =100.*(1-ttest(xAve,xVar,n1,yAve,yVar,n2,iflag,False))
     print("ok2---------------------------------------")
   copy_VarCoords(vary(0,:,:),prob(:,:))
         
  
   printVarSummary(prob)
;======================================        

   res=True
   res@gsnAddCyclic = False
   res@gsnMaximize             = True
  res@gsnDraw                 = False
  res@gsnFrame                = False
  
res@gsnLeftStringFontHeightF    = 0.025
res@gsnLeftStringOrthogonalPosF        =       0.05;标题在同一高度上
res@gsnLeftStringFontColor = "black"
res@mpOutlineBoundarySets ="NoBoundaries"

  res@gsnAddCyclic   = False   
  res@cnFillOn    = False
  res@mpDataSetName         = "Earth..4"   ; This new database contains
  res@mpDataBaseVersion      = "MediumRes"  ; Medium resolution database
  res@mpOutlineOn            = True         ; Turn on map outlines
; res@mpOutlineSpecifiers   = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)   

       res@mpMinLatF              =   15       ; Asia limits
  res@mpMaxLatF              = 55
  res@mpMinLonF              =  100
  res@mpMaxLonF              =  140
  res@tmXBValues               = ispan(100,140,10)
  res@tmXBLabels               = ispan(100,140,10)+"~F35~J~F25~E"
  res@tmYLValues               = ispan(15,55,10)
  res@tmYLLabels               = ispan(15,55,10)+"~F35~J~F25~N"
   res@tmXBMinorOn               = False
  res@tmYLMinorOn               = False
;  res@tmYROn                    = False  ; turn off Ylable Graduation
;  res@tmXTOn                    = False  ; turn off Xlable Graduation

;***************************************************************
  res@mpGeophysicalLineThicknessF = 2.      ; double the thickness of geophysical boundaries
  res@mpNationalLineThicknessF = 2.         ; double the thickness of national boundaries

  res@mpAreaMaskingOn = True   ;
  
  res@mpMaskAreaSpecifiers = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)   
  res@mpOceanFillColor = 0     ;
  
  res@mpInlandWaterFillColor = 0  ;
  res@mpLandFillColor = 0  ;"Transparent";=========================="white"

  res@cnFillOn       = True
  res@cnLinesOn      = False          ;绛夊

流浪的心

我也遇到类似问题，使程序错了吗?

Datian

同问啊·······，原因是什么呢》

水长东

同问……然而感觉不会有人

大气科学的花朵

啊啊啊，我也是，怎么解释这种原因呀

jiangbing

接下来的程序有吗

jiangbing

同问啊 ？？我也遇到了这个问题

眼里有星辰

不考虑程序问题的话，那就是差值大的地方本来变率就大呗