Python 3 compatibility

get_pdfseries.py

@@ -89,8 +89,8 @@ if __name__ == "__main__":
         pdf1 = norm(loc=x1s[i], scale=ss)
         pdf2 = norm(loc=x2s[i], scale=ss)
         samp = np.r_[
-                    pdf1.rvs(ngrid / 2),
-                    pdf2.rvs(ngrid / 2)
+                    pdf1.rvs(int(ngrid / 2)),
+                    pdf2.rvs(int(ngrid / 2))
                 ]
         timeseries_grid[:, i] = samp
         kde = gaussian_kde(samp)

plot_results.py

@@ -89,7 +89,8 @@ if __name__ == "__main__":
     # set up axes
     # for analysis of PDFs
     ax_a = fig.add_axes([lm1, bm1, wd1, ht1])
-    ax_b = ax_a.twinx()
+    #ax_b = ax_a.twinx()
+    ax_b = fig.add_axes([lm1, bm2, wd1, ht2])
     ax_b.set_position([lm1, bm2, wd1, ht2])
     ax_c = fig.add_axes([lm1, bm3, wd1, ht1])
     cx_a = fig.add_axes([lm2, bm1, wd2, ht1])