nhr=24
nday=52
wi_zp2_d=fltarr(nx,ny,nhr,nday)
ui_zp2_d=fltarr(nx,ny,nhr,nday)
nmf2_d=fltarr(nx,ny,nhr,nday)
hmf2_d=fltarr(nx,ny,nhr,nday)
tec_d=fltarr(nx,ny,nhr,nday)

for iday=0,51 do begin
   wi_zp2_d(*,*,0:23,iday) = wi_zp2(*,*,iday*24:(iday+1)*24-1)
   ui_zp2_d(*,*,0:23,iday) = ui_zp2(*,*,iday*24:(iday+1)*24-1)
   nmf2_d(*,*,0:23,iday) = nmf2(*,*,iday*24:(iday+1)*24-1)
   hmf2_d(*,*,0:23,iday) = hmf2(*,*,iday*24:(iday+1)*24-1)
   tec_d(*,*,0:23,iday) = tec(*,*,iday*24:(iday+1)*24-1)
endfor

wiav = fltarr(nx,ny,nhr)
uiav = fltarr(nx,ny,nhr)
nmf2av = fltarr(nx,ny,nhr)
hmf2av = fltarr(nx,ny,nhr)
tecav = fltarr(nx,ny,nhr)

wisdv = fltarr(nx,ny,nhr)
uisdv = fltarr(nx,ny,nhr)
nmf2sdv = fltarr(nx,ny,nhr)
hmf2sdv = fltarr(nx,ny,nhr)
tecsdv = fltarr(nx,ny,nhr)

for it=0,23 do begin
   for j=0,ny-1 do begin
      for i=0,nx-1 do begin
         wiav(i,j,it) = mean(reform(wi_zp2_d(i,j,it,*)))
         wisdv(i,j,it) = stddev(reform(wi_zp2_d(i,j,it,*)))
         uiav(i,j,it) = mean(reform(ui_zp2_d(i,j,it,*)))
         uisdv(i,j,it) = stddev(reform(ui_zp2_d(i,j,it,*)))
         nmf2av(i,j,it) = mean(reform(nmf2_d(i,j,it,*)))
         nmf2sdv(i,j,it) = stddev(reform(nmf2_d(i,j,it,*)))
         hmf2av(i,j,it) = mean(reform(hmf2_d(i,j,it,*)))
         hmf2sdv(i,j,it) = stddev(reform(hmf2_d(i,j,it,*)))
         tecav(i,j,it) = mean(reform(tec_d(i,j,it,*)))
         tecsdv(i,j,it) = stddev(reform(tec_d(i,j,it,*)))
      endfor
   endfor
endfor

end