pro temp

;IMPORTANT: run bin-ex.pro before temp.pro, 
   ;.run bin-ex.pro, then .run temp.pro

;converting ut times by first selecting just numbers after the decimal
;point, and then making them into whole numbers respective to the hour
;they represent (first hr=1, last hr=0)

ut=data.time
for k=0,23 do ut(k)=(ut(k)-floor(ut(k)))*24
      

 ;print, lon
 
;converting ut and lon into local times because want local times wrt
;to the earth than times wrt to the sun; next, looping through j=0-23 ut
;times and producing as many loct's as lon's

sloct=0
for j=0,23 do begin
loct= ut(j) + lon/15
 loct= loct mod(24)   ;override repeated hrs produced by ut wrt to 144 lon
  print, "ut=" , ut(j)
 loct= shift(loct, 6*(j+1)) 
  if loct(0) gt 0.0 then loct(0) = 0.0
   print, loct   ;24= 0th hr, 0= 1st hr

 sloct=sloct + loct  ;sum all local times

endfor

  aveloct= sloct/24  ;creates monthly averaged local time bin for day
   print, "aveloct=" , aveloct

;create new merwinds in 3D from the 4D array mer=(*,*,*,h) h=ut(0,23)

 smer=0
for k=0,23 do begin
 nmer=mer(*,*,*,k)     ;nmer=mer(nlon,nlat,nlev,nut), n=#
 nmer=shift(nmer, 6*(k+1), 0, 0)
 
  smer=smer + nmer
endfor
 avesmer= smer/24
  
;stop 
end