fixing tlusty Z determination

karllark · Nov 22, 2023 · 455722a · 455722a
1 parent 765e5db
commit 455722a
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/measure_extinction/utils/fit_model.py b/measure_extinction/utils/fit_model.py
@@ -162,7 +162,8 @@ def lnprob(params, obsdata, modeldata, fitinfo):
             information about the fitting
         """
         lnp = fitinfo.lnprior(params)
-        print(lnp, fitinfo.lnlike(params, obsdata, modeldata))
+        # print(params)
+        # print(lnp, fitinfo.lnlike(params, obsdata, modeldata))
         if lnp == lnp_bignnum:
             return lnp
         else:

diff --git a/measure_extinction/utils/make_all_tlusty_obsdata.py b/measure_extinction/utils/make_all_tlusty_obsdata.py
@@ -50,7 +50,10 @@ def decode_params(filename):
     gpos = filename.find("g", slashpos)
     vpos = filename.find("v", slashpos)
 
-    model_params["Z"] = float(filename[zpos + 1 : tpos]) * 0.01
+    if tpos - zpos > 4:
+        model_params["Z"] = float(filename[zpos + 1 : tpos]) * 0.001
+    else:
+        model_params["Z"] = float(filename[zpos + 1 : tpos]) * 0.01
     model_params["Teff"] = float(filename[tpos + 1 : gpos])
     model_params["logg"] = float(filename[gpos + 1 : vpos]) * 0.01
     model_params["vturb"] = float(filename[vpos + 1 : periodpos])