-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathbasin_prm_default.f90
63 lines (56 loc) · 5.37 KB
/
basin_prm_default.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
subroutine basin_prm_default
use basin_module
use hru_module, only : uptake
implicit none
if (bsn_prm%evlai < 1.e-6) bsn_prm%evlai = 3.0 !! leaf area index at which no evap occurs
if (bsn_prm%ffcb < 1.e-6) bsn_prm%ffcb = 0. !! intial soil water cont expressed as a fraction of fc
if (bsn_prm%surlag < 1.e-6) bsn_prm%surlag = 4.0 !! surface runoff lag time (days)
if (bsn_prm%adj_pkr < 1.e-6) bsn_prm%adj_pkr = 1.0 !! peak rate adjustment factor in the subbasin
if (bsn_prm%prf < 1.e-6) bsn_prm%prf = 1.0 !! peak rate adjustment factor for sediment routing in the channel
if (bsn_prm%cmn < 1.e-6) bsn_prm%cmn = 0.0003 !! rate factor for mineralization on active org N
if (bsn_prm%n_updis < 1.e-6) bsn_prm%n_updis = 20.0 !! nitrogen uptake dist parm
if (bsn_prm%p_updis < 1.e-6) bsn_prm%p_updis = 20.0 !! phosphorus uptake dist parm
if (bsn_prm%nperco < 1.e-6) bsn_prm%nperco = 0.1 !! nitrate perc coeff (0-1)
!! 0 = conc of nitrate in surface runoff is zero !! 1 = perc has same conc of nitrate as surf runoff
if (bsn_prm%pperco < 1.e-6) bsn_prm%pperco = 10.0 !! phos perc coeff (0-1)
!! 0 = conc of sol P in surf runoff is zero !! 1 = percolate has some conc of sol P as surf runoff
if (bsn_prm%phoskd < 1.e-6) bsn_prm%phoskd = 175.0 !! phos soil partitioning coef
if (bsn_prm%psp < 1.e-6) bsn_prm%psp = 0.40 !! phos availability index
if (bsn_prm%rsdco < 1.e-6) bsn_prm%rsdco = 0.05 !! residue decomposition coeff
if (bsn_prm%percop < 1.e-6) bsn_prm%percop = 0.5 !! pestcide perc coeff (0-1)
if (bsn_prm%msk_co1 < 1.e-6) bsn_prm%msk_co1 = 0.75 !! calibration coeff to control impact of the storage
!! time constant for the reach at bankfull depth
if (bsn_prm%msk_co2 < 1.e-6) bsn_prm%msk_co2 = 0.25 !!
if (bsn_prm%msk_x < 1.e-6) bsn_prm%msk_x = 0.20 !! weighting factor control relative importance of inflow rate
!! and outflow rate in determining storage on reach
if (bsn_prm%nperco_lchtile < 1.e-6) bsn_prm%nperco_lchtile = 0.50 !! n concentration coeff for tile flow and
!! leach from bottom layer
!if (bsn_prm%evrch < 1.e-6) bsn_prm%evrch = 0.60 !! reach evaporation adjustment factor
if (bsn_prm%cdn < 1.e-6) bsn_prm%cdn = 1.40 !! denitrification expoential rate coefficient
if (bsn_prm%sdnco < 1.e-6) bsn_prm%sdnco = 1.30 !! denitrification threshold frac of field cap
if (bsn_prm%bact_swf < 1.e-6) bsn_prm%bact_swf = 0.15 !! frac of manure containing active colony forming units
if (bsn_prm%tb_adj < 1.e-6) bsn_prm%tb_adj = 0. !! adjustment factor for subdaily unit hydrograph basetime
if (bsn_prm%cn_froz < 1.e-6) bsn_prm%cn_froz = 0.000862 !!
!if (bsn_prm%dorm_hr < 1.e-6) bsn_prm%dorm_hr = -1. !! time threshold used to define dormant (hrs)
if (bsn_prm%nfixmx < 1.e-6) bsn_prm%nfixmx = 20.0 !! max daily n-fixation (kg/ha)
if (bsn_prm%decr_min < 1.e-6) bsn_prm%decr_min = 0.01 !!
if (bsn_prm%rsd_covco < 1.e-6) bsn_prm%rsd_covco = 0.30 !! residue cover factor for computing frac of cover
if (bsn_prm%urb_init_abst < 1.e-6) bsn_prm%urb_init_abst = 0. !! PET adjustment (%) for Penman-Montieth and Preiestly-Taylor methods
if (bsn_prm%petco_pmpt < 0.5 .and. bsn_prm%petco_pmpt > 0.) bsn_prm%petco_pmpt = 0.0 !! reservoir sediment settling coeff
bsn_prm%petco_pmpt = (100. + bsn_prm%petco_pmpt) / 100. !! convert to fraction
if (bsn_prm%uhalpha < 1.e-6) bsn_prm%uhalpha = 1.0 !! alpha coeff for est unit hydrograph using gamma func
if (bsn_prm%eros_spl < 1.e-6) bsn_prm%eros_spl = 0. !! coeff of splash erosion varing 0.9-3.1
if (bsn_prm%rill_mult < 1.e-6) bsn_prm%rill_mult = 0. !!
if (bsn_prm%eros_expo < 1.e-6) bsn_prm%eros_expo = 0. !!
if (bsn_prm%c_factor < 1.e-6) bsn_prm%c_factor = 0. !!
if (bsn_prm%ch_d50 < 1.e-6) bsn_prm%ch_d50 = 0. !! median particle diameter of main channel (mm)
if (bsn_prm%co2 < 100.) bsn_prm%co2 = 400. !! assume co2 concentration = 400 ppm
!! set additional parameters
uptake%water_dis = 10.0 !! the uptake distribution for water is hardwired - users are not allowed to modify
uptake%water_norm = 1. - exp(-uptake%water_dis)
uptake%n_norm = 1. - exp(-bsn_prm%n_updis)
uptake%p_norm = 1. - exp(-bsn_prm%p_updis)
!! set maximum days to store hru, ru and channel output
bsn_prm%day_lag_mx = 2
return
end subroutine basin_prm_default