non negative sBranch

Author: ForModLabUHel

src/C_regionPrebas.f90

@@ -260,15 +260,15 @@ subroutine regionPrebas(siteOrder,HarvLim,minDharv,multiOut,nSites,areas,nClimID
 	 species = int(multiOut(siteX,ij,4,ijj,1))
 	 multiWood(siteX,ij,ijj,1) = multiWood(siteX,ij,ijj,2) / pCrobas(2,species)
 	 energyWood = energyWood + multiWood(siteX,ij,ijj,1) * areas(siteX)   !!energCuts !!!we are looking at volumes
-	 multiOut(siteX,ij,28,ijj,1) = (multiOut(siteX,ij,24,ijj,1)+multiOut(siteX,ij,51,ijj,1)) * &
+	 multiOut(siteX,ij,28,ijj,1) = max(0.,((multiOut(siteX,ij,24,ijj,1)+multiOut(siteX,ij,51,ijj,1)) * &
 		(1-energyRatio) + multiOut(siteX,ij,28,ijj,1) + multiOut(siteX,ij,32,ijj,1)*0.83 +&
-		 multiOut(siteX,ij,31,ijj,1)* (1-harvRatio) * (1-energyRatio)
+		 multiOut(siteX,ij,31,ijj,1)* (1-harvRatio) * (1-energyRatio)))
      multiOut(siteX,ij,29,ijj,1) = multiOut(siteX,ij,32,ijj,1)*0.17*(1-energyRatio)+ &
 			multiOut(siteX,ij,29,ijj,1) !0.1 takes into account of the stem residuals after clearcuts
 	else
-	 multiOut(siteX,ij,28,ijj,1) = multiOut(siteX,ij,24,ijj,1) + multiOut(siteX,ij,28,ijj,1) + &
+	 multiOut(siteX,ij,28,ijj,1) = max(0.,(multiOut(siteX,ij,24,ijj,1) + multiOut(siteX,ij,28,ijj,1) + &
 		multiOut(siteX,ij,24,ijj,1) + multiOut(siteX,ij,32,ijj,1)*0.83 + &
-		multiOut(siteX,ij,31,ijj,1)* (1-harvRatio)
+		multiOut(siteX,ij,31,ijj,1)* (1-harvRatio)))
      multiOut(siteX,ij,29,ijj,1)=multiOut(siteX,ij,32,ijj,1)*0.17+multiOut(siteX,ij,29,ijj,1) !0.1 takes into account of the stem residuals after clearcuts
 	endif
 !!energCuts