Skip to content
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
18 changes: 9 additions & 9 deletions examples/Notebooks/flopy3_mf6_A_simple-model.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -137,9 +137,9 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"Now that the overall simulation is setup, we can focus on building the groundwater flow model. The groundwater flow model will be built by adding packages to it that describe the model characteristics.\n",
"Now that the overall simulation is set up, we can focus on building the groundwater flow model. The groundwater flow model will be built by adding packages to it that describe the model characteristics.\n",
"\n",
"Define the discretization of the model. All layers are given equal thickness. The `bot` array is build from the `Hlay` values to indicate top and bottom of each layer, and `delrow` and `delcol` are computed from model size `L` and number of cells `N`. Once these are all computed, the Discretization file is built."
"Define the discretization of the model. All layers are given equal thickness. The `bot` array is build from `H` and the `Nlay` values to indicate top and bottom of each layer, and `delrow` and `delcol` are computed from model size `L` and number of cells `N`. Once these are all computed, the Discretization file is built."
]
},
{
Expand All @@ -149,7 +149,7 @@
"outputs": [],
"source": [
"# Create the discretization package\n",
"bot = np.linspace(-H/Nlay,-H,Nlay) \n",
"bot = np.linspace(-H/Nlay, -H, Nlay) \n",
"delrow = delcol = L/(N-1)\n",
"dis = flopy.mf6.modflow.mfgwfdis.ModflowGwfdis(gwf, pname='dis', nlay=Nlay, nrow=N, ncol=N,\n",
" delr=delrow,delc=delcol,top=0.0,\n",
Expand Down Expand Up @@ -191,14 +191,14 @@
"# must be entered as a tuple as the first entry.\n",
"# Remember that these must be zero-based indices!\n",
"chd_rec = []\n",
"chd_rec.append(((0, int(N / 4), int(N / 4)), 90.0))\n",
"chd_rec.append(((0, int(N / 4), int(N / 4)), h2))\n",
"for layer in range(0, Nlay):\n",
" for row_col in range(0, N):\n",
" chd_rec.append(((layer, row_col, 0), 100.0))\n",
" chd_rec.append(((layer, row_col, N - 1), 100.0))\n",
" chd_rec.append(((layer, row_col, 0), h1))\n",
" chd_rec.append(((layer, row_col, N - 1), h1))\n",
" if row_col != 0 and row_col != N - 1:\n",
" chd_rec.append(((layer, 0, row_col), 100.0))\n",
" chd_rec.append(((layer, N - 1, row_col), 100.0))\n",
" chd_rec.append(((layer, 0, row_col), h1))\n",
" chd_rec.append(((layer, N - 1, row_col), h1))\n",
"chd = flopy.mf6.modflow.mfgwfchd.ModflowGwfchd(gwf, pname='chd', maxbound=len(chd_rec), \n",
" stress_period_data=chd_rec, save_flows=True)"
]
Expand Down Expand Up @@ -1041,4 +1041,4 @@
},
"nbformat": 4,
"nbformat_minor": 2
}
}