doDepletion: doDepletion - depletes the stock to the declared initdep...

View source: R/aspm.r

doDepletionR Documentation

doDepletion - depletes the stock to the declared initdep level

Description

doDepletion - depletes the stock to the declared initdep level There is a printout to the screen of the final outcome. The procedure assumes you have the unfished NaA in column 1 of the stock$NaA object. The depletion is arrived at by literally searching for the first constant harvest rate that leads to the required depletion.

Usage

doDepletion(inR0, indepl, inprops, inglb, inc = 0.02, Numyrs = 50)

Arguments

inR0

the estimated log transformed unfished recruitment

indepl

the initial depletion level to be searched for

inprops

the props object from the data object

inglb

the globals object from the data object

inc

the starting value and step in harvest rate used when finding the selected depletion level; defaults to 0.02

Numyrs

the number of years of fishing while searching for the desired depletion. Some number between 40 - 50 seems to work well.

Value

A list containing the vector of numbers-at-age at the required depletion level, the Fvalue required to achieve that depletion, and the actual depletion level achieved

Examples

## Not run: 
data(dataspm)
fish <- dataspm$fish
glb <- dataspm$glb
props <- dataspm$props
dep <- doDepletion(glb$R0,indepl=0.4,props,glb,inc=0.02)
print(dep)

## End(Not run)  

haddonm/datalowSA documentation built on Nov. 5, 2023, 6:40 p.m.