From c9112f1dcaffa710a2cd3f857d32a2db1cfec756 Mon Sep 17 00:00:00 2001
From: Christopher Wipf <wipf@ligo.mit.edu>
Date: Thu, 9 Apr 2020 14:24:17 -0700
Subject: [PATCH] coatingthermal.py: formatting fixes
---
gwinc/noise/coatingthermal.py | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/gwinc/noise/coatingthermal.py b/gwinc/noise/coatingthermal.py
index 19c411c0..62bb4356 100644
--- a/gwinc/noise/coatingthermal.py
+++ b/gwinc/noise/coatingthermal.py
@@ -89,9 +89,10 @@ def coating_brownian(f, materials, wavelength, wBeam, dOpt):
zdir = -1
dcdp_z = zdir * dcdp # z-dir only matters here
- # layer contributions (b_j in PhysRevD.91.042002, errors corrected)
- brLayer = ((1-nN*dcdp_z/2)**2*((1-2*pratN)*(1+pratN)*Ysub)/((1-2*pratsub)*(1+pratsub)*yN) +
- ((1-2*pratsub)*(1+pratsub)*yN)/((1+pratN)*Ysub))/(1-pratN)
+ # layer contributions, b_j (eq 1) from doi:10.1103/PhysRevD.91.042002, errors corrected
+ brLayer = 1/(1-pratN) *
+ ( (1-nN*dcdp_z/2)**2 * (1-2*pratN)*(1+pratN)*Ysub / ((1-2*pratsub)*(1+pratsub)*yN) +
+ (1-2*pratsub)*(1+pratsub)*yN / ((1+pratN)*Ysub) )
# sum them up for total
w = 2 * pi * f
--
GitLab