addresses #41 & gets rid of N limitation of wood decay to struc in MI…

…MICS
wwieder · Nov 23, 2020 · 165716f · 165716f
1 parent e4ea7f3
commit 165716f
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 8 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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), &