The Bowman infrasound pressure model

gwinc/noise/newtonian.py

@@ -220,3 +220,17 @@ def atmois(f, ifo):
     n_if = coupling_if * psd_if
 
     return n_if * ifo.gwinc.sinc_sqr
+
+
+def atmoBowman(f):
+    """The Bowman infrasound model
+   
+    """ 
+    freq = np.array([
+                0.01, 0.0155, 0.0239, 0.0367, 0.0567,
+                0.0874, 0.1345, 0.2075, 0.32, 0.5,
+                0.76, 1.17, 1.8, 2.79, 4.3,
+                6.64, 10, 100])
+    pressure_asd = np.sqrt([22.8, 4, 0.7, 0.14, 0.027, 0.004, 0.0029, 0.0039, 7e-4,
+                                              1.44e-4, 0.37e-4, 0.12e-4, 0.56e-5, 0.35e-5, 0.26e-5, 0.24e-5, 2e-6, 2e-6])
+    return 10**(np.interp(np.log10(f), np.log10(freq), np.log10(pressure_asd)))