Rev 587

* Output hh data after NYSKIP only
* Fix accidental implicit variable use

Author: mlt

src/main.f

@@ -55,11 +55,11 @@ program main
 
       use parm
       implicit none
-      prog = "SWAT Jan 28 2013    VER 2012/Rev 585"
+      prog = "SWAT Feb 12 2013    VER 2012/Rev 587"
 
       write (*,1000)
  1000 format(1x,"               SWAT2012               ",/,             &
-     &          "               Rev. 585               ",/,             &
+     &          "               Rev. 587               ",/,             &
      &          "      Soil & Water Assessment Tool    ",/,             &
      &          "               PC Version             ",/,             &
      &          " Program reading from file.cio . . . executing",/)

src/print_hyd.f

@@ -25,7 +25,7 @@ subroutine print_hyd
       use parm
 
 !!    mauro/jerry whittaker hourly output file
-      if (iphr > 0) then 
+      if (iphr > 0 .and. curyr > nyskip) then 
          do ij = 1, nstep
           write (83,1000) iyr,i,ij,ihout,hhvaroute(2,ihout,ij)
          end do 
@@ -34,5 +34,4 @@ subroutine print_hyd
 !!    end hourly codes
 
       return
-      end
-
+      end

src/readrte.f

@@ -260,22 +260,22 @@ subroutine readrte
 !!  An estimate of channel bank erodibility coefficient from jet test if it is not available
 !!  Units of kd is (cm^3/N/s)
 !!  Base on Hanson and Simon, 2001
-      if (ch_bnk_kd(i) <= 1.e-6) then
-	  if (tc_bnk(i) <= 1.e-6) then
-	    ch_bnk_kd(i) = 0.2
+      if (ch_bnk_kd(irch) <= 1.e-6) then
+	  if (tc_bnk(irch) <= 1.e-6) then
+	    ch_bnk_kd(irch) = 0.2
 	  else 
-          ch_bnk_kd(i) = 0.2 / sqrt(tc_bnk(i))
+          ch_bnk_kd(irch) = 0.2 / sqrt(tc_bnk(irch))
 	  end if
 	end if
 
 !!  An estimate of channel bed erodibility coefficient from jet test if it is not available
 !!  Units of kd is (cm^3/N/s)
 !!  Base on Hanson and Simon, 2001
-      if (ch_bed_kd(i) <= 1.e-6) then
-	  if (tc_bed(i) <= 1.e-6) then
-	    ch_bed_kd(i) = 0.2
+      if (ch_bed_kd(irch) <= 1.e-6) then
+	  if (tc_bed(irch) <= 1.e-6) then
+	    ch_bed_kd(irch) = 0.2
 	  else 
-          ch_bed_kd(i) = 0.2 / sqrt(tc_bed(i))
+          ch_bed_kd(irch) = 0.2 / sqrt(tc_bed(irch))
 	  end if
 	end if