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