From 165716f0875875767f11c0a27aafae693d5bc4be Mon Sep 17 00:00:00 2001 From: wwieder Date: Mon, 23 Nov 2020 11:26:55 -0700 Subject: [PATCH] addresses #41 & gets rid of N limitation of wood decay to struc in MIMICS --- SOURCE_CODE/casa_inout.f90 | 6 ++++-- SOURCE_CODE/mimics_cycle_CN.f90 | 4 ++-- SOURCE_CODE/mimics_inout_CN.f90 | 4 ++-- SOURCE_CODE/mimics_variable_CN.f90 | 4 ++-- 4 files changed, 10 insertions(+), 8 deletions(-) diff --git a/SOURCE_CODE/casa_inout.f90 b/SOURCE_CODE/casa_inout.f90 index 99238288..ebb46240 100644 --- a/SOURCE_CODE/casa_inout.f90 +++ b/SOURCE_CODE/casa_inout.f90 @@ -1617,8 +1617,10 @@ SUBROUTINE biogeochem(iYrCnt,idoy,mdaily,nppScalar,cleaf2met,cleaf2str,croot2met ! Added isomModel argument to the casa_xkN2 function call. -mdh 3/23/2020. !!call casa_xkN2(xkNlimiting,casapool,casaflux,casamet,casabiome,veg) call casa_xkN2(xkNlimiting,casapool,casaflux,casamet,casabiome,veg,isomModel) - - DO j=1,mlitter + ! is this where wood decomp is slowed by xkNlimiting, even for MIMICS? + ! instead of looping through each litter, we'll just do this for leaves, roots + ! should only be done for MIMICS, as here, see also line 1565? - ww 11/23/20 + DO j=1,2 !mlitter casaflux%klitter(:,j) = casaflux%klitter(:,j)* xkNlimiting(:) ! write(*,*) 'biogeochem:' ! write(*,'(a24,3(f6.4,2x))') 'xkNlimiting(1) = ', xkNlimiting(1) diff --git a/SOURCE_CODE/mimics_cycle_CN.f90 b/SOURCE_CODE/mimics_cycle_CN.f90 index 92b37174..913eb9dd 100644 --- a/SOURCE_CODE/mimics_cycle_CN.f90 +++ b/SOURCE_CODE/mimics_cycle_CN.f90 @@ -1591,8 +1591,8 @@ SUBROUTINE mimics_soil_reverseMM_CN(mp,iYrCnt,idoy,mdaily,cleaf2met,cleaf2str,cr ! 0.3 tau_r(2) ! 0.00024 tau_k(1) ! 0.1 tau_k(2) - mimicsbiome%CN_r(npt) = mimicsbiome%CNr * SQRT(mimicsbiome%fmet(npt)/mimicsbiome%cnModDenom) - mimicsbiome%CN_k(npt) = mimicsbiome%CNk * SQRT(mimicsbiome%fmet(npt)/mimicsbiome%cnModDenom) + mimicsbiome%CN_r(npt) = mimicsbiome%CNr * SQRT(mimicsbiome%cnModNum/mimicsbiome%fmet(npt)) + mimicsbiome%CN_k(npt) = mimicsbiome%CNk * SQRT(mimicsbiome%cnModNum/mimicsbiome%fmet(npt)) mimicsbiome%tauR(npt) = mimicsbiome%tau_r(1) * & exp(mimicsbiome%tau_r(2) * mimicsbiome%fmet(npt)) * mimicsbiome%tauMod(npt) diff --git a/SOURCE_CODE/mimics_inout_CN.f90 b/SOURCE_CODE/mimics_inout_CN.f90 index 06a38698..ca96327c 100644 --- a/SOURCE_CODE/mimics_inout_CN.f90 +++ b/SOURCE_CODE/mimics_inout_CN.f90 @@ -306,12 +306,12 @@ SUBROUTINE mimics_readbiome(fname_mimicsbiome, mp, mvtype) read(101,*) mimicsbiome%CNr read(101,*) mimicsbiome%CNk read(101,*) mimicsbiome%fracDINavailMIC - read(101,*) mimicsbiome%cnModDenom + read(101,*) mimicsbiome%cnModNum write(*,'(2x,a5,2x,f10.6)') 'CNr=', mimicsbiome%CNr write(*,'(2x,a5,2x,f10.6)') 'CNk=', mimicsbiome%CNk write(*,'(2x,a16,2x,f10.6)') 'fracDINavailMIC=', mimicsbiome%fracDINavailMIC - write(*,'(2x,a13,2x,f10.6)') 'cnModDenom=', mimicsbiome%cnModDenom + write(*,'(2x,a13,2x,f10.6)') 'cnModNum=', mimicsbiome%cnModNum !! ---------------------------------------------------------------------- diff --git a/SOURCE_CODE/mimics_variable_CN.f90 b/SOURCE_CODE/mimics_variable_CN.f90 index 7abd9bdb..ab6491ed 100644 --- a/SOURCE_CODE/mimics_variable_CN.f90 +++ b/SOURCE_CODE/mimics_variable_CN.f90 @@ -67,7 +67,7 @@ MODULE mimicsvariable REAL(r_2), POINTER :: densDep, & fracDINavailMIC, & - cnModDenom, CNr, CNk + cnModNum, CNr, CNk REAL(r_2), DIMENSION(:), POINTER :: fmet, & ligninNratioAvg, & @@ -206,7 +206,7 @@ SUBROUTINE alloc_mimicsvariable(mp,mvtype,mplant) mimicsbiome%densDep, & mimicsbiome%CNr, & mimicsbiome%CNk, & - mimicsbiome%cnModDenom, & + mimicsbiome%cnModNum, & mimicsbiome%fracDINavailMIC) ALLOCATE(mimicsbiome%Vint(nRK), &