couple of AEP fixes to code

Rick-Methot-NOAA · Rick-Methot-NOAA · commit 7656e8b3bf23 · 2025-09-18T11:29:02.000-07:00
diff --git a/SS_benchfore.tpl b/SS_benchfore.tpl
@@ -2836,15 +2836,13 @@ FUNCTION void Get_Forecast()
                       (1.0) * (1.0 - join2); // scale right side
                   break;
                 }
-              //AEP AEP Add here (protection level)
+              }
+              // protection level; contributed by AEP
               if (H4010_prot > 0)
               {
-                join1 = 1.0 / (1. + mfexp(20.0 * (SSB_current - H4010_prot * temp)));
-                //cout << ABC_buffer(y) << endl;
+                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);
-                //cout  << SSB_current << " " << (SSB_current - H4010_prot * temp) << " " << SSB_current/temp << " " << join1 << " " << ABC_buffer(y) << endl;
               }
-            }
           } // end calc of ABC buffer
           else
           { //  ABC buffer remains at previously calculated value

Original file line number	Diff line number	Diff line change
`@@ -2836,15 +2836,13 @@ FUNCTION void Get_Forecast()`
`2836`	`2836`	`(1.0) * (1.0 - join2); // scale right side`
`2837`	`2837`	`break;`
`2838`	`2838`	`}`
`2839`		`- //AEP AEP Add here (protection level)`
	`2839`	`+ }`
	`2840`	`+ // protection level; contributed by AEP`
`2840`	`2841`	`if (H4010_prot > 0)`
`2841`	`2842`	`{`
`2842`		`- join1 = 1.0 / (1. + mfexp(20.0 * (SSB_current - H4010_prot * temp)));`
`2843`		`- //cout << ABC_buffer(y) << endl;`
	`2843`	`+ join1 = 1.0 / (1. + mfexp(20.0 * (SSB_current - H4010_prot * HCR_anchor)));`
`2844`	`2844`	`ABC_buffer(y) = 0.0001join1 + (1.0 - 0.0001)(1.0 - join1)*ABC_buffer(y);`
`2845`		`- //cout << SSB_current << " " << (SSB_current - H4010_prot * temp) << " " << SSB_current/temp << " " << join1 << " " << ABC_buffer(y) << endl;`
`2846`	`2845`	`}`
`2847`		`- }`
`2848`	`2846`	`} // end calc of ABC buffer`
`2849`	`2847`	`else`
`2850`	`2848`	`{ // ABC buffer remains at previously calculated value`