ADD projection component

[Craig44]

Recruitment

Options for projected recruit devs

• [x] Simulate from normal distribution with sigma R data$future_recruitment_type = 0
• [x] Empirically resample from input recruitment deviations over a given period data$future_recruitment_type = 1 use data$year_ndx_for_empirical_resampling to specify period to resample

Fishing mortality

• [x] User specifies F values during projection (unlikely useage) data$future_fishing_type = 0 populate data$future_fishing_inputs_trwl and data$future_fishing_inputs_fixed with F values
• [x] User specifies Catch values during projection (more likely) data$future_fishing_type = 1 populate data$future_fishing_inputs_trwl and data$future_fishing_inputs_fixed with future catches
• [x] The model will calculate a Fmsy based on long-term deterministic productivity and apply that during the projection period. Create an R function find_FMSY which will search over a range of fishing mortality rates to find the one that yields a desired projected biomass outcome.

[Craig44]

Created a function find_fref which will find an F that achieves some user defined global (summed across all regions) percent B0 target. See commit e4ef5e10ddd1ae4069a3b78849cfb782738a62a6.

We will probably be modifying these to allow for more complex conditions to solve for F's in different regions. Although this function will satisfy some overall SSB target, some regions may be at unacceptable levels.