From 172278a6f04e519b76586dd4f3e6261323a42603 Mon Sep 17 00:00:00 2001 From: craig-warren Date: Tue, 13 Oct 2015 21:28:37 +0300 Subject: [PATCH] Corrected bug where rough surfaces would not get created. --- gprMax/input_cmds_geometry.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/gprMax/input_cmds_geometry.py b/gprMax/input_cmds_geometry.py index 6599925d..7f4c75dd 100644 --- a/gprMax/input_cmds_geometry.py +++ b/gprMax/input_cmds_geometry.py @@ -1408,13 +1408,15 @@ def process_geometrycmds(geometry, G): for k in range(volume.zs, volume.zf + 1): if volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 1: numID = numIDx = numIDy = numIDz = volume.fractalvolume[i - volume.xs, j - volume.ys, k - volume.zs] + build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID) elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 2: waternumID = next(x.numID for x in G.materials if x.ID == 'water') numID = numIDx = numIDy = numIDz = waternumID + build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID) elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 3: grassnumID = next(x.numID for x in G.materials if x.ID == 'grass') numID = numIDx = numIDy = numIDz = grassnumID - build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID) + build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID) else: if volume.nbins == 1: