From 070abb186c8fc39d0171dcaf5697a686cacc7b03 Mon Sep 17 00:00:00 2001
From: Christopher Wipf <wipf@ligo.mit.edu>
Date: Fri, 10 Aug 2018 14:21:06 -0700
Subject: [PATCH] Fixes for cryo silicon suspensions Suspension type BQuad
changes the bottom stage material to silicon. Reverse the SUS stage list in
suspensionthermal.py to use the appropriate Temp. Addresses #23
---
gwinc/noise/suspensionthermal.py | 2 +-
gwinc/suspension.py | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/gwinc/noise/suspensionthermal.py b/gwinc/noise/suspensionthermal.py
index 91b00ef1..2dd31bfc 100644
--- a/gwinc/noise/suspensionthermal.py
+++ b/gwinc/noise/suspensionthermal.py
@@ -61,7 +61,7 @@ def susptherm(f, ifo):
##########################################################
dxdF = zeros(hForce.shape, dtype=complex)
- for n, stage in enumerate(ifo.Suspension.Stage):
+ for n, stage in enumerate(reversed(ifo.Suspension.Stage)):
# add up the contribution from each stage
# convert to beam line motion. theta is squared because
diff --git a/gwinc/suspension.py b/gwinc/suspension.py
index 5c3f9059..8ea03600 100644
--- a/gwinc/suspension.py
+++ b/gwinc/suspension.py
@@ -164,7 +164,7 @@ def suspQuad(f, ifo, material='Silica'):
if 'WireMaterial' in stage:
WireMaterial = stage.WireMaterial
elif n == last_stage:
- WireMaterial = 'Silica'
+ WireMaterial = material
else:
WireMaterial = 'C70Steel'
--
GitLab