-
Notifications
You must be signed in to change notification settings - Fork 17
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.
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
adjust warnings for 2DAR #617
adjust warnings for 2DAR #617
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me. Thanks for all the explorations.
@e-perl-NOAA Please reflect these revised warnings into the manual. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I am uncertain of the level of guidance that you should be giving the user. It seems like you are just warning them and then changing it behind the scenes, where I prefer that the model just stops if they try to estimate it.
} | ||
if (use_rho == 1) | ||
{ | ||
echoinput << "read two parameter lines for rho_yr and then rho_age (or length)" << endl; | ||
warnstream << "2DAR rho is incompletely implemented; it should only be used experimentally and never estimated"; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I am hesitant to allow users to even do this if it is incompletely implemented. I would prefer a full stop if use_rho = 1
. Especially because in this same if statement you set dtempvec(7) = -1
. Which, does that just change the values put in ss_new or does it actually turn off estimation without the user knowing?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The estimation has never been enabled. The value of rho is only used in the preliminary calcs section to create a cor matrix.
Concisely describe what has been changed/addressed in the pull request.
Find that the rho feature of 2DAR was incompletely implemented and a full fix is not feasible at this time. Provide new warnings and guidance to users.
What tests have been done?
Numerous runs with various permutations of estimating or fixing age-specific sigmasel values using hake 2DAR example.
Where are the relevant files?
What tests/review still need to be done?
none
Is there an input change for users to Stock Synthesis?
Additional information (optional).