Skip to content

Commit

Permalink
add fatal error when using MCMC with recdev option 1 (#590)
Browse files Browse the repository at this point in the history
- inspired by discussion: #566
  • Loading branch information
iantaylor-NOAA authored and e-perl-NOAA committed Nov 4, 2024
1 parent 95be4a1 commit 16cd878
Showing 1 changed file with 11 additions and 0 deletions.
11 changes: 11 additions & 0 deletions SS_readcontrol_330.tpl
Original file line number Diff line number Diff line change
Expand Up @@ -2230,6 +2230,17 @@

init_int do_recdev // 0=none; 1=devvector; 2=simple deviations; 3=dev from R0
!!echoinput<<do_recdev<<" do_recdev"<<endl;
// check for use of devvector with MCMC
LOCAL_CALCS
// clang-format on
if (do_recdev == 1 & mcmcFlag == 1)
{
warnstream << "do_recdev option 1=devvector should not be used with MCMC, recommend option 2=simple deviations. For more detail see https://github.com/admb-project/admb/issues/107.";
write_message (FATAL, 0);
}
// clang-format off
END_CALCS

init_int recdev_start;
!!echoinput<<recdev_start<<" recdev_start"<<endl;
init_int recdev_end;
Expand Down

0 comments on commit 16cd878

Please sign in to comment.