diff --git a/bfps/DNS.py b/bfps/DNS.py
index 6cc994cac2ffbd3a0e5812c2ddeeca5fd98fa334..909b5bb42accfc73b5325d6b0e3b40d8f3c7b86a 100644
--- a/bfps/DNS.py
+++ b/bfps/DNS.py
@@ -389,8 +389,8 @@ class DNS(_code):
                 self.statistics['kMeta' + suffix] *= 0.8
         self.statistics['Lint'] = ((np.pi /
                                     (2*self.statistics['Uint']**2)) *
-                                   np.nansum(self.statistics['energy(k)'] /
-                                                self.statistics['kshell']))
+                                   np.sum(self.statistics['energy(k)'][1:-2] /
+                                          self.statistics['kshell'][1:-2]))
         self.statistics['Re'] = (self.statistics['Uint'] *
                                  self.statistics['Lint'] /
                                  self.parameters['nu'])