Implementing new recruitment dev vector approach

Imported from redmine, Issue #57889
Opened by @RickMethot on 2018-11-23
Status when imported: In Progress

create new recruitment deviation vector approach
do_recdev==3

With this approach:
exp_rec(y,4) = R0_adj * mfexp(dev(y))
where R0_adj is R0 adjusted for blocks, etc. on R0;
and dev(y) is a simple vector of unconstrained deviations

Then in the objective function:
rather than getting the -logL from the dev vector, now create a new deviation which is dev=log(exp_rec(y,3)/exp_rec(y,4));
where exp_rec(y,4) is shown above and exp_rec(y,3) is the expected value of recruitment adjusted for spawn_recr curve, regimes and env, and bias-adjustment.

Initial test shows that the -logL is about same as using do_recrdev==2 which is simple devs in which exp_rec(y,4) = exp_rec(y,3) * mfexp(dev(y))