Numerical atmos test runs fixed to not go beyond atmosphere table limits. (#605)

* Numerical atmos test runs fixed to not go beyond atmosphere table limits.

Also made an output name change to allow simultaneous running of the 2 tests.

Author: thjsal

changelog.d/605.attribution.rst

@@ -0,0 +1 @@
+Tuomo Salmi

changelog.d/605.fixed.rst

@@ -0,0 +1 @@
+Fixed the example runs with numerical atmosphere data to not extrapolate beyond the atmosphere table limits.

examples/examples_modeling_tutorial/TestRun_Num.py

@@ -114,7 +114,10 @@ def from_response_files(cls, ARF, RMF, channel_edges, max_input,
 bounds = dict(super_colatitude = (None, None),
               super_radius = (None, None),
               phase_shift = (0.0, 0.1),
-              super_temperature = (5.1, 6.8))
+              super_temperature = (5.3, 6.8))
+#Note: The atmosphere table lower limit is 5.1, but we use 5.3 in this example,
+#because the secondary spot temperature is defined here to be always 0.2 below
+#the primary spot temperature (just for demonstrating how to derive parameters).
 
 primary = xpsi.HotRegion(bounds=bounds,
                         values={},

examples/examples_modeling_tutorial/TestRun_NumBeam.py

@@ -117,10 +117,13 @@ def from_response_files(cls, ARF, RMF, channel_edges, max_input,
 bounds = dict(super_colatitude = (None, None),
               super_radius = (None, None),
               phase_shift = (0.0, 0.1),
-              super_temperature = (5.1, 6.8),
+              super_temperature = (5.3, 6.8),
               super_abb = (-0.7, 0.7),
-              super_bbb = (-0.7, 0.7))              
-	                    
+              super_bbb = (-0.7, 0.7))
+#Note: The atmosphere table lower limit is 5.1, but we use 5.3 in this example,
+#because the secondary spot temperature is defined here to be always 0.2 below
+#the primary spot temperature (just for demonstrating how to derive parameters).
+
 from modules.CustomHotRegion_Beaming import CustomHotRegion
 
 primary =  CustomHotRegion(bounds=bounds,
@@ -539,7 +542,7 @@ def transform(self, p, **kwargs):
                   'importance_nested_sampling': False,
                   'multimodal': False,
                   'n_clustering_params': None,
-                  'outputfiles_basename': './run/run_Num',
+                  'outputfiles_basename': './run/run_NumBeam',
                   'n_iter_before_update': 50,
                   'n_live_points': 50,
                   'sampling_efficiency': 0.8,