ncl求助

sxing

gg=addfile("hgt.mon.mean.nc","r")
ymStrt=195001
ymLast=201512
yrStrt=ymStrt/100
yrLast=ymLast/100
tt=gg->time
TIME=ut_calendar(tt,-1)/100
iTIME=ind(TIME.ge.yrStrt.and.TIME.le.yrLast)
g2=gg->hgt(iTIME,:,{18:55},{73:135})
;printVarSummary(g2)
lev1=gg->level
print(lev1)
lat1=gg->lat
lon1=gg->lon

gh=g2(0:65,:,:,:)

  do i=0,779,12
    gh(floattointeger(i/12),:,:,:)=(g2(i+6,:,:,:)+g2(i+6,:,:,:)+g2(i+6,:,:,:))/3;  7月份
  end do
  copy_VarMeta(g2(0:65,:,:,:),gh)
printVarSummary(sph)
;====================== =======
hyear=(/1981,2000,2009,2012,2013/)            
lyear=(/1980,1983,1992,1997,2007/)
gh!0="time"
gh&time=ispan(1950,2015,1)
gh!1 = "lev1"
gh!2 = "lat"
gh!3 = "lon"

gh&lat = lat1    ；第41行
gh&lon = lon1
gg2=sph({hyear},5,:,:)
gg3=sph({lyear},5,:,:)

以上是我的部分相关的脚本，标注了出错行数的地方，