From 5ac3391500635ab2c05a4c821a443791fbea83f4 Mon Sep 17 00:00:00 2001 From: hposborn Date: Wed, 10 Jul 2024 13:57:49 +0200 Subject: [PATCH] fixing phase generation issue for trios --- MonoTools/fit.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/MonoTools/fit.py b/MonoTools/fit.py index 44e0943..d9e0119 100755 --- a/MonoTools/fit.py +++ b/MonoTools/fit.py @@ -1089,10 +1089,10 @@ def init_lc(self, oot_binsize=1/48, **kwargs): p=np.max(self.init_soln['per_'+pl]) if hasattr(self,'init_soln') else np.max(self.planets[pl]['period_aliases']) if self.model_t03_ttv: t0s= [self.init_soln['t0_'+pl],self.init_soln['t0_2_'+pl],self.init_soln['t0_3_'+pl]] if hasattr(self,'init_soln') else [self.planets[pl]['tcen'],self.planets[pl]['tcen_2'],self.planets[pl]['tcen_3']] - + phase=self.make_phase(self.lc.time,t0s,p) else: t0= self.init_soln['t0_'+pl] if hasattr(self,'init_soln') else self.planets[pl]['tcen_2'] - phase=self.make_phase(self.lc.time,t0s,p) + phase=self.make_phase(self.lc.time,t0,p) elif pl in self.duos: t0= self.init_soln['t0_'+pl] if hasattr(self,'init_soln') else self.planets[pl]['tcen'] p=abs(self.init_soln['t0_2_'+pl]-self.init_soln['t0_'+pl]) if hasattr(self,'init_soln') else abs(self.planets[pl]['tcen_2']-self.planets[pl]['tcen'])