From 4de2280c487dcb54c876e1a7ce611f9726249cfc Mon Sep 17 00:00:00 2001
From: Craig Warren <craig.warren@northumbria.ac.uk>
Date: Wed, 15 Dec 2021 13:45:31 +0000
Subject: [PATCH] Correct calculation of Er for Debye materials, used for
 estimating numerical dispersion.

---
 gprMax/materials.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/gprMax/materials.py b/gprMax/materials.py
index 835eacce..55f63b55 100644
--- a/gprMax/materials.py
+++ b/gprMax/materials.py
@@ -194,7 +194,7 @@ class DispersiveMaterial(Material):
         er = self.er
 
         w = 2 * np.pi * freq
-        er += self.se / (w * config.sim_config.em_consts['e0'])
+        er += self.se / (1j * w * config.sim_config.em_consts['e0'])
         if 'debye' in self.type:
             for pole in range(self.poles):
                 er += self.deltaer[pole] / (1 + 1j * w * self.tau[pole])