35
35
tdis_rc = [[1.0 , 1 , 1.0 ]]
36
36
porosity = 0.1
37
37
38
+
38
39
def chdhead (x ):
39
40
return x * 10.0 / 100.0
40
41
41
42
42
43
def get_tri (workspace , targets ) -> Triangle :
43
44
workspace .mkdir (exist_ok = True )
44
- tri = Triangle (angle = angle , maximum_area = max_area , model_ws = workspace , exe_name = targets .triangle )
45
+ tri = Triangle (
46
+ angle = angle ,
47
+ maximum_area = max_area ,
48
+ model_ws = workspace ,
49
+ exe_name = targets .triangle ,
50
+ )
45
51
tri .add_polygon (active_domain )
46
52
tri .build ()
47
53
return tri
@@ -53,9 +59,7 @@ def build_gwf_sim(name, ws, targets):
53
59
sim = flopy .mf6 .MFSimulation (
54
60
sim_name = name , version = "mf6" , exe_name = targets .mf6 , sim_ws = ws
55
61
)
56
- tdis = flopy .mf6 .ModflowTdis (
57
- sim , time_units = "DAYS" , perioddata = tdis_rc
58
- )
62
+ tdis = flopy .mf6 .ModflowTdis (sim , time_units = "DAYS" , perioddata = tdis_rc )
59
63
gwf = flopy .mf6 .ModflowGwf (sim , modelname = gwfname , save_flows = True )
60
64
ims = flopy .mf6 .ModflowIms (
61
65
sim ,
@@ -81,7 +85,10 @@ def build_gwf_sim(name, ws, targets):
81
85
cell2d = cell2d ,
82
86
)
83
87
npf = flopy .mf6 .ModflowGwfnpf (
84
- gwf , xt3doptions = [(True )], save_specific_discharge = True , save_saturation = True
88
+ gwf ,
89
+ xt3doptions = [(True )],
90
+ save_specific_discharge = True ,
91
+ save_saturation = True ,
85
92
)
86
93
ic = flopy .mf6 .ModflowGwfic (gwf )
87
94
chdlist = []
0 commit comments