catch: Catch at age

In vquennessen/densityratio: Using Density Ratios to Help Manage Fisheries In And Around Marine Reserves

Catch at age

Description

catch returns the catch at age vector. It is based on the continuous (Baranov) formulation of of selectivity, such that fishing mortality (FM) and natural mortality (M) are proportional to fish abundance (N) and act simultaneously and uniformly throughout the year, i.e. dN/dt = -(M+F)*N.

Usage

catch(t, cr, fdr, FM, Nat_mortality, N, A, Fb, E, Catch)

Arguments

t	temporary numeric value, the current time step.
cr	temporary numeric value, the current control rule.
fdr	temporary numeric value, the current final target density ratio.
FM	numeric array that corresponds to the fishing mortality at each age in each area, at each timestep, under all control rules, with all estimates of natural mortality.
Nat_mortality	numeric vector, the estimates of natural mortality.
N	numeric array, the number of individuals at each age, in each area, at each timestep, under each control rule, and for each estimate of natural mortality.
A	numeric value, the number of total areas in the model. Default value is 5.
Fb	numeric value, the historical fishing effort for the fished species.
E	numeric array, the relative fishing effort displayed in each area, at each time step, under each control rule, and for each natural mortality estimate.
Catch	numeric array, the number of individuals caught at each age, in each area, at each timestep, under each control rule, for each estimate of natural mortality.

Value

a numeric array with an updated vector of catch at ages given the current area, timestep, and control rule.

Examples

n = 34; A = 5; TimeT = 70; CR = 6; FDR = 4
FM <- array(rep(0.2, n*A*TimeT*CR*FDR), c(n, A, TimeT, CR, FDR))
N <- array(rep(10, n*A*TimeT*CR*FDR), c(n, A, TimeT, CR, FDR))
E <- array(rep(1, A*TimeT*CR*FDR), c(A, TimeT, CR, FDR))
Catch <- array(rep(2, n*A*TimeT*CR*FDR), c(n, A, TimeT, CR, FDR))
catch(t = 1, cr = 1, fdr = 1, FM, Nat_mortality = c(0.14, 0.09, 0.19),
   N, A = 5, Fb = 0.2, E, Catch)

vquennessen/densityratio documentation built on Aug. 28, 2022, 5:36 p.m.