OmegaThreeFisheries/runSizespectrum/YieldCalc.R

In nissandjac/commspectrum: Community size spectrum simulator

YieldCalc <- function(param,S){
  
  # Calculate the equilibrium yield from a size-spectrum run 
  
  tEnd<- param$tEnd/param$dt
  
  N <- S$N[tEnd,,]
  Yield <- matrix(0,param$nSpecies)  
  for (i in 1:param$nSpecies){
    
    Yield[i] <- sum(N[i,]*S$w*S$dw*S$Fin[i,])
    
  }
  
  return(Yield)
}