求取异常值（时次非12的倍数）求助！！！

aa22556688

对全球海冰（时间：1870.01~2015.02（dimsizes（time）=1742），多出两个月），calcMonAnomTLL要求必须是12的倍数：
所以先利用calcMonAnomTLL求出前1740个时次的异常值；
然后利用dim_avg_n求出第1741，1742个时次的异常值；
最后把他们放到一个ssta（1741，：，：）数组里面。但是画出来的图有问题。

主要在于copy_VarMeta(st,ssta) 的所放的位置。
当copy_VarMeta(st,ssta)放在下面第一个位置的时候，
画ssta（1739：1741，：：），图依次为下图：qian-1739（正常），qian-1740（有问题），qian-1741（有问题）。

当copy_VarMeta(st,ssta)放在下面第二个位置的时候，
画ssta（1739：1741，：：），图依次为下图：hou-1739（有问题），hou-1740（正常），hou-1741（正常）。

求高手指点，遇到这种时间不是12倍数的数据，例如2015年多出几个月，但想求其全部时次的异常值，应该如何求，NCL初学者，求指点！！！
latS             =  -89.5
latN             =  89.5
lonL             =  -179.5                 
lonR            =  179.5
yrStrt           = 1870
yrLast           = 2014
yrClimStrt      = 1961
yrClimLast     = 1990
dir    = "/cvfs01/disk1/wangq316/zc/d/"
fi1    = "HadISST_ice.nc"
fsst   = addfile(dir+fi1,"r")
st     = fsst->sic
ss     = fsst->sic(0,:,:)
lat    = fsst->latitude
lon    = fsst->longitude
time   = fsst->time
utc_date  = cd_calendar(time, 0)
YYYY       = tointeger(utc_date(:,0))  
istime      = ind(YYYY.ge.yrStrt.and.YYYY.le.yrLast)
cltime      = ind(YYYY.ge.yrClimStrt.and.YYYY.le.yrClimLast)
sst_selc   = st(istime,{latS:latN},{lonL:lonR})  
sst_clm    = st(cltime,{latS:latN},{lonL:lonR})
sstClm     = clmMonTLL(sst_clm)
ssta        = new((/dimsizes(time),dimsizes(lat),dimsizes(lon)/),"float",st@_FillValue)
copy_VarMeta(st,ssta)  ;;; 第一个位置
ssta(0:1739,:,:)     = calcMonAnomTLL(sst_selc(0:1739,:,:),sstClm)
ssta(1740,:,:)        = st(1740,:,:)-dim_avg_n(st(1092:1440:12,:,:),0)
ssta(1741,:,:)        = st(1741,:,:)-dim_avg_n(st(1093:1441:12,:,:),0)
copy_VarMeta(st,ssta)  ;;; 第二个位置
wks = gsn_open_wks("ps" ,"polar")  
;gsn_define_colormap(wks,"wgne15")            
res                                    = True                          
res@gsnPolar                     = "NH"   
res@cnFillOn                      = False          ; turns on the color
res@mpMinLatF                  = 40   
res@gsnSpreadColors          = True                ; use full colormap
res@gsnSpreadColorStart     = 2
res@gsnSpreadColorEnd      = -3        
plot = gsn_csm_contour_map_polar(wks,ssta(1739,:,:),res)

longlivehj

用缺测值补成12的倍数

aa22556688

longlivehj 发表于 2015-5-13 20:35
用缺测值补成12的倍数

哦哦，对对这确实是一个不错的方法，非常感谢！