From 3faf74e3335c591948dd18389d679a3fbd85083e Mon Sep 17 00:00:00 2001
From: Robert Izzard <r.izzard@surrey.ac.uk>
Date: Fri, 27 Sep 2019 07:44:14 +0100
Subject: [PATCH] change default lambda_ce to 0.5 because setting it to -1 is
 not valid over the whole parameter space, especially not at Z<0.02

---
 src/binary_c_parameters.h | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/src/binary_c_parameters.h b/src/binary_c_parameters.h
index d1d44db8e..781e34b7e 100644
--- a/src/binary_c_parameters.h
+++ b/src/binary_c_parameters.h
@@ -1093,9 +1093,14 @@
 
 /*
  * Lambda is the envelope binding energy parameter:
- * Ebind = -G M Menv / (R lambda) 
+ * Ebind = -G M Menv / (R lambda)
+ *
+ * Note: much as we'd like to set this to one of the 
+ *       three LAMBDA_* macros (for the Dewi/Tauris, 
+ *       Wang or Polytrope) these are not reliable over
+ *       the whole parameter space. Instead, use 0.5.
  */
-#define DEFAULT_LAMBDA_CE -1.0
+#define DEFAULT_LAMBDA_CE 0.5
 
 /* use Hurley et al. (2002) energy balance prescription by default */
 #define DEFAULT_COMENV_PRESCRIPTION COMENV_BSE
-- 
GitLab