Skip to content

augment F reporting #718

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 4 commits into from
Sep 9, 2025
Merged

augment F reporting #718

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
bfcb720
augment F reporting
Sep 4, 2025
5e3cd4d
better with year dim on the array
Sep 4, 2025
2723506
change name to Fdead
Sep 5, 2025
d1390d3
correct Fdead terminology
Rick-Methot-NOAA Sep 9, 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
1 change: 1 addition & 0 deletions SS_param.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -533,6 +533,7 @@ PARAMETER_SECTION
4darray sel_dead_bio(1,nseas,1,Nfleet,1,gmorph,0,nages); // sel * (retain + (1-retain)*discmort) * wt

4darray save_sel_num(styr-3*nseas,TimeMax_Fcast_std+nseas,1,Nfleet,1,gmorph,0,nages) // save sel_num (Asel_2) and save fecundity for output; +nseas covers no forecast setups
4darray save_F_dead(styr-3*nseas,TimeMax_Fcast_std+nseas,1,Nfleet,1,gmorph,0,nages) // save sel_num (Asel_2) and save fecundity for output; +nseas covers no forecast setups

4darray Sel_for_tag(TG_timestart*Do_TG,TimeMax*Do_TG,1,Nfleet,1,gmorph*Do_TG,0,nages)
vector TG_report(1,Nfleet*Do_TG);
Expand Down
4 changes: 2 additions & 2 deletions SS_selex.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2053,7 +2053,7 @@ FUNCTION void Make_FishSelex()
<< " " << f << " " << Wt_Age_t(tz, f, g) << " #wt_flt_" << f << endl;
}
} // end need to do it
save_sel_num(t, f, g) = sel_num(s, f, g); // save sel_num in save_fecundity array for output

save_sel_num(t, f, g) = sel_num(s, f, g); // save sel_num for output in age_selex report
save_F_dead(t, f, g) = sel_dead_num(s, f, g); // save for output
} // end fleet loop for mortality, retention
} // end Make_FishSelex
2 changes: 2 additions & 0 deletions SS_write_report.tpl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3096,6 +3096,7 @@ FUNCTION void write_bigoutput()
SS2out << "Asel_is_age_selectivity_alone" << endl;
SS2out << "Asel2_is_Asel*(selL*size_at_age(ALK)); Q and F parameters may appear higher than expected because Asel2 may have max < 1.0; " << endl;
SS2out << "Aret_is_age_retention" << endl;
SS2out << "F_is_F*Asel2;_Fdead_is_F*sel_dead_so_matches_the_mortality_calculations" << endl;
SS2out << "COMBINED_ALK*selL*selA*wtlen*ret*discmort_in_makefishsel_yr: " << makefishsel_yr << " With_MeanSel_From: " << Fcast_Sel_yr1 << " - " << Fcast_Sel_yr2; // SS_Label_380
SS2out << "; Year_styr-3_(" << styr - 3 << ")_stores_average_used_for_benchmark" << endl;

Expand Down Expand Up @@ -3164,6 +3165,7 @@ FUNCTION void write_bigoutput()

if (fleet_type(f) <= 2)
SS2out << "F " << f << " " << y << " " << s << " " << sx(g) << " " << g << " " << y << "_" << f << "_F" << Hrate(f, t) * save_sel_num(t, f, g) << endl;
SS2out << "Fdead " << f << " " << y << " " << s << " " << sx(g) << " " << g << " " << y << "_" << f << "_Fdead" << Hrate(f, t) * save_F_dead(t, f, g) << endl;
SS2out << "bodywt " << f << " " << y << " " << s << " " << sx(g) << " " << g << " " << y << "_" << f << "_bodywt" << Wt_Age_t(t, f, g) << endl;
}
}
Expand Down
Loading