c program to summarize the length and average temperature

c  of individual wetting periods

       dimension prd(500),pday(500),phr(500),ptemp(500),prain(500),
     +   pwet(500,4), wetsen(3)

       logical waswet
       THRESH = 0.2



      data i/0/,waswet/.false./,ltemp,lrh,lrain/0.,0.,0./

      data temp,rh,wet,rain/4*0.0/



      open(unit=5,file='wet1wpuc.prn')

      open(unit=6,file='wp4uc.txt', STATUS='UNKNOWN')
      write(6,*) 'This had better work'
      write(*,*) 'line 10'
      write(6,*) 'Day, hr, temp, hrs, rain '
c10    read(5,*,end=999) day,hr,temp,x,rh,xtemp,x2,xrh,wet,xwet,rain
     
10    read(5,*,end=999) day,hr,temp,rh, rain, wetsen(1), 
     +     wetsen(2), wetsen(3)
c      write(6,*) day, hr, temp, rh, wet, rain
c50    format(5f3.0,f5.0)
      wet = (wetsen(1) + wetsen(2) + wetsen(3))/3 
      if(wet.lt.0.0) wet=0.0

      if(rain.lt.0.0) rain=0.0

      if(temp.lt.-50.) temp=ltemp

      if(rh.lt.-50.) rh=lrh

      if(rain.lt.-50.) rain=lrain

c      wet=wet/100.

C     temp=(temp-32.)*0.5555

C     rain=rain*2.54/100.

      ltemp=temp

      lrh=rh

      lrain=rain



      if((wet.gt.THRESH).and.(.not.waswet)) then

        i=i+1

        prd(i)=0

        ptemp(i)=0

        stemp=0

        prain(i)=0

        pday(i)=day

        phr(i)=hr
        DO 25, j=1,3
        pwet(i,j) = 0    
25     continue 
      endif

      

      if (wet.gt.THRESH) then

        prd(i)=prd(i)+1

        stemp=stemp+temp

        prain(i)=prain(i)+rain
        Do 35, j=1,3
c        if (wet.gt.0.2) then
         pwet(i,j) = wetsen(j) + pwet(i,j)
c        else
c        pwet(i,j) = 0 + pwet(i,j)
35      continue
        
        waswet=.true.
        write(*,*) waswet
      endif

c      write(6,*) wet

      if((wet.le.THRESH).and.(waswet)) then
c        write(6,*) 'readme'
        waswet=.false.
c        write(6,*) 'readme'
        ptemp(i)=stemp/prd(i)
        write(6,60) 'red',pday(i),phr(i),ptemp(i),prd(i),prain(i), 
     +  pwet(i,1), pwet(i,2), pwet(i,3)
c        write(6,*) 'wetting period'
60      format(a4,1x,f5.0,1x,f5.0,1x,f6.2,1x,f6.0,1x,f6.2,1x,f4.1,
     +  1x,f4.1,1x,f4.1)
c        write(6,*) 'readme'
      endif



      goto 10



999   close(unit=5)

      close(unit=6)

      stop

      end