import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import griddata
from mpl_toolkits.mplot3d import Axes3D

# データ取得
data_x, data_y, data_phi = np.loadtxt("./output.dat", comments='!', unpack=True)

# grid data 生成
X1,X2 = np.meshgrid(np.unique(data_x),np.unique(data_y))
PHI = griddata((data_x, data_y), data_phi, (X1, X2))

# 図、枠（軸）の設定
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
#ax = Axes3D(fig)
# Axesのタイトルを設定
ax.set_title("Laplace", size = 20)

# 軸ラベルを設定
ax.set_xlabel("x", size = 14, color = "k")
ax.set_ylabel("y", size = 14, color = "k")
ax.set_zlabel(r"$\phi$", size = 14, color = "k")

#surface = ax.plot_surface(X1,X2,PHI)
#fig.colorbar(surface)
surface = ax.plot_wireframe(X1,X2,PHI)


plt.show()