HarvestControlRule Class Reference

A harvest control rule class. More...

#include <HarvestControlRule.h>

Collaboration diagram for HarvestControlRule:

Public Types
enum	enumHCR {   FORTY_TEN, FIXED_ESCAPEMENT, FIXED_ESCAPEMENT_CAP, FIXED_HARVEST_RATE,   CONDITIONAL_CONSTANT_CATCH, THIRTY_TWENTY, FIXED_HR_DELTA, FLOOR_THIRTY_TWENTY,   FAO_PA_COMPLIANT }

Public Member Functions
	HarvestControlRule (int &eHCR)
double	getTac (const double &bt, const double &fmsy, const double &msy, const double &bmsy, const double &bo, const double &delta, const double &ptac, const double &mintac)
	Return the TAC given estimates of current biomass and reference points. More...
double	FortyTen (const double &bt, const double &bo, const double &fmsy)
	Implement the 40:10 harvest control rule. More...
double	ThirtyTwenty (const double &bt, const double &bo, const double &fmsy)
	30:20 harvest control rule More...
double	FixedHarvestRate (const double &bt, const double &fmsy)
double	FixedEscapement (const double &bt, const double &bmsy)
double	FixedEscapementCap (const double &bt, const double &bmsy, const double &msy)
double	ConditionalConstantCatch (const double &bt, const double &bmsy, const double &msy, const double &fmsy)
	Implement the Conditional Constant Catch harvest control rule proposed by the IPHC. More...
double	FixedHarvestRateDelta (const double &bt, const double &fmsy, const double &delta, const double &ptac)
double	FloorThirtyTwenty (const double &bt, const double &bo, const double &fmsy, const double &mintac)

Detailed Description

A harvest control rule class.

Author

Steven Martell

Remarks

Contains an enum type that defines the various harvest control rules that can be used

Definition at line 10 of file HarvestControlRule.h.

Member Enumeration Documentation

enum HarvestControlRule::enumHCR

A enum type to select the appropriate harvest control rule.

Enumerator
FORTY_TEN	40:10 harvest contrl rule
FIXED_ESCAPEMENT	Fixed escapement rule
FIXED_ESCAPEMENT_CAP	Fixed escapement rule with a cap at MSY
FIXED_HARVEST_RATE	Fixed harvest rate rule based on Fmys
CONDITIONAL_CONSTANT_CATCH	IPHC's Conditional constant catch rule
THIRTY_TWENTY	30:20 harvest control rule
FIXED_HR_DELTA	Fixed catch with 15% maximum changed
FLOOR_THIRTY_TWENTY	30:20 rule with a minimum catch floor.
FAO_PA_COMPLIANT	Not implemented yet

Definition at line 19 of file HarvestControlRule.h.

        {
                FORTY_TEN,                      
                FIXED_ESCAPEMENT,                       
                FIXED_ESCAPEMENT_CAP,           
                FIXED_HARVEST_RATE,                     
                CONDITIONAL_CONSTANT_CATCH,     
                THIRTY_TWENTY,                          
                FIXED_HR_DELTA,             
                FLOOR_THIRTY_TWENTY,        
                FAO_PA_COMPLIANT                        
        };

Member Function Documentation

double HarvestControlRule::ConditionalConstantCatch	(	const double &	bt,
		const double &	bmsy,
		const double &	msy,
		const double &	fmsy
	)

Implement the Conditional Constant Catch harvest control rule proposed by the IPHC.

Author

Steve Martell

Parameters

bt	available biomass
bmsy	Biomass at MSY
msy	Maximum sustainable yield
fmsy	Fishing mortality rate that achieves MSY.

For this rule, fish at Fmsy if bt>0.8Bmsy and reduce the TAC to msy of the tac > MSY.

If the biomass is less than 0.8Bmsy and greater then 0.4Bmsy, then set fishing mortality rate as a liner function of depletion.

If the biomass is less than 0.4Bmsy, the set the tac=0

Definition at line 221 of file HarvestControlRule.cpp.

{
        double tac;
        double f = 0;
        if( bt >= 0.8*bmsy)
        {
                f = fmsy;
        }
        else if( bt < 0.8*bmsy && bt >= 0.4*bmsy )
        {
                f = (bt-0.4*bmsy)/(0.8-0.4)*fmsy;
        }
        tac = f * bt;
        if(tac > msy)
        {
                tac = msy;
        }
        return tac;
}

Here is the caller graph for this function:

double HarvestControlRule::FortyTen	(	const double &	bt,
		const double &	bo,
		const double &	fmsy
	)

Implement the 40:10 harvest control rule.

