update ECMmodel

Author: ForModLabUHel

src/A_routines.f90

@@ -1994,10 +1994,10 @@ subroutine rhoMcalc(rho_M, CN,pars)
 
 
 ! subroutine CUEcalc(ETS, st,r_r,W_RT,r_RT,rm_aut_roots,litt_RT,exud,normFactP,normFactETS,P_RT,pars)
-subroutine CUEcalc(ETS, st,r_r,W_RT,r_RT,rm_aut_roots,litt_RT,exud,normFactP,P_RT,pars)
+subroutine CUEcalc(ETS, st,r_r,W_RT,r_RT,rm_aut_roots,litt_RT,exud,P_RT,pars)
 
 implicit none
-	real(8),intent(in) :: ETS, st,r_r,W_RT, normFactP, pars(12)!, normFactETS
+	real(8),intent(in) :: ETS, st,r_r,W_RT, pars(12)!, normFactETS
 	real(8),intent(out) :: rm_aut_roots,litt_RT,exud,r_RT,P_RT
 
 	real(8) :: r_F, s_F, rho_M,CN
@@ -2028,7 +2028,7 @@ subroutine CUEcalc(ETS, st,r_r,W_RT,r_RT,rm_aut_roots,litt_RT,exud,normFactP,P_R
 
 ! Here, r_RT is the apparent maintenance respiration rate of fine roots when C input to the fungi has been taken into account.
 !Fine root maintenance respiration is replaced with the following:
-r_RT = (r_r+r_F * rho_M)/(1+rho_M) * normFactP
+r_RT = (r_r+r_F * rho_M)/(1+rho_M)
 
 ! When calculating the related loss of C to the atmosphere, we need to subtract the losses from this, i.e., use 
 !rm_aut_roots  losses of C from roots & ECM to 
@@ -2039,14 +2039,14 @@ subroutine CUEcalc(ETS, st,r_r,W_RT,r_RT,rm_aut_roots,litt_RT,exud,normFactP,P_R
 !normalise parameters according to P and ETS factors
 ! s_H = s_H * normFactETS
 ! phi_M = phi_M * normFactETS
-rm_aut_roots = rm_aut_roots * normFactP
+rm_aut_roots = rm_aut_roots
 
 ! To Yasso as root litter (or later maybe special composition as hyphal litter)
-litt_RT = (rho_M * h_M * s_H)/(1+rho_M) * W_RT * normFactP
+litt_RT = (rho_M * h_M * s_H)/(1+rho_M) * W_RT
 ! To Yasso as sugar (W) exudates
-exud = (rho_M *ksi_M)/(1+rho_M) * W_RT * normFactP
+exud = (rho_M *ksi_M)/(1+rho_M) * W_RT
 
-P_RT =  (rho_M *ksi_M*phi_M)/(1+rho_M) * W_RT * normFactP
+P_RT =  (rho_M *ksi_M*phi_M)/(1+rho_M) * W_RT
 
 end subroutine
 

src/B_prebas.f90

@@ -111,7 +111,7 @@ subroutine prebas(nYears,nLayers,nSp,siteInfo,pCrobas,initVar,thinning,output, &
 real (kind=8) :: Nmort, BAmort
 !!ECMmodelling
  real (kind=8) :: r_RT, rm_aut_roots, litt_RT, exud(nLayers), P_RT
- real (kind=8) :: normFactP, Cost_m !normFactETS,!!Cost_m is the "apparent maintenance respiration" rate of fine roots when C input to the fungi has been taken into account.
+ real (kind=8) :: Cost_m !normFactP,normFactETS,!!Cost_m is the "apparent maintenance respiration" rate of fine roots when C input to the fungi has been taken into account.
 
 !fix parameters
  real (kind=8) :: qcTOT0,Atot,fAPARprel(365)
@@ -728,9 +728,9 @@ subroutine prebas(nYears,nLayers,nSp,siteInfo,pCrobas,initVar,thinning,output, &
 			! ECM modelling
 			if(ECMmod==1) then !!!ECMmodelling
 			    ! normFactETS = 1. + par_aETS * (ETS-ETS_ref)/ETS_ref
-				normFactP = p0 / p0_ref
+				! normFactP = p0 / p0_ref
 				! call CUEcalc(ETS, sitetype,par_mr0,W_froot,r_RT,rm_aut_roots,litt_RT,exud(ij),normFactP,normFactETS,P_RT,pECMmod) !!!ECMmodelling
-				call CUEcalc(ETS, sitetype,par_mr0,W_froot,r_RT,rm_aut_roots,litt_RT,exud(ij),normFactP,P_RT,pECMmod) !!!ECMmodelling
+				call CUEcalc(ETS, sitetype,par_mr,W_froot,r_RT,rm_aut_roots,litt_RT,exud(ij),P_RT,pECMmod) !!!ECMmodelling
 				modOut((year+1),45,ij,1) = P_RT  !add priming to heterotrophic respiration
 				Respi_m = par_mf * wf_STKG + par_mw * W_wsap + rm_aut_roots * W_froot  !!!ECMmodelling
 				Cost_m  = par_mf * wf_STKG + par_mw * W_wsap + r_RT * W_froot  !!!ECMmodelling