#74 Reformat codes using non-beta Black

deshima_sensitivity/atmosphere.py

@@ -44,8 +44,8 @@ def eta_atm_func(
         Atmospheric tranmsmission. Units: None.
 
     """
-    if np.average(F) > 10.0 ** 9:
-        F = F / 10.0 ** 9
+    if np.average(F) > 10.0**9:
+        F = F / 10.0**9
 
     # give F a length if it is an integer.
     if not hasattr(F, "__len__"):

deshima_sensitivity/filter.py

@@ -80,9 +80,9 @@ def eta_filter_lorentzian(
         Units: None.
     """
 
-    if np.average(F) < 10.0 ** 9:
-        F = F * 10.0 ** 9
-        FWHM = FWHM * 10.0 ** 9
+    if np.average(F) < 10.0**9:
+        F = F * 10.0**9
+        FWHM = FWHM * 10.0**9
 
     # give F a length if it is an integer.
     if not hasattr(F, "__len__"):
@@ -165,8 +165,8 @@ def eta_filter_csv(
 
     F_int = eta_filter_df.columns.values.astype(float)
 
-    if np.average(F_int) < 10.0 ** 9:
-        F_int = F_int * 10.0 ** 9
+    if np.average(F_int) < 10.0**9:
+        F_int = F_int * 10.0**9
         eta_filter_df.columns = F_int
 
     # Fit to lorentzian model
@@ -265,9 +265,9 @@ def expand_F(
 
     """
 
-    if np.average(F) < 10.0 ** 9:
-        F = F * 10.0 ** 9
-        FHWM = FHWM * 10 ** 9
+    if np.average(F) < 10.0**9:
+        F = F * 10.0**9
+        FHWM = FHWM * 10**9
 
     N = len(F)
     if N > 1:

deshima_sensitivity/galaxy.py

@@ -114,4 +114,4 @@ def co_luminosity(z, d_l, f_obs, int):
     c2 = 1.04e-03
 
     # Lp  = c1 * int * d_l**2 / (f_obs**2 * (1 + z)**3 )
-    return c2 * int * d_l ** 2 * f_obs
+    return c2 * int * d_l**2 * f_obs

deshima_sensitivity/instruments.py

@@ -12,7 +12,7 @@
 
 
 # constants
-Delta_Al = 188.0 * 10 ** -6 * e  # gap energy of Al
+Delta_Al = 188.0 * 10**-6 * e  # gap energy of Al
 eta_pb = 0.4  # Pair breaking efficiency
 eta_Al_ohmic_850 = 0.9975  # Ohmic loss of an Al surface at 850 GHz.
 # Shitov+, ISSTT2008. https://www.nrao.edu/meetings/isstt/papers/2008/2008263266.pdf

deshima_sensitivity/physics.py

@@ -11,9 +11,9 @@
 
 
 # constants
-h = 6.62607004 * 10 ** -34  # Planck constant
-k = 1.38064852 * 10 ** -23  # Boltzmann constant
-e = 1.60217662 * 10 ** -19  # electron charge
+h = 6.62607004 * 10**-34  # Planck constant
+k = 1.38064852 * 10**-23  # Boltzmann constant
+e = 1.60217662 * 10**-19  # electron charge
 c = 299792458.0  # velocity of light
 
 

deshima_sensitivity/plotting.py

@@ -104,7 +104,7 @@ def MDLF_simple(
     ax.set_ylabel(r"Minimum Detectable Line Flux ($\mathrm{W\, m^{-2}}$)")
     ax.set_yscale("log")
     ax.set_xlim(200, 460)
-    ax.set_ylim([10 ** -20, 10 ** -17])
+    ax.set_ylim([10**-20, 10**-17])
     ax.tick_params(direction="in", which="both")
     ax.grid(True)
     ax.set_title(
@@ -197,7 +197,7 @@ def MS_simple(F: ArrayLike, pwv: float = 0.5, EL: float = 60.0) -> HTML:
     ax.set_ylabel(r"Mapping Speed ($\mathrm{arcmin^2\, mJy^{-2}\, h^{-1}}$)")
     ax.set_yscale("log")
     ax.set_xlim(200, 460)
-    ax.set_ylim([10 ** -5, 10 ** -2])
+    ax.set_ylim([10**-5, 10**-2])
     ax.tick_params(direction="in", which="both")
     ax.grid(True)
     ax.set_title(

deshima_sensitivity/simulator.py

@@ -340,7 +340,7 @@ def spectrometer_sensitivity(
 
     # Collect efficiencies at the same temperature
     eta_M1 = eta_M1_ohmic * eta_M1_spill
-    eta_wo = eta_Al_ohmic ** n_wo_mirrors * eta_wo_spill
+    eta_wo = eta_Al_ohmic**n_wo_mirrors * eta_wo_spill
 
     # Forward efficiency: does/should not include window loss
     # because it is defined as how much power out of

tests/test_physics.py

@@ -3,7 +3,7 @@
 
 
 def test_johnson_nyquist_psd():
-    expected = 1.34778 * 10 ** -21
+    expected = 1.34778 * 10**-21
     output = physics.johnson_nyquist_psd(1e11, 100)
     assert isclose(output, expected, rel_tol=0.01)