diff --git a/tools/plotting/fastnnlo_scaleunc.py b/tools/plotting/fastnnlo_scaleunc.py
index f9acfefa58dddde5883c3567259f8c748d162f3a..3f43ab3d5450c638d908f841b9f80a5741332939 100755
--- a/tools/plotting/fastnnlo_scaleunc.py
+++ b/tools/plotting/fastnnlo_scaleunc.py
@@ -156,8 +156,8 @@ def plotting(x_axis, xmin, xmax, xs_all, rel_scale_unc, abs_scale_unc, dxsr_cn,
     #        ax1.get_xaxis().set_minor_formatter(NullFormatter())
     #    ax1.set_xscale('log', nonposx='clip')
     #    ax1.set_yscale('log', nonposy='clip')
-    ax1.set_xscale('log', nonposx='clip') # Will become nonpositive= in matplotlib 3.3.x
-    ax1.set_yscale('log', nonposy='clip') # Will become nonpositive= in matplotlib 3.3.x
+    ax1.set_xscale('log', nonpositive='clip') # Will become nonpositive= in matplotlib 3.3.x
+    ax1.set_yscale('log', nonpositive='clip') # Will become nonpositive= in matplotlib 3.3.x
     # Set label on x axis on coupled axis from lower plot
     #        ax1.set_xlabel(r'%s' %xlabel, horizontalalignment='right', x=1.0, verticalalignment='top', y=1.0)
     ax1.set_ylabel(r'%s' % ylabel, horizontalalignment='right', x=1.0,
@@ -445,6 +445,7 @@ def main():
 
         # Creating x-axis
         bin_bounds = np.array(fnlo.GetObsBinsBounds(0))
+        #        bin_bounds = np.array(fnlo.GetDim1BinBounds(0)) Mods needed for 2D tables ...
         if verb:
             print('[fastnnlo_scaleunc]: bin_bounds.T: \n', bin_bounds.T, '\n')
             print('[fastnnlo_scaleunc]: bin_bounds.flatten()',