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,ndmax=900)
      dimension x(imax,jmax),xvar(imax,jmax),cnt(imax,jmax)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file=path//'irwin-SA-daily-anoms-boxes-nov1980-2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (51,file=path//'irwin-SA-daily-anoms-var-boxes-nov1980-2009'
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      do 1 i=1,imax
      do 2 j=1,jmax
      xvar(i,j)=0.
      cnt(i,j)=0.
 2    continue
 1    continue
      do 4 iday=1,ndmax
      read(11,rec=iday)x
      do 5 i=1,imax
      do 6 j=1,jmax
      if(x(i,j).lt.0.)go to 6
      cnt(i,j)=cnt(i,j)+1.
      xvar(i,j)=xvar(i,j)+x(i,j)*x(i,j)
 6    continue
 5    continue
 4    continue
      do 7 i=1,imax
      do 8 j=1,jmax
      xvar(i,j)=xvar(i,j)/cnt(i,j)
 8    continue
 7    continue
      write(51,rec=1)xvar
      stop
      end