diff --git a/gwinc/noise/suspensionthermal.py b/gwinc/noise/suspensionthermal.py
index 91b00ef11febb9de85edb20248cc9a2d25af6c65..2dd31bfc844806400444d41010194d165bf4e9c0 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 5c3f9059baf84ae5255cf0d6d86730451a9d230a..8ea036007899e35fdd8309ff10c628512bfce5fd 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'