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