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=365)
      dimension x(imax,jmax),xcli(imax,jmax),xa(imax,jmax)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file=path//'irwin-SA-daily-ave-1.75x1.75deg-boxes-y2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (12,file=path//'irwin-daily-clim365-boxes-y1980-2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (51,file=path//'irwin-SA-daily-anoms-boxes-y2009'
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      do 4 iday=1,ndmax
      read(11,rec=iday)x
      read(12,rec=iday)xcli
      do 5 i=1,imax
      do 6 j=1,jmax
      if(x(i,j).lt.0.)go to 60
      if(x(i,j).gt.310.)go to 60
      xa(i,j)=x(i,j)-xcli(i,j)
      go to 6
 60   xa(i,j)=-999.
 6    continue
 5    continue
      write(51,rec=iday)xa
 4    continue
      stop
      end