Fixed temperature mortality stress function to have the right type of…

… arguments.
NGEET · May 30, 2018 · a36a79d · a36a79d
1 parent 0b6261d
commit a36a79d
Showing 1 changed file with 5 additions and 6 deletions.
diff --git a/biogeochem/EDMortalityFunctionsMod.F90 b/biogeochem/EDMortalityFunctionsMod.F90
@@ -60,10 +60,9 @@ subroutine mortality_rates( cohort_in,bc_in,cmort,hmort,bmort,frmort )
     real(r8) :: frac  ! relativised stored carbohydrate
     real(r8) :: b_leaf ! target leaf biomass kgC
     real(r8) :: hf_sm_threshold    ! hydraulic failure soil moisture threshold 
-    real(r8) :: temp_dep           ! Temp. function (freezing mortality)
+    real(r8) :: temp_dep_fraction  ! Temp. function (freezing mortality)
     real(r8) :: temp_in_C          ! Daily averaged temperature in Celcius
-    real(r8),parameter :: frost_mort_buffer = 5.0_r8  ! 5deg buffer for freezing mortality
-
+    real(r8), parameter :: frost_mort_buffer = 5.0_r8  ! 5deg buffer for freezing mortality
     logical, parameter :: test_zero_mortality = .false. ! Developer test which
                                                         ! may help to debug carbon imbalances
                                                         ! and the like
@@ -107,9 +106,9 @@ subroutine mortality_rates( cohort_in,bc_in,cmort,hmort,bmort,frmort )
     !           doi: 10.1111/j.1365-2486.2006.01254.x                                    
 
     temp_in_C = bc_in%t_veg24_si - tfrz
-    temp_dep  = max(0.0,min(1.0,1.0 - (temp_in_C - &
-                EDPftvarcon_inst%freezetol(cohort_in%pft))/frost_mort_buffer) )
-    frmort    = EDPftvarcon_inst%mort_scalar_coldstress(cohort_in%pft) * temp_dep
+    temp_dep_fraction  = max(0.0_r8, min(1.0_r8, 1.0_r8 - (temp_in_C - &
+                         EDPftvarcon_inst%freezetol(cohort_in%pft))/frost_mort_buffer) )
+    frmort    = EDPftvarcon_inst%mort_scalar_coldstress(cohort_in%pft) * temp_dep_fraction
 
 
     !mortality_rates = bmort + hmort + cmort