From c8e2544c71007851da3aaad3c0ca13d395da81da Mon Sep 17 00:00:00 2001 From: Noah Franz Date: Fri, 2 Aug 2024 09:12:09 -0500 Subject: [PATCH] modify input kwargs to allow for arguments to be passed to Lnu_of_nu in the MQ24 code --- src/syncfit/_version.py | 2 +- src/syncfit/models/mq_model.py | 8 ++------ 2 files changed, 3 insertions(+), 7 deletions(-) diff --git a/src/syncfit/_version.py b/src/syncfit/_version.py index e19434e..334b899 100644 --- a/src/syncfit/_version.py +++ b/src/syncfit/_version.py @@ -1 +1 @@ -__version__ = "0.3.3" +__version__ = "0.3.4" diff --git a/src/syncfit/models/mq_model.py b/src/syncfit/models/mq_model.py index 1fd96a4..9d7ca02 100644 --- a/src/syncfit/models/mq_model.py +++ b/src/syncfit/models/mq_model.py @@ -37,14 +37,10 @@ def __init__(self, prior=None, p=None, t=None): self.t = t def SED(self, nu, p, log_bG_sh, log_Mdot, log_epsilon_T, log_epsilon_e, log_epsilon_B, - t, lum_dist, **kwargs): + t, lum_dist, ell_dec = 1.0, f = 3.0/16.0, **kwargs): # set microphysical and geometric parameters - # log_epsilon_e = -1 - # log_epsilon_B = log_epsilon_e # assume equipartition delta = 10**log_epsilon_e/10**log_epsilon_T - f = 3.0/16.0 - ell_dec = 1.0 t = (t*u.day).to(u.s).value @@ -53,7 +49,7 @@ def SED(self, nu, p, log_bG_sh, log_Mdot, log_epsilon_T, log_epsilon_e, log_epsi Lnu = Lnu_of_nu( 10**log_bG_sh, Mdot_over_vw, nu, t, p=p, epsilon_T=10**log_epsilon_T, epsilon_B=10**log_epsilon_B, epsilon_e=10**log_epsilon_e, - f=f,ell_dec=ell_dec,radius_insteadof_time=False + f=f,ell_dec=ell_dec, **kwargs ) * u.erg / (u.s * u.Hz) lum_dist_cm = lum_dist*u.cm # give it units so the conversion works well