aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2015-08-11 11:08:22 +0200
committerVolker Hoffmann <volker@cheleb.net>2015-08-11 11:08:22 +0200
commit913df3c914e54ed6c62c3dc93bee9417702b379e (patch)
tree6bc6d1b022ef66b8d06ef163339debb758cacc42
parent65b10bef6144dfcb1398f08427a8f6cf3c0443d3 (diff)
fix: calibrate only to data after blowup (if no blowup happened in the model [solar2/solar2_hd], we still don't have data before this)
-rw-r--r--Helpers/impact_helpers.py13
1 files changed, 12 insertions, 1 deletions
diff --git a/Helpers/impact_helpers.py b/Helpers/impact_helpers.py
index 86fea15..dec516d 100644
--- a/Helpers/impact_helpers.py
+++ b/Helpers/impact_helpers.py
@@ -356,8 +356,19 @@ def load_all():
"Morby", "Morby/HD", \
"Blowup", "Blowup2", "Blowup2/HD" ]
+ #
# Time Range for Calibration
- tlo_all = np.array([0.0, 0.0, \
+ #
+ # NB:
+ # For Solar2/Solar2_HD runs, the lower time limit for calibration is the
+ # time after the blow up of the model. ALthough there is no blow-up in
+ # this model, we only have actual evidence for the time after this. So
+ # there's no point in calibrating to the time before.
+ #
+ # For Morby/Morby_HD runs, the lower time limit is a dummy values. All
+ # calibrations for these are obtained from Blowup2/Blowup2_HD runs
+ #
+ tlo_all = np.array([tblowup, tblowup, \
0.0, 0.0, \
tblowup, tblowup, tblowup ])
thi_all = np.ones_like(tlo_all) * np.nan