c  written by Vern Kousky 10/29/13
c  Program to compute daily ave. cloud brightness temps. for 1.75 deg boxes
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (imax=40,jmax=40,nobs=2920,ndmax=365)
      dimension x(imax,jmax),xave(imax,jmax,nobs),xdave(imax,jmax)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file=path//'irwin-SA-ave-1.75x1.75deg-boxes-y2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (51,file=path//'irwin-SA-daily-ave-1.75x1.75deg-boxes-y2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      do 1 n=1,nobs
      read(11,rec=n)x
      do 2 i=1,imax
      do 3 j=1,jmax
      xave(i,j,n)=x(i,j)
 3    continue
 2    continue
 1    continue
      min=1
      max=8
      do 4 iday=1,ndmax
      do 5 i=1,imax
      do 6 j=1,jmax
      sum=0.
      cnt=0.
      do 7 n=min,max 
      if(xave(i,j,n).lt.0.)go to 7
      sum=sum+xave(i,j,n)
      cnt=cnt+1.
 7    continue
      if(cnt.eq.0.)go to 77
      xdave(i,j)=sum/cnt
      go to 6
 77   xdave(i,j)=-999.
 6    continue
 5    continue
      write(51,rec=iday)xdave
      min=min+8
      max=max+8
 4    continue
      stop
      end