From e0927c6423fd4819d70dfe389cc18c679110783c Mon Sep 17 00:00:00 2001 From: craig-warren Date: Thu, 8 Dec 2016 08:32:34 +0000 Subject: [PATCH] Fixed dispersion analysis for 2D simulations with non-cubic cells. --- gprMax/grid.py | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) diff --git a/gprMax/grid.py b/gprMax/grid.py index 550e0581..b35e8d88 100644 --- a/gprMax/grid.py +++ b/gprMax/grid.py @@ -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