Fixed dispersion analysis for 2D simulations with non-cubic cells.

这个提交包含在:
craig-warren
2016-12-08 08:32:34 +00:00
父节点 a99c0ed14f
当前提交 e0927c6423

查看文件

@@ -224,7 +224,15 @@ def dispersion_analysis(G):
minwavelength = minvelocity / results['maxfreq']
# Maximum spatial step
delta = max(G.dx, G.dy, G.dz)
if G.dimension == '3D':
delta = max(G.dx, G.dy, G.dz)
elif G.dimension == '2D':
if G.nx == 1:
delta = max(G.dy, G.dz)
elif G.ny == 1:
delta = max(G.dx, G.dz)
elif G.nz == 1:
delta = max(G.dx, G.dy)
# Courant stability factor
S = (c * G.dt) / delta