Skip to content

Commit

Permalink
fix(flopy3_MT3DMS_examples.ipynb): match constant heads to original p…
Browse files Browse the repository at this point in the history
…roblem (#809)

Problem 10 in the MT3DMS manual specifies constant heads along the top and bottom boundaries, not just the sides
  • Loading branch information
emorway-usgs committed Feb 11, 2020
1 parent 207cd1e commit 0d15251
Showing 1 changed file with 8 additions and 1 deletion.
9 changes: 8 additions & 1 deletion examples/Notebooks/flopy3_MT3DMS_examples.ipynb
Expand Up @@ -1900,8 +1900,15 @@
" botm=botm,\n",
" perlen=perlen_mf)\n",
" ibound = np.ones((nlay, nrow, ncol), dtype=np.int)\n",
"\n",
" # left side\n",
" ibound[:, :, 0] = -1\n",
" # right side\n",
" ibound[:, :, -1] = -1\n",
" # top\n",
" ibound[:, 0, :] = -1\n",
" # bottom\n",
" ibound[:, -1, :] = -1\n",
"\n",
" f = open(os.path.join(datadir, 'p10shead.dat'))\n",
" s0 = np.empty((nrow * ncol), dtype=np.float)\n",
Expand Down Expand Up @@ -2212,7 +2219,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
"version": "3.5.3"
}
},
"nbformat": 4,
Expand Down

0 comments on commit 0d15251

Please sign in to comment.