figure 9 quadratic scatter plots

modelsusceptlown.py

@@ -58,11 +58,12 @@ def plot_si_diags(ax, s, data, alphax, alphay, xthresh, ythresh, cell_name):
     l = linregress(six, siy)
     x = np.linspace(0, 10, 10)
     ax.set_visible(True)
+    ax.set_aspect('equal')
     ax.set_title(f'$c$={100*alphay:g}\\,\\%', fontsize='medium')
     ax.plot(x, x, **s.lsLine)
     ax.plot(x, l.slope*x + l.intercept, **s.lsGrid)
     ax.axhline(ythresh, **s.lsLine)
-    ax.axvline(xthresh, 0, 0.5, **s.lsLine)
+    ax.axvline(xthresh, 0, 0.65, **s.lsLine)
     if alphax == 0:
         mask = datax['triangle'] > 0.5
         ax.plot(six[mask], siy[mask], zorder=30, label='strong', **s.psA1m)
@@ -137,8 +138,8 @@ if __name__ == '__main__':
     xthresh = 1.2
     ythresh = 1.8
     s = plot_style()
-    fig, axs = plt.subplots(3, 4, cmsize=(s.plot_width, 0.6*s.plot_width),
-                            height_ratios=[3, 3, 0, 2])
+    fig, axs = plt.subplots(3, 4, cmsize=(s.plot_width, 0.62*s.plot_width),
+                            height_ratios=[3, 3, 0, 3])
     fig.subplots_adjust(leftm=7, rightm=9, topm=2, bottomm=4,
                         wspace=1, hspace=0.6)
     for ax in axs.flat: