Add plot of surface in test AHF Kerr

msc/test_AHF_Kerr_a00.py

@@ -38,3 +38,15 @@ ahf = AHF(lmax = 4, mass_tol = 1e-4,
           outputdir = outputdir)
 ahf.run([dfile])
 
+
+# Horizon surface plot
+# --------------------
+
+x,y,z = np.loadtxt('test_AHF_Kerr_a00/horizon00_surface.txt',unpack=True)
+
+import matplotlib.pyplot as plt
+
+fig = plt.figure()
+ax = fig.add_subplot(projection='3d')
+ax.plot_surface(x.reshape(8,16), y.reshape(8,16), z.reshape(8,16), alpha=0.8)
+plt.show()