subgrid only geometry view for voxels

user_models/sub-gridding/subgrid_basic.py

@@ -81,21 +81,29 @@ gv = gprMax.GeometryView(p1=(0, 0, 0),
                          output_type='f',
                          multi_grid=True)
 
-# create a geometry view of the sub grid
+# create a geometry view of the sub grid only
 gv_sg = gprMax.GeometryView(p1=sg_p0,
                          p2=sg_p1,
                          dl=(1e-3, 1e-3, 1e-3),
                          filename=fn.with_suffix('').parts[-1] + '_only',
                          output_type='f')
 
+# create a geometry view of the sub grid only
+gv_sg_normal = gprMax.GeometryView(p1=sg_p0,
+                         p2=sg_p1,
+                         dl=(1e-3, 1e-3, 1e-3),
+                         filename=fn.with_suffix('').parts[-1] + '_only',
+                         output_type='n')
+
 # create a geometry view of the main grid and the sub grid stitched together
 gv_normal = gprMax.GeometryView(p1=(0, 0, 0),
                          p2=domain.props.p1,
                          dl=dl,
                          filename=fn.with_suffix('').parts[-1] + '_voxels',
                          output_type='n')
-#scene.add(gv)
+scene.add(gv)
 sg.add(gv_sg)
+sg.add(gv_sg_normal)
 #scene.add(gv_normal)
 
 gprMax.run(scenes=[scene], n=1, geometry_only=False, outputfile=fn, subgrid=True, autotranslate=True)