diff --git a/orthofinder/scripts/scnorm.py b/orthofinder/scripts/scnorm.py
index 4b482ab74b1ee9a220f84b8a5910f9cbc7517735..6d8b4484906831dab0cf17fcc221ff8e9a84bca6 100644
--- a/orthofinder/scripts/scnorm.py
+++ b/orthofinder/scripts/scnorm.py
@@ -31,6 +31,7 @@ import numpy as np
 import matplotlib as mpl                        # install
 mpl.use('Agg')
 import matplotlib.pyplot as plt                 # install
+import sys
 
 
 counter=0
@@ -103,13 +104,13 @@ class scnorm:
     @staticmethod
     def DrawLogLinear(params, x):
         global counter
+        f = []
         global f
         nbvals=10000
         maxBS=2500
         minBS=0
         x = np.linspace(minBS, maxBS, nbvals)
-        plt.plot(x, scnorm.loglinear(x, *params), label='species '+str(counter))
-        f = plt.plot(x, scnorm.loglinear(x, *params), label='species '+str(counter))
+        f.append(plt.plot(x, scnorm.loglinear(x, *params), label='species '+str(counter)))
         plt.legend()
         plt.xlabel('Bits Scores')
         plt.ylabel('Normalized Scores')