c     Program written 9/6/10 by Vern Kousky
c     program to calculate the average daily precip
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (imax=157,jmax=226,maxd=1380)
      dimension x(imax,jmax),sum1(imax,jmax),cnt1(imax,jmax)
      character path*34
      data path/'/cpc/home/wd52vk/samer/precip/jas/'/
      open (11,file=path//'oi-t382-daily-precip-SA-jas1996-2010',
     *access='direct',form='unformatted',recl=imax*jmax*4)
      open (51,file=path//'oi-t382-avedaily-precip-SA-jas1996-2010',
     *access='direct',form='unformatted',recl=imax*jmax*4)
      do 99 i=1,imax
      do 98 j=1,jmax
      sum1(i,j)=0.
      cnt1(i,j)=0.
 98   continue
 99   continue
      do 998 iday=1,maxd
      read(11,rec=iday)x
      do 1 i=1,imax
      do 2 j=1,jmax
      if(x(i,j).lt.0.)go to 2
      sum1(i,j)=sum1(i,j)+x(i,j)
      cnt1(i,j)=cnt1(i,j)+1.
 2    continue
 1    continue
 998  continue
      do 3 i=1,imax
      do 4 j=1,jmax
      if(cnt1(i,j).eq.0.)go to 5
      sum1(i,j)=sum1(i,j)/cnt1(i,j)
      go to 4
 5    sum1(i,j)=-9.99e+08
 4    continue
 3    continue
      write(6,200)sum1(74,113),cnt1(74,113)
 200  format(2f10.1)
      write(51,rec=1)sum1
      stop
      end