initialize share_delta_neff back to False when used from python,e.g. …

…in case starting from copy of .ini file
cmbant · Nov 20, 2024 · c525b16 · c525b16
1 parent ecd375b
commit c525b16
Show file tree

Hide file tree

Showing 4 changed files with 4 additions and 3 deletions.
diff --git a/camb/__init__.py b/camb/__init__.py
@@ -8,7 +8,7 @@
 __author__ = "Antony Lewis"
 __contact__ = "antony at cosmologist dot info"
 __url__ = "https://camb.readthedocs.io"
-__version__ = "1.5.8"
+__version__ = "1.5.9"
 
 from . import baseconfig
 

diff --git a/camb/dark_energy.py b/camb/dark_energy.py
@@ -102,7 +102,7 @@ def validate_params(self):
     def set_w_a_table(self, a, w):
         # check w array has elements that do not cross -1
         if np.sign(1 + np.max(w)) - np.sign(1 + np.min(w)) == 2:
-            raise ValueError('fluid dark energy model does not support w crossing -1')
+            raise CAMBError('fluid dark energy model does not support w crossing -1')
         super().set_w_a_table(a, w)
 
 

diff --git a/fortran/config.f90 b/fortran/config.f90
@@ -3,7 +3,7 @@ module config
     use constants, only: const_twopi
     implicit none
 
-    character(LEN=*), parameter :: version = '1.5.8'
+    character(LEN=*), parameter :: version = '1.5.9'
 
     integer :: FeedbackLevel = 0 !if >0 print out useful information about the model
 

diff --git a/fortran/model.f90 b/fortran/model.f90
@@ -252,6 +252,7 @@ subroutine CAMBparams_SetNeutrinoHierarchy(this, omnuh2, omnuh2_sterile, nnu, ne
         return
     end if
     this%Nu_mass_eigenstates=0
+    this%share_delta_neff = .false.
     if (omnuh2 > omnuh2_sterile) then
         normal_frac =  (omnuh2-omnuh2_sterile)/omnuh2
         if (neutrino_hierarchy == neutrino_hierarchy_degenerate) then