dr_ru.f90 Source File

      subroutine dr_ru
    
      use hydrograph_module
      use hru_lte_module
      use ru_module
      use hru_module, only : ihru, tconc
      
      implicit none

      integer :: ii = 0             !none          |counter
      integer :: ielem = 0          !none          |counter
      real :: rto = 0.              !none          |cloud cover factor
      
     ! compute delivery ratio for each hru in the sub
      do iru = 1, sp_ob%ru
        do ii = 1, ru_def(iru)%num_tot
          ielem = ru_def(iru)%num(ii)
          ihru = ru_elem(ielem)%obtypno

          if (ru_elem(ielem)%dr_name == "calc" .or. ru_elem(ielem)%dr_name == "0") then
            select case (ru_elem(ielem)%obtyp)
            case ("hru")
              rto = tconc(ihru) / ru_tc(iru)
            case ("hlt")
              rto = (hlt_db(ihru)%tc / 3600.) / ru_tc(iru)
            case ("sdc")
              rto = 1.
            case ("ru")
              rto = 1.
            end select
            
            rto = amin1(1.0, rto ** .5)
            ru_elem(ielem)%dr = rto .add. hz
            ru_elem(ielem)%dr%flo = 1.
          end if
          if (ru_elem(ielem)%dr_name == "full" .or. ru_elem(ielem)%dr_name == "1") then
            ru_elem(ielem)%dr = 1. .add. hz
          end if
        end do
      end do
      
      return
      end subroutine dr_ru