output_store_hydrotope_values Subroutine

  subroutine output_store_hydrotope_values(id, array)
    ! Store hydrotope values at the end of timestep.
    integer, intent(in) :: id
    real, dimension(:, :), intent(in) :: array
    integer s, ns, nh

    if (id > 0) then
      ! Unsparse hydrotope data
      do s = 1, output_nsubbasins
        ns = output_subbasin_hydrotope_stix(s)
        nh = output_subbasin_nhydrotopes(s)
        output_storage_hydrotope(id, ns:ns+nh-1, output_store_day_counter) = array(s, :nh)
      end do
    end if
  end subroutine output_store_hydrotope_values