pep8

wafo/misc.py

@@ -2487,7 +2487,8 @@ def _make_bars(limits, bin_):
 
 
 # pylint: disable=redefined-builtin
-def _histogram(data, bins=None, range=None, normed=False, weights=None, density=None):  # @ReservedAssignment
+def _histogram(data, bins=None, range=None, normed=False, weights=None,  # @ReservedAssignment
+               density=None):
     """
     Example
     -------

wafo/spectrum/core.py

@@ -3722,8 +3722,8 @@ class SpecData1D(PlotData):
         # Cov(Hm0,Tm01) Cov(Hm0,Tm02) Cov(Tm01,Tm02)
         S0 = r_[2. / (sqrt(m[0]) * m[1]) * (mij[0] - m[0] * mij[1] / m[1]),
                 1. / sqrt(m[2]) * (mij[0] / m[0] - mij[2] / m[2]),
-                1. / (2 * m[1]) * sqrt(m[0] / m[2]) * (mij[0] / m[0] - mij[2] /
+                1. / (2 * m[1]) * sqrt(m[0] / m[2])
-                                                       m[2] - mij[1] / m[1] + m[0] * mij[3] / (m[1] * m[2]))]
+                * (mij[0] / m[0] - mij[2] / m[2] - mij[1] / m[1] + m[0] * mij[3] / (m[1] * m[2]))]
 
         R1 = ones((15, 15))
         R1[:, :] = nan

wafo/spectrum/models.py

@@ -1313,8 +1313,7 @@ class Wallop(Bretschneider):
             wp = 2. * pi / Tp
             kp = w2k(wp, 0, inf)[0]  # wavenumber at peak frequency
             Lp = 2. * pi / kp  # wave length at the peak frequency
-            N = np.abs((log(2. * pi ** 2.) + 2 * log(Hm0 / 4) -
+            N = np.abs((log(2. * pi ** 2.) + 2 * log(Hm0 / 4) - 2.0 * log(Lp)) / log(2))
-                     2.0 * log(Lp)) / log(2))
 
         super(Wallop, self).__init__(Hm0, Tp, N, M, chk_seastate)
 
@@ -2135,5 +2134,6 @@ def main():
     else:
         test_docstrings()
 
+
 if __name__ == '__main__':
     main()