Author : stefan.liersch@pik-potsdam.de Date : 2010-02-25 modified: 2010-02-26
PURPOSE : Assigning individual subbasin parameters read from subcatch.prm
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | mb |
subroutine catchment_assign_subcatch(mb)
!-------------------------------------------------------------------------------
! Author : stefan.liersch@pik-potsdam.de
! Date : 2010-02-25
! modified: 2010-02-26
!
! PURPOSE : Assigning individual subbasin parameters read from subcatch.prm
!
! CALLED : from subroutine read_subcatch_def
!-------------------------------------------------------------------------------
use groundwater, only : &
gw_abf, &
gw_bff, &
gw_delay, &
gw_rchrgc, &
gw_revapc, &
gw_revapmn, &
abf, &
bff, &
delay, &
gwht, &
gwq, &
syld, &
rchrgc, &
revapc, &
revapmn
use soil, only: bsn_cncor, bsn_sccor, cncor, sccor
use river, only: bsn_roc2, bsn_roc4, roc2, roc4
use snow, only: &
bsn_gmrate, &
bsn_smrate, &
bsn_tmelt, &
bsn_tsnfall, &
smrate, &
gmrate, &
tmelt, &
tsnfall
use evapotranspiration, only: bsn_thc, bsn_ecal, thc, ecal
integer, intent(in) :: mb
integer :: i, si
do i = 1, mb
si = subcatch_id(i)
if (si /= 0) then
! groundwater parameters
bff(i) = gw_bff(subcatch_id(i))
gwht(i) = .5 ! gw_gwht(subcatch_id(i))
gwq(i) = .5 ! gw_gwq(subcatch_id(i))
abf(i) = gw_abf(subcatch_id(i))
syld(i) = .003 ! gw_syld(subcatch_id(i))
delay(i) = gw_delay(subcatch_id(i))
revapc(i) = gw_revapc(subcatch_id(i))
rchrgc(i) = gw_rchrgc(subcatch_id(i))
revapmn(i) = gw_revapmn(subcatch_id(i))
! bsn parameters
ecal(i) = bsn_ecal(subcatch_id(i))
thc(i) = bsn_thc(subcatch_id(i))
sccor(i) = bsn_sccor(subcatch_id(i))
roc2(i) = bsn_roc2(subcatch_id(i))
roc4(i) = bsn_roc4(subcatch_id(i))
cncor(i) = bsn_cncor(subcatch_id(i))
tsnfall(i) = bsn_tsnfall(subcatch_id(i))
tmelt(i) = bsn_tmelt(subcatch_id(i))
smrate(i) = bsn_smrate(subcatch_id(i))
gmrate(i) = bsn_gmrate(subcatch_id(i))
end if
end do
end subroutine catchment_assign_subcatch