Skip to content

first commit for protection level #729

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
Rick-Methot-NOAA merged 2 commits into from
Sep 18, 2025
Merged

first commit for protection level #729

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a80a78d
first commit for protection level
Rick-Methot-NOAA Sep 12, 2025
7b92162
couple of AEP fixes to code
Rick-Methot-NOAA Sep 15, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 7 additions & 1 deletion SS_benchfore.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2836,7 +2836,13 @@ FUNCTION void Get_Forecast()
(1.0) * (1.0 - join2); // scale right side
break;
}
}
}
// protection level; contributed by AEP
if (H4010_prot > 0)
{
join1 = 1.0 / (1. + mfexp(20.0 * (SSB_current - H4010_prot * HCR_anchor)));
ABC_buffer(y) = 0.0001*join1 + (1.0 - 0.0001)*(1.0 - join1)*ABC_buffer(y);
}
} // end calc of ABC buffer
else
{ // ABC buffer remains at previously calculated value
Expand Down
15 changes: 13 additions & 2 deletions SS_readdata_330.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4029,7 +4029,9 @@
int Rebuild_Yinit;
int HarvestPolicy; // 0=none; 1=west coast adjust catch; 2=AK to adjust F
number H4010_top_rd;
number H4010_bot_rd;
number H4010_bot;
number H4010_prot;
number H4010_scale;
number H4010_scale_rd;
int Do_Impl_Error;
Expand Down Expand Up @@ -4283,8 +4285,17 @@
echoinput << HarvestPolicy << " # echoed HarvestPolicy " << endl;
*(ad_comm::global_datafile) >> H4010_top_rd; // as fraction of HCR_anchor; use -1 as legacy approach to set H4010_top to Bmsy/SSB_unf
echoinput << H4010_top_rd << " # echoed control rule inflection" << endl;
*(ad_comm::global_datafile) >> H4010_bot;
echoinput << H4010_bot << " # echoed control rule cutoff " << endl;
*(ad_comm::global_datafile) >> H4010_bot_rd;
echoinput << H4010_bot_rd << " # echoed control rule cutoff " << endl;
H4010_bot = H4010_bot_rd;
H4010_prot = -1000;
if (H4010_bot < 0)
{
H4010_bot *= -1;
echoinput << H4010_bot << " # echoed control rule cutoff " << endl;
*(ad_comm::global_datafile) >> H4010_prot;
echoinput << H4010_prot << " # echoed protection level cutoff " << endl;
}
*(ad_comm::global_datafile) >> H4010_scale_rd;
H4010_scale = H4010_scale_rd;
echoinput << H4010_scale << " # echoed control rule scalar " << endl;
Expand Down
6 changes: 5 additions & 1 deletion SS_write_ssnew.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1645,7 +1645,11 @@ FUNCTION void write_nucontrol()
NuFore << HarvestPolicy << " # Control rule method (0: none; 1: ramp does catch=f(SSB), buffer on F; 2: ramp does F=f(SSB), buffer on F; 3: ramp does catch=f(SSB), buffer on catch; 4: ramp does F=f(SSB), buffer on catch) " << endl;
NuFore << "# values for top, bottom and buffer required, but not used when Policy=0" << endl;
NuFore << H4010_top_rd << " # Control rule inflection for constant F (as frac of HCR_anchor, see below); must be > control rule cutoff" << endl;
NuFore << H4010_bot << " # Control rule cutoff for no F (as frac of HCR_anchor, e.g. 0.10) " << endl;
NuFore << H4010_bot_rd << " # Control rule cutoff for no F (as frac of HCR_anchor, e.g. 0.10); negative value to also invoke read of protection level " << endl;
if (H4010_bot_rd < 0)
{ NuFore << H4010_prot << " # AK protection level: control rule level below which F goes to 0.0001" << endl;}
else
{ NuFore << " # placeholder for AK protection level: control rule level below which F goes to 0.0001" << endl;}
NuFore << H4010_scale_rd << " # Buffer: enter Control rule target as fraction of Flimit (e.g. 0.75), negative value invokes list of [year, scalar]. -year fills from year to YrMax " << endl;
NuFore << "# Also see HCR_anchor below to use virgin vs benchmark SSB or Bmsy as basis for inflection and cutoff" << endl;
if (H4010_scale_rd < 0)
Expand Down
Loading