From 16cd87824c6eff20a1ea47da30ff93b9be700fdb Mon Sep 17 00:00:00 2001
From: Ian Taylor <4992918+iantaylor-NOAA@users.noreply.github.com>
Date: Fri, 24 May 2024 10:43:20 -0400
Subject: [PATCH] add fatal error when using MCMC with recdev option 1 (#590)

- inspired by discussion: https://github.com/nmfs-ost/ss3-source-code/discussions/566
---
 SS_readcontrol_330.tpl | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/SS_readcontrol_330.tpl b/SS_readcontrol_330.tpl
index 0e4dfab2..7c9568ef 100644
--- a/SS_readcontrol_330.tpl
+++ b/SS_readcontrol_330.tpl
@@ -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;