basin_prm_default.f90 [src] Revision: default  Date:
       subroutine basin_prm_default
    
       use basin_module
       use hru_module, only :  uptake
       
       implicit none
       
       real :: ubw      !                                        |the uptake distribution for water is hardwired
       real :: uobw     !(not used in this model) none           |water uptake normalization parameter 
       real :: uobn     !(not used in this model) none           |nitrogen uptake normalization parameter 
       real :: uobp     !(not used in this model) none           |phosphorus uptake normalization parameter 
      
         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%spcon < 1.e-6) bsn_prm%spcon = 0.0001        !! linear parm for calc sed reentrained in channel sed routin
         if (bsn_prm%spexp < 1.e-6) bsn_prm%spexp = 1.0           !! exponent parameter for calc sed reentrained in channel sed routing
         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 = 20.0        !! 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%trnsrch < 1.e-6) bsn_prm%trnsrch = 0.        !! fraction of transmission losses from main channel that enter
                                                                  !!  deep aquifer
         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%vcrit < 1.e-6) bsn_prm%vcrit = 0.            !!
         if (bsn_prm%res_stlr_co < 1.e-6) bsn_prm%res_stlr_co = 0.184 !! reservoir sediment settling coeff
         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%sig_g < 1.e-6) bsn_prm%sig_g = 0.            !! geometric std dev of part sizes for the main channel

         !! 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)
         
         return
          
      end subroutine basin_prm_default