Actually include phase effect in Hy

lsst-sssc · Sep 27, 2024 · 4437f63 · 4437f63
1 parent 4a1dfac
commit 4437f63
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/models/activity.py b/models/activity.py
@@ -81,4 +81,9 @@ def Hy(H, y, rh, delta, phase):
 
     """
 
-    return H + 5 * np.log10(rh * delta) - (2.5 * y) * np.log10(rh)
+    return (
+        H
+        + 5 * np.log10(rh * delta)
+        - (2.5 * y) * np.log10(rh)
+        - 2.5 * np.log10(schleicher_marcus(phase))
+    )