新手提问，ncl如何对两个文件中的数据做平均

ql455020268

ncl新手
想问下如何对两个文件中的数据做平均，分别是1号和2号的数据，变量为s1和s2，想求s1和s2的平均，s1，s2都为三维数组（时间，经度，纬度），时间固定（选择第二个“1”），经纬度全选，求大神指导
脚本如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
    f1=addfile("/mnt/hgfs/f/download/grib201306/cdas1.20130601.sfluxgrbf.grb2","r")
    f2=addfile("/mnt/hgfs/f/download/grib201306/cdas1.20130602.sfluxgrbf.grb2","r")
    ;time=f1->time
    s1=f1->SHTFL_P8_L1_GGA0_avg
    s2=f2->SHTFL_P8_L1_GGA0_avg
    aves=(s1(1,:,:)+s2(1,:,:))/2

    wks = gsn_open_wks("png","/mnt/hgfs/f/aves1&2-1")

    res = True ; 建立源变量
    res@mpMinLonF = 90. ; 设定源变量属性 mpMinLonF 指定地图最小经度
    res@mpMaxLonF = 110. ; 设定源变量属性 mpMaxLonF 指定地图最大经度
    res@mpMinLatF = 30. ; 设定源变量属性 mpMinLatF 指定地图最小纬度
    res@mpMaxLatF = 44. ; 设定源变量属性 mpMaxLatF 指定地图最大纬度
    res@cnFillOn=True ;加填充色
    res@lbOrientation="Vertical" ;色标垂直

    ;加地图
    res@mpDataSetName              = "$NCARG_ROOT/lib/ncarg/database/Earth..4"
    res@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"
    res@mpAreaMaskingOn            = True
    res@mpMaskAreaSpecifiers       = (/"China"/)
    res@mpOutlineSpecifiers        = (/"China","China:Provinces"/)

    plot = gsn_csm_contour_map(wks, aves, res)
end