@@ -874,18 +874,27 @@ subroutine fuse_cohorts(patchptr, bc_in)
874874 nextc% n* nextc% lmort_infra)/ newn
875875
876876 ! npp diagnostics
877- currentCohort% npp_leaf = (currentCohort% n* currentCohort% npp_leaf + nextc% n* nextc% npp_leaf)/ newn
878- currentCohort% npp_froot = (currentCohort% n* currentCohort% npp_froot + nextc% n* nextc% npp_froot)/ newn
879- currentCohort% npp_bsw = (currentCohort% n* currentCohort% npp_bsw + nextc% n* nextc% npp_bsw)/ newn
880- currentCohort% npp_bdead = (currentCohort% n* currentCohort% npp_bdead + nextc% n* nextc% npp_bdead)/ newn
881- currentCohort% npp_bseed = (currentCohort% n* currentCohort% npp_bseed + nextc% n* nextc% npp_bseed)/ newn
882- currentCohort% npp_store = (currentCohort% n* currentCohort% npp_store + nextc% n* nextc% npp_store)/ newn
877+ currentCohort% npp_leaf = (currentCohort% n* currentCohort% npp_leaf + nextc% n* nextc% npp_leaf) &
878+ / newn
879+ currentCohort% npp_froot = (currentCohort% n* currentCohort% npp_froot + nextc% n* nextc% npp_froot) &
880+ / newn
881+ currentCohort% npp_bsw = (currentCohort% n* currentCohort% npp_bsw + nextc% n* nextc% npp_bsw) &
882+ / newn
883+ currentCohort% npp_bdead = (currentCohort% n* currentCohort% npp_bdead + nextc% n* nextc% npp_bdead) &
884+ / newn
885+ currentCohort% npp_bseed = (currentCohort% n* currentCohort% npp_bseed + nextc% n* nextc% npp_bseed) &
886+ / newn
887+ currentCohort% npp_store = (currentCohort% n* currentCohort% npp_store + nextc% n* nextc% npp_store) &
888+ / newn
883889
884890 ! biomass and dbh tendencies
885891 currentCohort% ddbhdt = (currentCohort% n* currentCohort% ddbhdt + nextc% n* nextc% ddbhdt)/ newn
886- currentCohort% dbalivedt = (currentCohort% n* currentCohort% dbalivedt + nextc% n* nextc% dbalivedt)/ newn
887- currentCohort% dbdeaddt = (currentCohort% n* currentCohort% dbdeaddt + nextc% n* nextc% dbdeaddt)/ newn
888- currentCohort% dbstoredt = (currentCohort% n* currentCohort% dbstoredt + nextc% n* nextc% dbstoredt)/ newn
892+ currentCohort% dbalivedt = (currentCohort% n* currentCohort% dbalivedt + nextc% n* nextc% dbalivedt) &
893+ / newn
894+ currentCohort% dbdeaddt = (currentCohort% n* currentCohort% dbdeaddt + nextc% n* nextc% dbdeaddt) &
895+ / newn
896+ currentCohort% dbstoredt = (currentCohort% n* currentCohort% dbstoredt + nextc% n* nextc% dbstoredt) &
897+ / newn
889898
890899 do i= 1 , nlevleaf
891900 if (currentCohort% year_net_uptake(i) == 999._r8 .or. nextc% year_net_uptake(i) == 999._r8 ) then
0 commit comments