update example 02

GeoStat-Framework · Aug 16, 2019 · c3107f9 · c3107f9
1 parent e057571
commit c3107f9
Showing 1 changed file with 2 additions and 8 deletions.
diff --git a/examples/02_pump_test_het_3D.py b/examples/02_pump_test_het_3D.py
@@ -12,7 +12,7 @@
 model.gli.generate("radial", dim=3, angles=64, rad_out=100, z_size=-10)
 model.gli.add_polyline("pwell", [[0, 0, 0], [0, 0, -10]])
 model.msh.generate(
-    "radial", dim=3, angles=64, rad=np.arange(101), z_arr=-np.arange(11),
+    "radial", dim=3, angles=64, rad=np.arange(101), z_arr=-np.arange(11)
 )
 cond = np.exp(srf.mesh(model.msh))
 mpd = MPD(model.task_id)
@@ -36,20 +36,14 @@
     GEO_TYPE=["POLYLINE", "pwell"],
     DIS_TYPE=["CONSTANT_NEUMANN", -1.0e-3],
 )
-model.ic.add_block(  # set the initial condition
-    PCS_TYPE="GROUNDWATER_FLOW",
-    PRIMARY_VARIABLE="HEAD",
-    GEO_TYPE="DOMAIN",
-    DIS_TYPE=["CONSTANT", 0.0],
-)
 model.mmp.add_block(  # permeability, storage and porosity
     GEOMETRY_DIMENSION=3,
     STORAGE=[1, 1.0e-4],
     PERMEABILITY_DISTRIBUTION=model.task_id + ".mpd",
     POROSITY=0.2,
 )
 model.num.add_block(  # numerical solver
-    PCS_TYPE='GROUNDWATER_FLOW',
+    PCS_TYPE="GROUNDWATER_FLOW",
     LINEAR_SOLVER=[2, 5, 1.0e-14, 1000, 1.0, 100, 4],
 )
 model.out.add_block(  # set the outputformat