pro nino34,nyr,ts,n34

ilon1 = 76
ilon2 = 96
ilat1 = 45
ilat2 = 50

nmnth = nyr*12
tr34 = fltarr(12,nyr)
for iyr=0,nyr-1 do begin
   for imn=0,11 do begin
      tr34(imn,iyr) = mean(ts(ilon1:ilon2,ilat1:ilat2,imn+iyr*12))
   endfor
endfor


tr34clm = fltarr(12)

for i=0,11 do begin
   tr34clm(i) = total(tr34(i,*))/nyr
endfor

tr34anom=fltarr(nmnth)

for iyr=0,nyr-1 do begin
   for imn=0,11 do begin
      tr34anom(imn+iyr*12) = tr34(imn,iyr)-tr34clm(imn)
   endfor
endfor

n34 = tr34anom

for imn=2,nmnth-3 do begin
   n34(imn) = total(tr34anom(imn-2:imn+2))/5.
endfor

return
end