Calculate the Cape Town Procedure TAC

Wrapper for the R/RTMB translation of the Cape Town Procedure (CTP). The function calculates the gene-tagging and CPUE components, fits the CKMR-driven adult population model, and applies the CTP harvest control rule. The internal CKMR optimization must return convergence code zero, a finite objective, and a finite final gradient whose maximum absolute component does not exceed 5e-4. Its exact final RTMB Hessian must also be finite, symmetric, and positive definite, with every eigenvalue exceeding sqrt(.Machine$double.eps).

Usage

CTP(ctpobj)

Arguments

ctpobj

A list containing the CTP inputs: mpdat, gtN, gtR, cpue, cpuedat, popdat, hspdat, paly, paly_first_yr, tacold, first_yr, current_yr, and implementation_yr.

Value

A list containing the bounded updated TAC, the constrained deltaTAC, independently auditable GT, CPUE, CKMR, and unbounded-TAC intermediate values, optimizer convergence code and objective, max_abs_gradient, gradient_tolerance, gradient_refinements, the final optimizer route and any deterministic numerical nudge, Hessian eigenvalue and condition-number diagnostics, and fitted parameter vector par. If the final zero or maximum-TAC bound is active, deltaTAC need not equal TAC - ctpobj$tacold.