7-10月的平均SSTA的全球分布

小宝19950118

我想画多个年份的每年7-10月的SSTA的全球的分布，遇到一点小问题，想了一天没有解决，请教各位

我的ncl脚本：
sst=addfile("E:/NCL/corr/sst.mnmean.v4.nc","r")

;提取sst.mnmean.v4.nc中的经纬度信息
st= sst->sst
lat=sst->lat
lon=sst->lon

lat@units = "degrees_north"
lon@units = "degrees_east"

latS   =  -88.0
latN   =  88.0
lonL   =  0.0
lonR   =  358.0

yrClimStrt = 1979
yrClimLast = 2015

utc_date                = cd_calendar(sst->time, 0)
YYYY                        = tointeger(utc_date(:,0))
istime2                        = ind(YYYY.ge.1979 .and. YYYY.le.1999)
istime3                        = ind(YYYY.ge.2000 .and. YYYY.le.2016)
cltime                        = ind(YYYY.ge.yrClimStrt .and. YYYY.le.yrClimLast)
st!0                        = "time"
st!1                        = "lat"
st!2                        = "lon"
sst_selc2                = st(istime2 , {latS:latN}, {lonL:lonR})
sst_selc3                = st(istime3 , {latS:latN}, {lonL:lonR})
sst_clm                        = st(cltime , {latS:latN}, {lonL:lonR})
slon                        = lon({lonL:lonR})
slat                        = lat({latS:latN})
sstClm                        = clmMonTLL(sst_clm)
ssta2                        = calcMonAnomTLL (sst_selc2, sstClm)
ssta3                        = calcMonAnomTLL (sst_selc3, sstClm)


我已经算出了全球SSTA的月分布（SSTA2和SSTA3），SSTA2表示的是1979-1999的每月SSTA的分布，维度是（时间，经度，维度），要求每年7-10月4个月的平均


我想问两个问题：
1.NCL有没有特定的函数可以解决这个问题？在论坛和ncl官网上找了，但是只有求一个或多个维度的平均（dim_avg_n_Wrap函数），这样的话，只能求全年（1-12月）的平均，不能求特定月份（7-10月）的。


2.如果找不到现成的函数，可以将7-10月的数据提取出来存放在一个新的数组中，再用dim_avg_n_Wrap函数求第一维度的平均。我觉得这个方法可行，但是刚接触ncl，不知道应该如何实现。需不需要写一个循环？


不想做伸手党，实在是不知道怎么办才来发帖子的！求论坛大佬们不吝赐教！

小宝19950118

yeah... 发表于 2017-4-25 08:34
帮你顶上去

谢谢！我已经解决啦！脚本是：
ssta3=new((/68,89,180/),float)

k=0
do i=0,203
do j=0,20
if((i.eq.(12*j+6)).or.(i.eq.(12*j+7)).or.(i.eq.(12*j+8)).or.(i.eq.(12*j+9)))then
ssta3(k,:,:)=ssta33(i,:,:)
k=k+1
end if
end do
end do

yeah...

帮你顶上去

Tobyright

谢谢楼主，刚好遇到同样的问题

HDdd

ncl官网有气候函数可以用