Author

Steve Martell

Parameters

<bt>	biomass available at time t
<bo>	theoretical unfished biomass
<fmsy>	Fishing mortality rate that generates Maximum Sustainable Yield.

This function implements the 40:10 harvest control rule that is commonly used by the Pacific Fisheries Management Council.

Definition at line 84 of file HarvestControlRule.cpp.

{
        double dt = bt/bo;
        double ft = fmsy;
        if( dt <= 0.1 )
        {
                ft = 0; 
        } 
        else if( dt > 0.1 && dt <= 0.4 )
        {
                ft = fmsy * (dt-0.1)/(0.4-0.1);
        }
        else if( dt > 0.4 )
        {
                ft = fmsy;
        }

        return bt*(1.-exp(-ft));
}

Here is the caller graph for this function:

double HarvestControlRule::getTac	(	const double &	bt,
		const double &	fmsy,
		const double &	msy,
		const double &	bmsy,
		const double &	bo,
		const double &	delta,
		const double &	ptac,
		const double &	mintac
	)

Return the TAC given estimates of current biomass and reference points.

Author

Steve Martell

Remarks

Uses a swithc statement based on the type of harvest control rule to call the appropriate function to calculate the TAC.

Parameters

bt	preseason biomass forecast
fmsy	estimate of Fmsy
msy	estimate of MSY
bmsy	estimate of Bmsy
bo	estimate of unfished spawning biomass

Definition at line 10 of file HarvestControlRule.cpp.

{
        double tac = 0;
        switch (m_enum)
        {
                case FORTY_TEN:
                        tac = FortyTen(bt, bo, fmsy);
                        cout<<"Forty Ten tac = "<<tac<<endl;
                        break;
                case FIXED_ESCAPEMENT:
                        tac = FixedEscapement(bt,bmsy);
                        cout<<"Fixed Escapement tac = "<<tac<<endl;
                        break;
                case FIXED_ESCAPEMENT_CAP:
                        tac = FixedEscapementCap(bt,bmsy,msy);
                        cout<<"Fixed Escapement cap tac = "<<tac<<endl;
                        break;
                case FIXED_HARVEST_RATE:
                        tac = FixedHarvestRate(bt,fmsy);
                        cout<<"Fixed harvest rate tac = "<<tac<<endl;
                        break;
                case CONDITIONAL_CONSTANT_CATCH:
                        tac = ConditionalConstantCatch(bt,bmsy,msy,fmsy);
                        cout<<"Conditional constant catch = "<<tac<<endl;
                        break;
                case THIRTY_TWENTY:
                        tac = ThirtyTwenty(bt,bo,fmsy);
                        cout<<"Thirty Twenty tac = "<<tac<<endl;
                        break;
                case FIXED_HR_DELTA:
                        tac = FixedHarvestRateDelta(bt,fmsy,delta,ptac);
                        cout<<"Fixed Harvest Rate Delta tac = "<<tac<<endl;
                        break;
                case FAO_PA_COMPLIANT:
                        cout<<"FAO_PA_COMPLIANT"<<endl;
                        break;
                case FLOOR_THIRTY_TWENTY:
                        cout<<"FLOOR_THIRTY_TWENTY"<<endl;
                        tac = FloorThirtyTwenty(bt,bo,fmsy,mintac);
                        break;

        }

        // Put in SUFD adjustment here
        // 50% decrease, 33% increase

        // Evaluation framework, need to be able
        // to tune the MP to achieve the best 
        // tradeoff in the performance measures.


        // Floor 
        return tac;
}

Here is the call graph for this function:

Here is the caller graph for this function:

double HarvestControlRule::ThirtyTwenty	(	const double &	bt,
		const double &	bo,
		const double &	fmsy
	)

30:20 harvest control rule

Author

Steve Martell

Remarks

Parameters

bt	preseson available biomass forecast.

Definition at line 110 of file HarvestControlRule.cpp.

{
        double dt = bt/bo;
        double ft = fmsy;
        if( dt <= 0.2 )
        {
                ft = 0; 
        } 
        else if( dt > 0.2 && dt <= 0.3 )
        {
                ft = fmsy * (dt-0.2)/(0.3-0.2);
        }
        else if( dt > 0.3 )
        {
                ft = fmsy;
        }

        return bt*(1.-exp(-ft));
}

Here is the caller graph for this function:

---

The documentation for this class was generated from the following files:

• /Users/stevenmartell1/Documents/IPHC/MSEdemo/src/HarvestControlRule.h
• /Users/stevenmartell1/Documents/IPHC/MSEdemo/src/HarvestControlRule.cpp