def plot_Reynolds_number(Fr, ReyNum, edge):
# figr = Fr.plot_fracture("complete", "footPrint")
# ax = figr.axes[0]
figr = plt.figure()
ax = figr.add_subplot(111)
ReMesh = np.resize(ReyNum[edge, :], (Fr.mesh.ny, Fr.mesh.nx))
x = np.linspace(-Fr.mesh.Lx, Fr.mesh.Lx, Fr.mesh.nx)
y = np.linspace(-Fr.mesh.Ly, Fr.mesh.Ly, Fr.mesh.ny)
xv, yv = np.meshgrid(x, y)
# cax = ax.contourf(xv, yv, ReMesh, levels=[0, 100, 2100, 10000])
cax = ax.matshow(ReMesh)
figr.colorbar(cax)
plt.title("Reynolds number")
plt.show()
return figr
#-----------------------------------------------------------------------------------------------------------------------
评论列表
文章目录