Update proposal distribution for kinetics

OpenThermochronology · Apr 14, 2024 · 5d82fdd · 5d82fdd
1 parent a28b662
commit 5d82fdd
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/src/inversion.jl b/src/inversion.jl
@@ -882,10 +882,10 @@
 
         # Standard deviations of Gaussian proposal ("jumping") distributions
         # for kinetics
-        σaD0 = adm.D0L_logsigma
-        σaEa = adm.EaL_logsigma
-        σzD0 = (zdm.DzD0_logsigma+zdm.DN17D0_logsigma)
-        σzEa = (zdm.DzEa_logsigma+zdm.DN17Ea_logsigma)
+        σaD0 = adm.D0L_logsigma/10
+        σaEa = adm.EaL_logsigma/10
+        σzD0 = (zdm.DzD0_logsigma+zdm.DN17D0_logsigma)/20
+        σzEa = (zdm.DzEa_logsigma+zdm.DN17Ea_logsigma)/20
         # for temperature and time
         σⱼt = tinit/60
         σⱼT = Tinit/60
@@ -1278,10 +1278,10 @@
 
         # Standard deviations of Gaussian proposal ("jumping") distributions
         # for kinetics
-        σaD0 = adm.D0L_logsigma
-        σaEa = adm.EaL_logsigma
-        σzD0 = (zdm.DzD0_logsigma+zdm.DN17D0_logsigma)
-        σzEa = (zdm.DzEa_logsigma+zdm.DN17Ea_logsigma)
+        σaD0 = adm.D0L_logsigma/10
+        σaEa = adm.EaL_logsigma/10
+        σzD0 = (zdm.DzD0_logsigma+zdm.DN17D0_logsigma)/20
+        σzEa = (zdm.DzEa_logsigma+zdm.DN17Ea_logsigma)/20
         # for temperature and time
         σⱼt = tinit/60
         σⱼT = Tinit/60