画中国两个时期夏季降水差值图，打点区域是对的，可是差值填色图却很奇怪

一颗梅子糖

    dy_ind = ind(year.le.1993.and.year.ge.1979)
    ppre2 =f->pre(dy_ind,:,:)
    opt = True
    opt@nval_crit = 12   ; require at least 12 values for the "avg" is calculated.            
    pre_mon = calculate_monthly_values (ppre2, "avg", 0,opt)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    dy_ind2 = ind(year.le.2008.and.year.ge.1994)
    pp2 =f->pre(dy_ind2,:,:)
    opt = True
    opt@nval_crit = 12   ; require at least 12 values for the "avg" is calculated.            
    ppre_mon = calculate_monthly_values (pp2, "avg", 0,opt)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

pre1 = new((/15,163,283/), float)
do i = 1994,2008
pre1(i-1994,:,:) = (pre_mon(5+(i-1994)*12,:,:)+pre_mon(6+(i-1994)*12,:,:)+pre_mon(7+(i-1994)*12,:,:))/3
end do
copy_VarMeta(pre_mon(0,:,:),pre1(0,:,:))
pre1!0 = "year"
pre1&year = ispan(1994, 2008, 1)
comp1 = dim_avg_n_Wrap(pre1(:,:,:), 0)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
pre2 = new((/15,163,283/), float)
do i = 1994,2008
pre2(i-1994,:,:) = (ppre_mon(5+(i-1994)*12,:,:)+ppre_mon(6+(i-1994)*12,:,:)+ppre_mon(7+(i-1994)*12,:,:))/3
end do
copy_VarMeta(ppre_mon(0,:,:),pre2(0,:,:))
pre2!0 = "year"
pre2&year = ispan(1994, 2008, 1)
comp2 = dim_avg_n_Wrap(pre2(:,:,:), 0)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
comp = comp1 - comp2
copy_VarMeta(comp1, comp)

massage

pre1和pre2有没有经纬度属性？我看着应该是没有的

来自苹果APP客户端

一颗梅子糖

massage 发表于 2022-7-5 23:58
pre1和pre2有没有经纬度属性？我看着应该是没有的

谢谢你啊，画出来了