Merge pull request #290 from nmfs-stock-synthesis/fix_inityr_catage_r…

…eport fix reporting of catage in equilibrium and initial year report
nmfs-ost · Apr 13, 2022 · d120ea5 · d120ea5
2 parents 1486a3d + b0a7b2e
commit d120ea5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/SS_popdyn.tpl b/SS_popdyn.tpl
@@ -1819,10 +1819,10 @@ FUNCTION void Do_Equil_Calc(const prevariable& equ_Recr)
                  equ_catch_fleet(2,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_dead_bio(s,f,g))*Zrate2(p,g);      // dead catch bio
                  equ_catch_fleet(5,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_dead_num(s,f,g))*Zrate2(p,g);      // deadfish catch numbers
                  equ_catch_fleet(3,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_ret_bio(s,f,g))*Zrate2(p,g);      // retained catch bio
-                 equ_catage(s,f,g)=elem_prod(elem_prod(equ_numbers(s,p,g)(0,nages),sel_dead_num(s,f,g)) , Zrate2(p,g));
                  equ_catch_fleet(1,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_bio(s,f,g))*Zrate2(p,g);      // encountered catch bio
                  equ_catch_fleet(4,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_num(s,f,g))*Zrate2(p,g);      // encountered catch bio
                  equ_catch_fleet(6,s,f)+=Hrate(f,t)*elem_prod(equ_numbers(s,p,g)(0,nages),sel_ret_num(s,f,g))*Zrate2(p,g);      // retained catch numbers
+                 equ_catage(s,f,g)=Hrate(f,t)*elem_prod(elem_prod(equ_numbers(s,p,g)(0,nages),sel_dead_num(s,f,g)) , Zrate2(p,g));
                }
              }
              else  // F_method=1