diff --git a/00_model.ipynb b/00_model.ipynb index 709b9ce..322090f 100755 --- a/00_model.ipynb +++ b/00_model.ipynb @@ -7,11 +7,12 @@ "source": [ "# Background\n", "\n", - "This series of tutorials are based on the is based on the GMDSI Open Loop Low Temperature Geothermal System (non-python) tutorials, which can be found here: https://gmdsi.org/blog/gwe_slideshow/\n", + "This series of tutorials are based on the GMDSI Open Loop Low Temperature Geothermal System (non-python) tutorials, which can be found here:\n", "\n", - "The purpose of these `notebooks` is to demonstrate how to:\n", + "https://gmdsi.org/blog/gwe_slideshow/\n", "\n", - " - go about programatically settting up a `PEST++` interface with an unstructured grid model;\n", + "The purpose of these `notebooks` is to demonstrate how to:\n", + " - go about programmatically setting up a `PEST++` interface with an unstructured grid model;\n", " - parameterize parameter fields with non-stationary geostatistics informed by conceptual knowledge; and\n", " - undertake surrogate-based data assimilation and predictive uncertainty analysis using data space inversion" ] @@ -43,9 +44,9 @@ "\n", "The image below is extracted from the [slideshow](https://gmdsi.org/wp-content/uploads/2025/04/gwf_gwe_model.pptx) file that accompanies the original GMDSI tutorial. It shows the model grid and boundary conditions.\n", "\n", - "The model represents a shallow, channelized alluviual deposits close to a river (North boundary), which pinch out where bedrock outcrops to the South. The system is represented using a single layer model. North boundary is simulated using the RIV package. East and Western boundaries are specified as general head boundaries. \n", + "The model represents a shallow, channelized alluvial deposits close to a river (North boundary), which pinch out where bedrock outcrops to the South. The system is represented using a single layer model. North boundary is simulated using the RIV package. East and Western boundaries are specified as general head boundaries. \n", "\n", - "The model is used to simulate a (past) historical period and a (future) forecast period. In the forecast period an open loop shalow geothermal system is active. The flow model is simulated as steady-state for both historical and forecast periods. The heat transport model is transient. In the foreacst period it simulates periods of seasonaly varying injection temperature. The background temperature is 16°C.\n", + "The model is used to simulate a (past) historical period and a (future) forecast period. In the forecast period an open loop shallow geothermal system is active. The flow model is simulated as steady-state for both historical and forecast periods. The heat transport model is transient. The forecast period simulates periods of seasonally varying injection temperature. The background temperature is 16°C.\n", "\n", "\n", "![image info](./data/Pmodel.png)\n" @@ -74,7 +75,7 @@ "id": "ae9cdc79", "metadata": {}, "source": [ - "Run the next cell to call a helper function we have prepared, just to amke sure the correct `modflow6` and `pest` executables are in the model filder:" + "Run the next cell to call a helper function we have prepared, just to make sure the correct `modflow6` and `pest` executables are in the model folder:" ] }, { @@ -93,7 +94,7 @@ "metadata": {}, "source": [ "## Run the model once\n", - "Unusualy for our tutorials, this model actualy takes a while to run. We need to run it once to make sure all the model output files are generated.\n", + "Unusually for our tutorials, this model actually takes a while to run. We need to run it once to make sure all the model output files are generated.\n", "\n", "### Warning: this might take a few minutes. \n", "For context, it takes ~5min on a MacBook Pro.\n" @@ -251,7 +252,7 @@ "id": "2d371601", "metadata": {}, "source": [ - "A look at everyones favourite paramter..." + "A look at everyone's favorite parameter..." ] }, { @@ -261,7 +262,7 @@ "metadata": {}, "outputs": [], "source": [ - "gwf.npf.k.plot()" + "gwf.npf.k.plot(colorbar=True)" ] }, { @@ -271,13 +272,13 @@ "source": [ "# Conceptual model\n", "\n", - "A key aspect of this case is that, we have some conceptual knowledge about the geostatistics of the hydraulic conducvitiy. \n", + "A key aspect of this case is that, we have some conceptual knowledge about the geostatistics of the hydraulic conductivity. \n", "\n", "We sat down with the project geologist and translated their conceptual knowledge on \"how K is expected to vary in space\" into a set of hyperparameters assigned to pilot points. These include:\n", "\n", " - mean K value\n", " - variance\n", - " - spatial correlatin length\n", + " - spatial correlation length\n", " - anisotropy ratio\n", " - bearing\n", "\n", @@ -301,9 +302,9 @@ "id": "8de5b284", "metadata": {}, "source": [ - "In the plot below, at each pilot point, we have plotted the mean K value (color of the point), the line of bearing, the length of which reflects the correlation length (longer lines means greater correlation), and the perpendicular line showing the transversal bearing&correlation length.\n", + "In the plot below, at each pilot point, we have plotted the mean K value (color of the point), the line of bearing, the length of which reflects the correlation length (longer lines means greater correlation), and the perpendicular line showing the transversal bearing and correlation length.\n", "\n", - "Thi sinformation will be used later to generate random feilds of K which respect these patterns of heterogeinty. Strap in!" + "This information will be used later to generate random fields of K which respect these patterns of heterogeneity. Strap in!" ] }, { @@ -353,11 +354,19 @@ "\n", "fig.tight_layout();" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "52fb9a67-d853-46a3-8f33-d088a5591243", + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { "kernelspec": { - "display_name": "gwe", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, diff --git a/01_setup_pstfrom.ipynb b/01_setup_pstfrom.ipynb index d41eaa0..d7a1a03 100755 --- a/01_setup_pstfrom.ipynb +++ b/01_setup_pstfrom.ipynb @@ -27,9 +27,11 @@ "source": [ "# Getting started\n", "\n", - "These notebooks will be less verbose than other GMDSI tutorial notebooks. We assume the reader is already familiar with many of the topics. We focus specificaly on aspects that relate to:\n", + "These notebooks will be less verbose than other GMDSI tutorial notebooks. We assume the reader is already familiar with many of the topics. \n", + "\n", + "We focus specifically on aspects that relate to:\n", "- `PstFrom` with a `mf6` `disv` grid\n", - "- use of hyper parameters\n", + "- use of hyperparameters with `pypestutils`\n", "- passing in prior knowledge to hyperparameters\n", "- (subsequent notebooks) use of DSI\n", "\n", @@ -60,8 +62,10 @@ { "cell_type": "code", "execution_count": null, - "id": "cd6a229b", - "metadata": {}, + "id": "e781d5bb-56a0-4889-b6e2-fe29b88ed21b", + "metadata": { + "scrolled": true + }, "outputs": [], "source": [ "# load simulation\n", @@ -73,29 +77,21 @@ { "cell_type": "code", "execution_count": null, - "id": "f60f390d", - "metadata": {}, - "outputs": [], - "source": [ - "tdis = sim.tdis\n", - "perioddata = tdis.perioddata.get_data()\n", - "\n", - "perioddata['nstp'] = 1\n", - "\n", - "perioddata\n", - "tdis.perioddata.set_data(perioddata)\n", - "tdis.write()" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "f0eef6a3", - "metadata": {}, + "id": "a2348537-8f67-45f8-af43-2338e980a79d", + "metadata": { + "scrolled": true + }, "outputs": [], "source": [ - "# run the model once to make sure it works\n", - "pyemu.os_utils.run(\"mf6\",cwd=tmp_d)" + "# because this model is a little slow... only run if starting fresh (hopefully we ran in the last notebook?):\n", + "if 'mfsim.lst' in os.listdir(tmp_d):\n", + " with open(os.path.join(tmp_d, 'mfsim.lst'), 'r') as fp:\n", + " lines = fp.readlines()\n", + "else:\n", + " lines = ['no']\n", + "if \"Normal termination\" not in lines[-1]:\n", + " # run the model once to make sure it works\n", + " pyemu.os_utils.run(\"mf6\", cwd=tmp_d)" ] }, { @@ -103,7 +99,7 @@ "id": "8ca97268", "metadata": {}, "source": [ - "This model is a `DISV` grid. `PstFrom` is going to require the `flopy` model grid object to setup pilot points and spatialy varying covariance:" + "This model is a `DISV` grid. `PstFrom` is going to require the `flopy` model grid object to setup pilot points and spatially varying covariance:" ] }, { @@ -153,7 +149,7 @@ "id": "e7fa6d8a", "metadata": {}, "source": [ - "We are going to keep things super simple for this tutorial. We are only going to parameterize hydrualic conductivity. As we are all sofisticated, we understand that in a real-world application other parameters and boundary conditions would likely also be important aspects to consider for data assimilation and uncertianty analysis...." + "We are going to keep things super simple for this tutorial. We are only going to parameterize hydraulic conductivity. As we are all sofisticated, we understand that in a real-world application other parameters and boundary conditions would likely also be important aspects to consider for data assimilation and uncertainty analysis...." ] }, { @@ -175,7 +171,7 @@ "id": "9da02034", "metadata": {}, "source": [ - "Unfortunatley `flopy` doesnt write tidy model input files...so we need to fix them..." + "Unfortunately `flopy` doesn't write tidy model input files...so we need to fix them..." ] }, { @@ -208,7 +204,7 @@ "id": "5666c6b4", "metadata": {}, "source": [ - "No one specifgied the `idomain` in the original model setup, so lets just create a \"zone array\". `PstFrom` expects this when we setup pilot points and so on later. Note that the shape of `ib` is the same as the shape of the `k`:" + "No one specified the `idomain` in the original model setup, so lets just create a \"zone array\". `PstFrom` expects this when we setup pilot points and so on later. Note that the shape of `ib` is the same as the shape of the `k`:" ] }, { @@ -247,7 +243,7 @@ "id": "980b482b", "metadata": {}, "source": [ - "However, `PstFrom` and `pypestutils` expect a strict format in terms of column anmes and information. So first we need to make that. A pilot point file must have the following columns: `['name','zone','x','y','parval1']`" + "However, `PstFrom` and `pypestutils` expect a strict format in terms of column names and information. So first we need to make that. A pilot point file must have the following columns: `['name','zone','x','y','parval1']`" ] }, { @@ -291,7 +287,7 @@ "id": "dca6384e", "metadata": {}, "source": [ - "We spatialy varying mean values of K in the conceptaul pilot points..." + "We spatially varying mean values of K in the conceptual pilot points..." ] }, { @@ -358,7 +354,7 @@ "id": "0d530d9c", "metadata": {}, "source": [ - "Now, setup a goestatistical structure to pass to `PstFrom`. This is the geostatiscs for the \"hyper parameter\". Inception much...\n", + "Now, setup a geostatistical structure to pass to `PstFrom`. This is the geostatiscs for the \"hyper parameter\". Inception much...\n", "\n", "Lets just use the median from the conceptual pilot points:" ] @@ -433,7 +429,7 @@ " pp_options={\"prep_hyperpars\":True, \"pp_space\":ppfname}\n", "```\n", "\n", - "Where we have specfied that `PstFrom` shoudl setup hyper parameters, and use the `ppfname` file for pilot point locations." + "Where we have specified that `PstFrom` should setup hyper parameters, and use the `ppfname` file for pilot point locations." ] }, { @@ -565,7 +561,7 @@ "id": "be4a3d0b", "metadata": {}, "source": [ - "OK, so lets go through each of the hyper parameters an dparameterize them using pilot points...same as before..." + "OK, so lets go through each of the hyper parameters and parameterize them using pilot points...same as before..." ] }, { @@ -573,7 +569,7 @@ "id": "f5bd86e6", "metadata": {}, "source": [ - "Lets start with anisotropy. We are going to make these hyperparameters all \"additive\" type paraemters. " + "Lets start with anisotropy. We are going to make these hyperparameters all \"additive\" type parameters. " ] }, { @@ -601,22 +597,6 @@ "ppdf" ] }, - { - "cell_type": "code", - "execution_count": null, - "id": "3a0ce04f", - "metadata": {}, - "outputs": [], - "source": [] - }, - { - "cell_type": "code", - "execution_count": null, - "id": "08794faa", - "metadata": {}, - "outputs": [], - "source": [] - }, { "cell_type": "code", "execution_count": null, @@ -624,12 +604,8 @@ "metadata": {}, "outputs": [], "source": [ - "\n", "afile = tag+'.aniso.dat'\n", - "\n", - "\n", "tidy_array(os.path.join(template_ws,afile))\n", - "\n", "atag = afile.split('.')[0].replace(\"_\",\"-\")+\"-aniso\"\n", "_df = pf.add_parameters(afile,\n", " par_type=\"pilotpoints\",\n", @@ -648,14 +624,6 @@ "_ = pf.add_observations(afile, prefix=atag, obsgp=atag)" ] }, - { - "cell_type": "code", - "execution_count": null, - "id": "4567fbd4", - "metadata": {}, - "outputs": [], - "source": [] - }, { "cell_type": "markdown", "id": "03471fa7", @@ -732,9 +700,6 @@ "metadata": {}, "outputs": [], "source": [ - "\n", - "\n", - "\n", "afile = tag+'.corrlen.dat'\n", "tidy_array(os.path.join(template_ws,afile))\n", "atag = afile.split('.')[0].replace(\"_\",\"-\")+\"-corrlen\"\n", @@ -782,7 +747,7 @@ "source": [ "# Check the prior K fields\n", "\n", - "An extremly useful check can be done now. Note that we have not yet added a model run to the `PstFrom`. All that is in the forward run workflow up to now is the interpolation from pilot points to the model grid. In other words, the interpolation is our \"forward run\" at the moment. \n", + "An extremely useful check can be done now. Note that we have not yet added a model run to the `PstFrom`. All that is in the forward run workflow up to now is the interpolation from pilot points to the model grid. In other words, the interpolation is our \"forward run\" at the moment. \n", "\n", "Build the pest control file and forward run .py file to see:" ] @@ -812,11 +777,11 @@ "id": "127d561c", "metadata": {}, "source": [ - "This is super powerful. We can generate a prior ensemble of pilot point values, run that ensemble and collate all the generater model input parameter fields. We can then check them and make sure they make sense and look pretty :) \n", + "This is super powerful. We can generate a prior ensemble of pilot point values, run that ensemble and collate all the generated model input parameter fields. We can then check them and make sure they make sense and look pretty :) \n", "\n", "All we need to do is generate the ensemble and run it once. Because we are not running modflow, this will be super fast!\n", "\n", - "Lets just make sure wverything is working first:" + "Lets just make sure everything is working first:" ] }, { @@ -986,7 +951,7 @@ "id": "639b24a1", "metadata": {}, "source": [ - "Read in the results of the `pestpp-ies` prior monte carlo. We can use some the `Pst` inbuilt helpers for this. Start b reading the .pst file form the mater dir:" + "Read in the results of the `pestpp-ies` Prior Monte Carlo. We can use some the `Pst` inbuilt helpers for this. Start by reading the .pst file form the mater dir:" ] }, { @@ -1004,7 +969,7 @@ "id": "7ba450ff", "metadata": {}, "source": [ - "Now, if parmeter or observation ensemble files are avialable in the folder, `pyemu` will try and load those:" + "Now, if parameter or observation ensemble files are available in the folder, `pyemu` will try and load those:" ] }, { @@ -1100,7 +1065,7 @@ "id": "c24eb449", "metadata": {}, "source": [ - "This provides a practical way of checking that all the plumibg works...that there arent silly mistakes with parameter values and bounds etc...and that you are happy with how the prior knowledge is being expressed through parameterisation. It can be quite helpfull to show results at this stage to stakeholders for example, to ensure that everyone agrees on the \"reasonableness\" of parameter values and distirbutions." + "This provides a practical way of checking that all the plumbing works...that there aren't silly mistakes with parameter values and bounds etc...and that you are happy with how the prior knowledge is being expressed through parameterisation. It can be quite helpful to show results at this stage to stakeholders for example, to ensure that everyone agrees on the \"reasonableness\" of parameter values and distributions." ] }, { @@ -1110,9 +1075,9 @@ "source": [ "# Finishing up the PEST setup\n", "\n", - "Once we are happy with that, we can go through the process of adding in the rest of the PEST setup. Such as observcations, other paramters...and importnatly the mf6 model run.\n", + "Once we are happy with that, we can go through the process of adding in the rest of the PEST setup. Such as observations, other parameters...and importantly the mf6 model run.\n", "\n", - "To make our lives easier later on when we use DSI, we are simply going to track the model simualted heads and temperatures in all cells and all stressperiods. We have prepared a utulity function that processes `mf6` outpfiles and writes observations to clean .txt files (see `herebedragons.py` for details; its pretty simple). \n", + "To make our lives easier later on when we use DSI, we are simply going to track the model simulated heads and temperatures in all cells and all stressperiods. We have prepared a utility function that processes `mf6` output files and writes observations to clean .txt files (see `herebedragons.py` for details; it's pretty simple). \n", "\n", "Lets just call it and run here:" ] @@ -1265,9 +1230,9 @@ "id": "eaeff0cf", "metadata": {}, "source": [ - "Now we can draw the prior. Lets draw a large number of reals so we can play around with DSI later. \n", + "Now we can draw the prior. Let's draw a large number of reals so we can play around with DSI later. \n", "\n", - "And that is it...pest setup ready to run a prior monte carlo." + "And that is it...pest setup ready to run a Prior Monte Carlo." ] }, { @@ -1290,11 +1255,19 @@ "print(pe.shape,pf.pst.npar,pf.pst.npar_adj)\n", "assert pe.shape[1] == pf.pst.npar_adj" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "9db9494d-d251-4f09-a108-1a2ab4eca552", + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { "kernelspec": { - "display_name": "gwe", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, diff --git a/02_priomc.ipynb b/02_priomc.ipynb index 3fde6be..0ba699c 100644 --- a/02_priomc.ipynb +++ b/02_priomc.ipynb @@ -24,10 +24,10 @@ "source": [ "# Intro\n", "\n", - "This notebook runs a prior monte carlo using the pest setup generated in the `setup_pstfrom` notebook. We will run the prior parameter ensemble once and stop. We will take a look at some of the results.\n", + "This notebook runs a Prior Monte Carlo using the pest setup generated in the `setup_pstfrom` notebook. We will run the prior parameter ensemble once and stop. We will take a look at some of the results.\n", "\n", "### Warning\n", - "This notebook can take some time to run. That is the cost of a >5min model run time...imagine if your model takes longer than that (#suffering). For context, this takes about 60min on a MacBook Pro. Expect longer on Windows. We recommend setting this to run over night or when you go out for lunch, before proceeding to the next notebooks. " + "This notebook can take some time to run. That is the cost of a >5 min model run time...imagine if your model takes longer than that (#suffering). For context, this takes about 600 min on a MacBook Pro. Expect longer on Windows. We recommend setting this to run over night or when you go out for lunch, before proceeding to the next notebooks. " ] }, { @@ -117,7 +117,7 @@ "id": "98fd3310", "metadata": {}, "source": [ - "### Warning: set num_workers according to your reserouces\n", + "### Warning: set num_workers according to your resources\n", "\n", "The model took ~1 min to run... so that would be about 16.7h for the full 100 runs in serial. Gonna need to run them in parallel to be reasonable. \n", "\n", @@ -144,11 +144,11 @@ "source": [ "\n", "pyemu.os_utils.start_workers(t_d, # the folder which contains the \"template\" PEST dataset\n", - " 'pestpp-ies', #the PEST software version we want to run\n", - " 'pest.pst', # the control file to use with PEST\n", - " num_workers=num_workers, #how many agents to deploy\n", - " worker_root='.', #where to deploy the agent directories; relative to where python is running\n", - " master_dir=m_d, #the manager directory\n", + " 'pestpp-ies', #the PEST software version we want to run\n", + " 'pest.pst', # the control file to use with PEST\n", + " num_workers=num_workers, #how many agents to deploy\n", + " worker_root='.', #where to deploy the agent directories; relative to where python is running\n", + " master_dir=m_d, #the manager directory\n", " )" ] }, @@ -210,7 +210,7 @@ "source": [ "## Make some figures\n", "\n", - "Lets look at some of the simualted max temperature fields." + "Lets look at some of the simulated max temperature fields." ] }, { @@ -255,11 +255,19 @@ " plt.show()\n", " plt.close();" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "30ca6c2a-225c-4407-809a-9f246f95f10f", + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { "kernelspec": { - "display_name": "gwe", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, diff --git a/03_dsi.ipynb b/03_dsi.ipynb index c30c216..4c0b2ed 100755 --- a/03_dsi.ipynb +++ b/03_dsi.ipynb @@ -25,17 +25,17 @@ "source": [ "# Intro\n", "\n", - "We saw that this awfully long running model (~5min) is going to be a major pain to history match because it takes too long and we dont have a huge (and free) cluster at our disposal...so what can we do? Surrogate models to the rescue!\n", + "We saw that this awfully long running model (~5min) is going to be a major pain to history match because it takes too long and we don't have a huge (and free) cluster at our disposal...so what can we do? Surrogate models to the rescue!\n", "\n", "We wont go into the details on data space inversion (DSI). See the various GMDSI videos and tutorial notebooks for more technical details and background. Here we focus on showing how to implement it with `pyemu`.\n", "\n", "What we need:\n", - " - results from a prior monte carlo\n", + " - results from a Prior Monte Carlo\n", " - some data to assimilate\n", " - some prediction we care about\n", "\n", "## New things in this tutorial\n", - " - We introduce `pestpp` new (as of Nov 2025) external run manager. This allows the user to handle paralel forward runs without relying on pest \"workers\"." + " - We introduce `pestpp` new (as of Nov 2025) external run manager. This allows the user to handle parallel forward runs without relying on pest \"workers\"." ] }, { @@ -79,7 +79,7 @@ "source": [ "## Choose the truth\n", "\n", - "Lets choose one of the realizations as the rtuth, then remove if form the training data. We are going to select an inconvenient truth: one in which the max temperature plume extends a bit farther than most." + "Lets choose one of the realizations as the truth, then remove if form the training data. We are going to select an inconvenient truth: one in which the max temperature plume extends a bit farther than most." ] }, { @@ -136,7 +136,7 @@ "id": "31d09618", "metadata": {}, "source": [ - "Drop the truth form the training data set" + "Drop the truth from the training data set" ] }, { @@ -328,7 +328,7 @@ "source": [ "# Start DSI\n", "\n", - "For DSI, we dont need all the observations we have been tracking. We only really need the non-zero weighted obs and the predictions of interest. " + "For DSI, we don't need all the observations we have been tracking. We only really need the non-zero weighted obs and the predictions of interest. " ] }, { @@ -350,7 +350,7 @@ "id": "534a3b6e", "metadata": {}, "source": [ - "Lets use normal score transformation (usualy a good idea...)" + "Lets use normal score transformation (usually a good idea...)" ] }, { @@ -396,9 +396,9 @@ "source": [ "Nothing new so far. But now we can introduce some new functionality in `pyemu.DSI`.\n", "\n", - "Recall that the dsi \"model\" is really just a matrix multiplication. We are passing in the `pval` array of parameter values and mulplitying it by the projection matrix. This is easly vectorizable...\n", + "Recall that the dsi \"model\" is really just a matrix multiplication. We are passing in the `pval` array of parameter values and multiplying it by the projection matrix. This is easily vectorizable...\n", "\n", - "Instread of passing each realization of `pvals` once at a time, we can do the same matrix multiplication on an ensemble of `pvals` - all in one go. Traditionaly, `pest` and `pestpp` were not designed for this. However, `pestpp` now allows a user to handle the \"paralelisation\" externaly. \n", + "Instead of passing each realization of `pvals` once at a time, we can do the same matrix multiplication on an ensemble of `pvals` - all in one go. Traditionally, `pest` and `pestpp` were not designed for this. However, `pestpp` now allows a user to handle the \"parallelisation\" externally. \n", "\n", "Lets quickly go through what vectorization of `DSI` looks like, then demo how to run the external manager with `pestpp-ies`." ] @@ -408,7 +408,7 @@ "id": "258084e4", "metadata": {}, "source": [ - "First, genera te a dummy ensemble of DSi parameter values:" + "First, generate a dummy ensemble of DSi parameter values:" ] }, { @@ -480,7 +480,7 @@ "id": "6a01ae4a", "metadata": {}, "source": [ - "key detail now, spceify the `use_runstor=True` argument to setup the forward run using the external manager:" + "key detail now, specify the `use_runstor=True` argument to setup the forward run using the external manager:" ] }, { @@ -541,7 +541,7 @@ "metadata": {}, "source": [ "And we are good to go...sheesh that was hard.\n", - "Lets just run for 2 iterations. Why? Because DSI will usualy get a good fit pretty fast. Feel free to experiment with more or less if you like.\n", + "Lets just run for 2 iterations. Why? Because DSI will usually get a good fit pretty fast. Feel free to experiment with more or less if you like.\n", "\n", "Let's go a bit higher with the number of DSi reals that we had for the training data set:" ] @@ -614,7 +614,7 @@ "id": "cae57f5b", "metadata": {}, "source": [ - "Now we will run `pestpp-ies` with the external run manager option. We do this by calling `pestpp-ies [controlfile].pst /e`. The `/e` trigegrs the external run manager option.\n", + "Now we will run `pestpp-ies` with the external run manager option. We do this by calling `pestpp-ies [controlfile].pst /e`. The `/e` triggers the external run manager option.\n", "\n", "Note that this run is in \"serial\" as far as `pestpp` is concerned. We are handling the paralelization of the dsi forward runs.\n", "\n", @@ -937,11 +937,19 @@ "ax.axvline(obs.loc[target_col].obsval, color='r')\n", "fig.tight_layout();" ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "5c51f949-4316-4a8e-af41-e0061037b6d6", + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { "kernelspec": { - "display_name": "gwe", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, diff --git a/04_dsiae.ipynb b/04_dsiae.ipynb index a77f321..f7d2d12 100644 --- a/04_dsiae.ipynb +++ b/04_dsiae.ipynb @@ -93,9 +93,9 @@ "\n", "## what we are going to do\n", "\n", - "In this noteook we are going to throught the mechanics of using the `pyemu.DSIAE` class and some of the caveats. This is a vanilla AE (think non-linear PCA), without any fancy LSTM or convolutional layers.\n", + "In this notebook we are going to through the mechanics of using the `pyemu.DSIAE` class and some of the caveats. This is a vanilla AE (think non-linear PCA), without any fancy LSTM or convolutional layers.\n", "\n", - "THe first section mirrors the same steps as in the DSI notebook.\n" + "The first section mirrors the same steps as in the DSI notebook.\n" ] }, { @@ -139,7 +139,7 @@ "source": [ "## Choose the truth\n", "\n", - "Lets choose one of the realizations as the rtuth, then remove if form the training data. We are going to select an inconvenient truth: one in which the max temperature plume extends a bit farther than most." + "Lets choose one of the realizations as the truth, then remove if form the training data. We are going to select an inconvenient truth: one in which the max temperature plume extends a bit farther than most." ] }, { @@ -388,7 +388,7 @@ "source": [ "# Start DSI-AE\n", "\n", - "For DSI, we dont need all the observations we have been tracking. We only really need the non-zero weighted obs and the prediction sof interest. " + "For DSI, we dont need all the observations we have been tracking. We only really need the non-zero weighted obs and the predictions of interest. " ] }, { @@ -410,7 +410,7 @@ "id": "534a3b6e", "metadata": {}, "source": [ - "Lets use normal score transformation (usualy a good choice..)" + "Lets use normal score transformation (usually a good choice..)" ] }, { @@ -454,12 +454,12 @@ "metadata": {}, "source": [ "## Training the AE\n", - "Here is where things start to differ. For `DSI`, \"fitting\" the model just meant doing SVD. For `DSIAE` it means training the neural net to encode/decode the data set. There is a bit more user input and options for hyperparameter tunning available.\n", + "Here is where things start to differ. For `DSI`, \"fitting\" the model just meant doing SVD. For `DSIAE` it means training the neural net to encode/decode the data set. There is a bit more user input and options for hyperparameter tuning available.\n", "\n", "### Machine learning dependency\n", "The `DSIAE` class uses `tensorflow` in the background. `Tensorflow` is an open-source machine learning framework developed by Google, widely used for building and training neural networks, including autoencoders. It provides flexible tools for both research and production, supporting CPUs, GPUs, and TPUs for scalable computation. While powerful, `tensorflow`s performance and installation can vary across operating systems—GPU acceleration, for instance, often requires careful setup of compatible CUDA and cuDNN libraries on Windows or Linux. Additionally, version compatibility between `tensorflow`, Python, and hardware drivers can sometimes pose challenges in maintaining stable environments.\n", "\n", - "It is not our intention here to provide traning on machine learning (there are many btter resources out there to lear from...). We assume the reader has some knowledge on machine learning, and the use of tensorflow. Here are some basic insghts:\n", + "It is not our intention here to provide training on machine learning (there are many better resources out there to lear from...). We assume the reader has some knowledge on machine learning, and the use of tensorflow. Here are some basic insights:\n", "\n", "* **validation_split=0.1** – Reserves 10% of the data for validation, helping monitor overfitting during training. Typical values range from 0.1 to 0.2; larger splits may be used when data are abundant.\n", "* **hidden_dims=(128, 64)** – Defines the number and size of hidden layers in the autoencoder. More layers or larger dimensions increase model capacity but may cause overfitting; smaller architectures are preferable for simpler data.\n", @@ -493,7 +493,7 @@ "id": "30acb4d0", "metadata": {}, "source": [ - "If you wish, pyemu has a built in hyperparameetr tuner (again, just using tensorflow in the background) that undertakes a grid serach to do hyperparameter tunning. Uncomment the code below if you wish to exeriment. It can take a few minutes..." + "If you wish, pyemu has a built in hyperparametr tuner (again, just using tensorflow in the background) that undertakes a grid search to do hyperparameter tuning. Uncomment the code below if you wish to experiment. It can take a few minutes..." ] }, { @@ -521,7 +521,7 @@ "source": [ "## Encoding/Decoding\n", "\n", - "A convenient aspect of DSI-AE versus vanilla DSI, is that we can directly project/encode the physics-based model outputs directly to latent space and back again. This allows us to estimate the \"error\" we incurr from simplification across data space. \n", + "A convenient aspect of DSI-AE versus vanilla DSI, is that we can directly project/encode the physics-based model outputs directly to latent space and back again. This allows us to estimate the \"error\" we incur from simplification across data space. \n", "\n", "Here is how we do it:" ] @@ -560,7 +560,7 @@ "id": "6e9abcb6", "metadata": {}, "source": [ - "We can project these back to data space by \"predciting\" (aka decoding) with the dsiae object:" + "We can project these back to data space by \"predicting\" (aka decoding) with the dsiae object:" ] }, { @@ -616,7 +616,7 @@ "source": [ "And setup the DSI pest dir...\n", "\n", - "key detail now, spceify the `use_runstor=True` argument to setup the forward run using the external manager:" + "key detail now, specify the `use_runstor=True` argument to setup the forward run using the external manager:" ] }, { @@ -660,7 +660,7 @@ "source": [ "## AE latent space prior\n", "\n", - "In vanilla DSI, the latent space prior is straightforward: Guasiian normal distributions with mean 0 and stdv 1 for all latent space parameters.\n", + "In vanilla DSI, the latent space prior is straightforward: Gaussian normal distributions with mean 0 and stdv 1 for all latent space parameters.\n", "\n", "In DSI-AE things get a bit more complicated. Autoencoder-based approaches learn a nonlinear latent space, and the resulting latent parameter distribution is generally non-Gaussian and data-driven. This flexibility allows autoencoders to capture more complex, multimodal relationships in the data, though it can make interpretation and sampling from the latent space less straightforward than in PCA-based methods.\n", "\n", @@ -686,11 +686,11 @@ "id": "357759c2", "metadata": {}, "source": [ - "If we are happy using the same number of realizations with DSIAE as we had realizations in the training data set, then we are fine. We can explicilty pass the encoded FOM prior values as the dsiae-prior parameter ensemble to pestpp-ies.\n", + "If we are happy using the same number of realizations with DSIAE as we had realizations in the training data set, then we are fine. We can explicitly pass the encoded FOM prior values as the dsiae-prior parameter ensemble to pestpp-ies.\n", "\n", - "If we want to use a larger ensemble (usualy good ideia if we can aford it..), then we need to do a bit more work to empirically sample the dsi-ae latent parameter prior distribution...\n", + "If we want to use a larger ensemble (usually good idia if we can aford it..), then we need to do a bit more work to empirically sample the dsi-ae latent parameter prior distribution...\n", "\n", - "Lets do that now, using some tricks built into pyemu. First, load the latent space prior as a `ParameterEnsemble` object. This was pre-preared when you called `.prepare_pespp()`" + "Lets do that now, using some tricks built into pyemu. First, load the latent space prior as a `ParameterEnsemble` object. This was pre-prepared when you called `.prepare_pespp()`" ] }, { @@ -760,7 +760,7 @@ "id": "4f87536c", "metadata": {}, "source": [ - "Lets comapre the distirbution of one of the parameters to see what this looks like:" + "Lets compare the distribution of one of the parameters to see what this looks like:" ] }, { @@ -808,7 +808,7 @@ "id": "633da096", "metadata": {}, "source": [ - "And write the extended parametr ensemble back to the pest dir, whislt also updateing the pespp options." + "And write the extended parameter ensemble back to the pest dir, whilst also updating the pestpp options." ] }, { @@ -830,7 +830,7 @@ "source": [ "And we are good to go...sheesh that was hard.\n", "\n", - "Lets run for 3 iterations. WThis is more that we diod for DSI. Why? Becasue the relationship is not as linear as for DSI, getting a good fit is more chalenging." + "Lets run for 3 iterations. This is more that we did for DSI. Why? Because the relationship is not as linear as for DSI, getting a good fit is more challenging." ] }, { @@ -913,9 +913,9 @@ "id": "cae57f5b", "metadata": {}, "source": [ - "Now we will run `pestpp-ies` with the external run manager option. We do this by calling `pestpp-ies [controlfile].pst /e`. The `/e` trigegrs the external run manager option.\n", + "Now we will run `pestpp-ies` with the external run manager option. We do this by calling `pestpp-ies [controlfile].pst /e`. The `/e` triggers the external run manager option.\n", "\n", - "Note that this run is in \"serial\" as far as `pestpp` is concerned. We are handling the paralelization of the dsi forward runs.\n", + "Note that this run is in \"serial\" as far as `pestpp` is concerned. We are handling the parallelization of the dsi forward runs.\n", "\n", "Here we go!" ] @@ -953,7 +953,7 @@ "id": "af837509", "metadata": {}, "source": [ - "See how it take slonger to get better fits?" + "See how it takes longer to get better fits?" ] }, { @@ -1248,7 +1248,7 @@ ], "metadata": { "kernelspec": { - "display_name": "gwe", + "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, diff --git a/model/mfsim.lst b/model/mfsim.lst deleted file mode 100644 index 6357846..0000000 --- a/model/mfsim.lst +++ /dev/null @@ -1,3088 +0,0 @@ - MODFLOW 6 - U.S. GEOLOGICAL SURVEY MODULAR HYDROLOGIC MODEL - VERSION 6.7.0.dev3 (preliminary) 10/16/2025 - ***DEVELOP MODE*** - - MODFLOW 6 compiled Oct 16 2025 14:05:53 with Intel(R) Fortran Intel(R) 64 - Compiler Classic for applications running on Intel(R) 64, Version 2021.7.0 - Build 20220726_000000 - -This software is preliminary or provisional and is subject to -revision. It is being provided to meet the need for timely best -science. The software has not received final approval by the U.S. -Geological Survey (USGS). No warranty, expressed or implied, is made -by the USGS or the U.S. Government as to the functionality of the -software and related material nor shall the fact of release -constitute any such warranty. The software is provided on the -condition that neither the USGS nor the U.S. Government shall be held -liable for any damages resulting from the authorized or unauthorized -use of the software. - - -As a work of the United States Government, this USGS product is -in the public domain within the United States. You can copy, -modify, distribute, and perform the work, even for commercial -purposes, all without asking permission. Additionally, USGS -waives copyright and related rights in the work worldwide -through CC0 1.0 Universal Public Domain Dedication -(https://creativecommons.org/publicdomain/zero/1.0/). - -The following GNU Lesser General Public License (LGPL) libraries -are used in this USGS product: - - SPARSKIT version 2.0 - ilut, luson, and qsplit - (https://www-users.cse.umn.edu/~saad/software/SPARSKIT/) - - RCM - Reverse Cuthill McKee Ordering - (https://people.math.sc.edu/Burkardt/f_src/rcm/rcm.html) - - BLAS - Basic Linear Algebra Subprograms Level 1 - (https://people.math.sc.edu/Burkardt/f_src/blas1_d/blas1_d.html) - - SPARSEKIT - Sparse Matrix Utility Package - amux, dperm, dvperm, rperm, and cperm - (https://people.sc.fsu.edu/~jburkardt/f77_src/sparsekit/sparsekit.html) - -The following BSD-3 License libraries are used in this USGS product: - - Modern Fortran DAG Library - Copyright (c) 2018, Jacob Williams - All rights reserved. - (https://github.com/jacobwilliams/daglib) - -MODFLOW 6 compiler options: -Isrc/libmf6_external.a.p -Isrc/libmf6core.a.p --Isrc -I../src -warn general -warn truncated_source -stand=f18 -O2 -fpe0 --no-heap-arrays -traceback -diag-disable:7416 -diag-disable:7025 --diag-disable:5268 -module src/libmf6core.a.p --gen-dep=src/libmf6core.a.p/Utilities_compilerversion.F90.o --gen-depformat=make -o src/libmf6core.a.p/Utilities_compilerversion.F90.o -c - -System command used to initiate simulation: -./mf6 - -MODFLOW was compiled using uniform precision. - -Real Variables - KIND: 8 - TINY (smallest non-zero value): 2.225074-308 - HUGE (largest value): 1.797693+308 - PRECISION: 15 - SIZE IN BITS: 64 - -Integer Variables - KIND: 4 - HUGE (largest value): 2147483647 - SIZE IN BITS: 32 - -Long Integer Variables - KIND: 8 - HUGE (largest value): 9223372036854775807 - SIZE IN BITS: 64 - -Logical Variables - KIND: 4 - SIZE IN BITS: 32 - - - OPENED mfsim.nam - FILE TYPE:NAM6 UNIT 1001 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - OPENED mfsim.tdis - FILE TYPE:TDIS6 UNIT 1002 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - OPENED gwf.nam - FILE TYPE:GWF6 UNIT 1003 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - OPENED ./gwf.disv - FILE TYPE:DISV6 UNIT 1004 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - - OPENED gwf.disv_vertices.txt - FILE TYPE:OPEN/CLOSE UNIT 1007 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwf.disv_cell2d.txt - FILE TYPE:OPEN/CLOSE UNIT 1008 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwf.npf - FILE TYPE:NPF6 UNIT 1009 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 17:56:37. - - OPENED ./gwf.ic - FILE TYPE:IC6 UNIT 1012 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - - OPENED ./gwf.oc - FILE TYPE:OC6 UNIT 1014 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwf.obs - FILE TYPE:OBS6 UNIT 1015 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwf.inj.wel - FILE TYPE:WEL6 UNIT 1017 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - AUXILIARY VARIABLE: TEMPERATURE - - OPENED ./gwf.ext.wel - FILE TYPE:WEL6 UNIT 1019 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - AUXILIARY VARIABLE: TEMPERATURE - - OPENED ./gwf.riv - FILE TYPE:RIV6 UNIT 1021 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - AUXILIARY VARIABLE: TEMPERATURE - - OPENED ./gwf.rch - FILE TYPE:RCH6 UNIT 1023 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - AUXILIARY VARIABLE: TEMPERATURE - - OPENED ./gwf.chd - FILE TYPE:CHD6 UNIT 1025 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - AUXILIARY VARIABLE: TEMPERATURE - - OPENED gwe.nam - FILE TYPE:GWE6 UNIT 1026 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:47. - # FMI6 gwe.fmi fmi - - OPENED ./gwe.disv - FILE TYPE:DISV6 UNIT 1027 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED gwe.disv_vertices.txt - FILE TYPE:OPEN/CLOSE UNIT 1030 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwe.disv_cell2d.txt - FILE TYPE:OPEN/CLOSE UNIT 1031 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwe.ic - FILE TYPE:IC6 UNIT 1032 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED ./gwe.est - FILE TYPE:EST6 UNIT 1034 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED ./gwe.adv - FILE TYPE:ADV6 UNIT 1038 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED ./gwe.cnd - FILE TYPE:CND6 UNIT 1039 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED ./gwe.ssm - FILE TYPE:SSM6 UNIT 1044 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:48. - - OPENED gwe.ssm_sources.txt - FILE TYPE:OPEN/CLOSE UNIT 1045 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwe.ssm_fileinput.txt - FILE TYPE:OPEN/CLOSE UNIT 1046 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwe.oc - FILE TYPE:OC6 UNIT 1047 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED ./gwe.obs - FILE TYPE:OBS6 UNIT 1048 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED mfsim.gwfgwe - FILE TYPE:GWF6-GWE6 UNIT 1049 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - READING SIMULATION OPTIONS - MAXIMUM NUMBER OF ERRORS THAT WILL BE STORED IS 1000 - END OF SIMULATION OPTIONS - - READING SIMULATION TIMING - - TDIS -- TEMPORAL DISCRETIZATION PACKAGE, - VERSION 1 : 11/13/2014 - INPUT READ FROM MEMPATH: __INPUT__/SIM/TDIS - PROCESSING TDIS OPTIONS - SIMULATION TIME UNIT IS SECONDS - END OF TDIS OPTIONS - PROCESSING TDIS DIMENSIONS - 61 STRESS PERIOD(S) IN SIMULATION - END OF TDIS DIMENSIONS - PROCESSING TDIS PERIODDATA - - - STRESS PERIOD LENGTH TIME STEPS MULTIPLIER FOR DELT - ---------------------------------------------------------------------------- - 1 1.000000 1 1.000 - 2 1.8408600E+07 14 1.000 - 3 1.3149000E+07 10 1.000 - 4 1.8408600E+07 14 1.000 - 5 1.3149000E+07 10 1.000 - 6 1.8408600E+07 14 1.000 - 7 1.3149000E+07 10 1.000 - 8 1.8408600E+07 14 1.000 - 9 1.3149000E+07 10 1.000 - 10 1.8408600E+07 14 1.000 - 11 1.3149000E+07 10 1.000 - 12 1.8408600E+07 14 1.000 - 13 1.3149000E+07 10 1.000 - 14 1.8408600E+07 14 1.000 - 15 1.3149000E+07 10 1.000 - 16 1.8408600E+07 14 1.000 - 17 1.3149000E+07 10 1.000 - 18 1.8408600E+07 14 1.000 - 19 1.3149000E+07 10 1.000 - 20 1.8408600E+07 14 1.000 - 21 1.3149000E+07 10 1.000 - 22 1.8408600E+07 14 1.000 - 23 1.3149000E+07 10 1.000 - 24 1.8408600E+07 14 1.000 - 25 1.3149000E+07 10 1.000 - 26 1.8408600E+07 14 1.000 - 27 1.3149000E+07 10 1.000 - 28 1.8408600E+07 14 1.000 - 29 1.3149000E+07 10 1.000 - 30 1.8408600E+07 14 1.000 - 31 1.3149000E+07 10 1.000 - 32 1.8408600E+07 14 1.000 - 33 1.3149000E+07 10 1.000 - 34 1.8408600E+07 14 1.000 - 35 1.3149000E+07 10 1.000 - 36 1.8408600E+07 14 1.000 - 37 1.3149000E+07 10 1.000 - 38 1.8408600E+07 14 1.000 - 39 1.3149000E+07 10 1.000 - 40 1.8408600E+07 14 1.000 - 41 1.3149000E+07 10 1.000 - 42 1.8408600E+07 14 1.000 - 43 1.3149000E+07 10 1.000 - 44 1.8408600E+07 14 1.000 - 45 1.3149000E+07 10 1.000 - 46 1.8408600E+07 14 1.000 - 47 1.3149000E+07 10 1.000 - 48 1.8408600E+07 14 1.000 - 49 1.3149000E+07 10 1.000 - 50 1.8408600E+07 14 1.000 - 51 1.3149000E+07 10 1.000 - 52 1.8408600E+07 14 1.000 - 53 1.3149000E+07 10 1.000 - 54 1.8408600E+07 14 1.000 - 55 1.3149000E+07 10 1.000 - 56 1.8408600E+07 14 1.000 - 57 1.3149000E+07 10 1.000 - 58 1.8408600E+07 14 1.000 - 59 1.3149000E+07 10 1.000 - 60 1.8408600E+07 14 1.000 - 61 1.3149000E+07 10 1.000 - END OF TDIS PERIODDATA - END OF SIMULATION TIMING - - READING SIMULATION MODELS - GWF6 model 1 will be created - GWE6 model 2 will be created - END OF SIMULATION MODELS - - READING SIMULATION EXCHANGES - GWF6-GWE6 exchange 1 will be created to connect model 1 with model 2 - END OF SIMULATION EXCHANGES - - READING SOLUTIONGROUP - - Creating solution: SLN_1 - - OPENED gwf.ims - FILE TYPE:IMS UNIT 1052 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - Creating solution: SLN_2 - - OPENED gwe.ims - FILE TYPE:IMS UNIT 1053 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - END OF SOLUTIONGROUP - -PROCESSING MODEL CONNECTIONS -END OF MODEL CONNECTIONS - - OPENED gwf.inj.wel_stress_period_data_2.txt - FILE TYPE:OPEN/CLOSE UNIT 1055 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwf.ext.wel_stress_period_data_2.txt - FILE TYPE:OPEN/CLOSE UNIT 1056 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwf.riv_stress_period_data_1.txt - FILE TYPE:OPEN/CLOSE UNIT 1057 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - OPENED gwf.chd_stress_period_data_1.txt - FILE TYPE:OPEN/CLOSE UNIT 1058 STATUS:OLD - FORMAT:FORMATTED ACCESS:SEQUENTIAL - ACTION:READ - - - IMS -- ITERATIVE MODEL SOLUTION PACKAGE, VERSION 6, 4/28/2017 - INPUT READ FROM UNIT 1052 - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - PROCESSING IMS OPTIONS - MODERATE OPTION: - DEFAULT SOLVER INPUT VALUES REFLECT MODERATELY NONLINEAR MODEL - END OF IMS OPTIONS - - PROCESSING IMS NONLINEAR - END OF IMS NONLINEAR DATA - **UNDER-RELAXATION WILL BE USED*** - - - PROCESSING LINEAR DATA - END OF LINEAR DATA - ***IMS LINEAR SOLVER WILL BE USED*** - - IMSLINEAR -- UNSTRUCTURED LINEAR SOLUTION PACKAGE, VERSION 8, 04/28/2017 - A symmetric matrix will be solved - - OUTER ITERATION CONVERGENCE CRITERION (DVCLOSE) = 0.100000E-03 - MAXIMUM NUMBER OF OUTER ITERATIONS (MXITER) = 200 - SOLVER PRINTOUT INDEX (IPRIMS) = 0 - NONLINEAR ITERATION METHOD (NONLINMETH) = 3 - LINEAR SOLUTION METHOD (LINMETH) = 1 - UNDER-RELAXATION WEIGHT REDUCTION FACTOR (THETA) = 0.900000E+00 - UNDER-RELAXATION WEIGHT INCREASE INCREMENT (KAPPA) = 0.100000E-03 - UNDER-RELAXATION PREVIOUS HISTORY FACTOR (GAMMA) = 0.000000E+00 - UNDER-RELAXATION MOMENTUM TERM (AMOMENTUM) = 0.000000E+00 - - SOLUTION BY THE CONJUGATE-GRADIENT METHOD - ------------------------------------------------------------------ - MAXIMUM OF 200 CALLS OF SOLUTION ROUTINE - MAXIMUM OF 50 INTERNAL ITERATIONS PER CALL TO SOLUTION ROUTINE - LINEAR ACCELERATION METHOD = CG - MATRIX PRECONDITIONING TYPE = MOD. INCOMPLETE LU - MATRIX SCALING APPROACH = NO SCALING - MATRIX REORDERING APPROACH = ORIGINAL ORDERING - NUMBER OF ORTHOGONALIZATIONS = 0 - HEAD CHANGE CRITERION FOR CLOSURE = 0.10000E-03 - RESIDUAL CHANGE CRITERION FOR CLOSURE = 0.10000E-01 - RESIDUAL CONVERGENCE OPTION = 0 - RESIDUAL CONVERGENCE NORM = INFINITY NORM - RELAXATION FACTOR = 0.97000E+00 - - - - - IMS -- ITERATIVE MODEL SOLUTION PACKAGE, VERSION 6, 4/28/2017 - INPUT READ FROM UNIT 1053 - # File generated by Flopy version 3.10.0.dev5 on 10/29/2025 at 13:54:46. - - PROCESSING IMS OPTIONS - MODERATE OPTION: - DEFAULT SOLVER INPUT VALUES REFLECT MODERATELY NONLINEAR MODEL - END OF IMS OPTIONS - - PROCESSING IMS NONLINEAR - END OF IMS NONLINEAR DATA - **UNDER-RELAXATION WILL BE USED*** - - - PROCESSING LINEAR DATA - END OF LINEAR DATA - ***IMS LINEAR SOLVER WILL BE USED*** - - IMSLINEAR -- UNSTRUCTURED LINEAR SOLUTION PACKAGE, VERSION 8, 04/28/2017 - An asymmetric matrix will be solved - - OUTER ITERATION CONVERGENCE CRITERION (DVCLOSE) = 0.100000E-05 - MAXIMUM NUMBER OF OUTER ITERATIONS (MXITER) = 1000 - SOLVER PRINTOUT INDEX (IPRIMS) = 0 - NONLINEAR ITERATION METHOD (NONLINMETH) = 3 - LINEAR SOLUTION METHOD (LINMETH) = 1 - UNDER-RELAXATION WEIGHT REDUCTION FACTOR (THETA) = 0.700000E+00 - UNDER-RELAXATION WEIGHT INCREASE INCREMENT (KAPPA) = 0.800000E-01 - UNDER-RELAXATION PREVIOUS HISTORY FACTOR (GAMMA) = 0.500000E-01 - UNDER-RELAXATION MOMENTUM TERM (AMOMENTUM) = 0.000000E+00 - MAXIMUM NUMBER OF BACKTRACKS (NUMTRACK) = 20 - BACKTRACKING TOLERANCE FACTOR (BTOL) = 0.200000E+01 - BACKTRACKING REDUCTION FACTOR (BREDUC) = 0.200000E+00 - BACKTRACKING RESIDUAL LIMIT (RES_LIM) = 0.500000E-03 - - SOLUTION BY THE BICONJUGATE-GRADIENT STABILIZED METHOD - ------------------------------------------------------------------ - MAXIMUM OF 1000 CALLS OF SOLUTION ROUTINE - MAXIMUM OF 100 INTERNAL ITERATIONS PER CALL TO SOLUTION ROUTINE - LINEAR ACCELERATION METHOD = BCGS - MATRIX PRECONDITIONING TYPE = INCOMPLETE LUT - MATRIX SCALING APPROACH = NO SCALING - MATRIX REORDERING APPROACH = ORIGINAL ORDERING - NUMBER OF ORTHOGONALIZATIONS = 2 - HEAD CHANGE CRITERION FOR CLOSURE = 0.10000E-05 - RESIDUAL CHANGE CRITERION FOR CLOSURE = 0.10000E-05 - RESIDUAL CONVERGENCE OPTION = 1 - RESIDUAL CONVERGENCE NORM = INFINITY NORM S - RELAXATION FACTOR = 0.00000E+00 - NUMBER OF LEVELS = 8 - DROP TOLERANCE = 0.10000E-02 - - - - Solving: Stress period: 1 Time step: 1 - -1 - STRESS PERIOD NO. 1, LENGTH = 1.000000 - ------------------------------------------ - NUMBER OF TIME STEPS = 1 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1.000000 - - Solving: Stress period: 2 Time step: 1 - -1 - STRESS PERIOD NO. 2, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 2 Time step: 2 - - - Solving: Stress period: 2 Time step: 3 - - - Solving: Stress period: 2 Time step: 4 - - - Solving: Stress period: 2 Time step: 5 - - - Solving: Stress period: 2 Time step: 6 - - - Solving: Stress period: 2 Time step: 7 - - - Solving: Stress period: 2 Time step: 8 - - - Solving: Stress period: 2 Time step: 9 - - - Solving: Stress period: 2 Time step: 10 - - - Solving: Stress period: 2 Time step: 11 - - - Solving: Stress period: 2 Time step: 12 - - - Solving: Stress period: 2 Time step: 13 - - - Solving: Stress period: 2 Time step: 14 - - - Solving: Stress period: 3 Time step: 1 - -1 - STRESS PERIOD NO. 3, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 3 Time step: 2 - - - Solving: Stress period: 3 Time step: 3 - - - Solving: Stress period: 3 Time step: 4 - - - Solving: Stress period: 3 Time step: 5 - - - Solving: Stress period: 3 Time step: 6 - - - Solving: Stress period: 3 Time step: 7 - - - Solving: Stress period: 3 Time step: 8 - - - Solving: Stress period: 3 Time step: 9 - - - Solving: Stress period: 3 Time step: 10 - - - Solving: Stress period: 4 Time step: 1 - -1 - STRESS PERIOD NO. 4, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 4 Time step: 2 - - - Solving: Stress period: 4 Time step: 3 - - - Solving: Stress period: 4 Time step: 4 - - - Solving: Stress period: 4 Time step: 5 - - - Solving: Stress period: 4 Time step: 6 - - - Solving: Stress period: 4 Time step: 7 - - - Solving: Stress period: 4 Time step: 8 - - - Solving: Stress period: 4 Time step: 9 - - - Solving: Stress period: 4 Time step: 10 - - - Solving: Stress period: 4 Time step: 11 - - - Solving: Stress period: 4 Time step: 12 - - - Solving: Stress period: 4 Time step: 13 - - - Solving: Stress period: 4 Time step: 14 - - - Solving: Stress period: 5 Time step: 1 - -1 - STRESS PERIOD NO. 5, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 5 Time step: 2 - - - Solving: Stress period: 5 Time step: 3 - - - Solving: Stress period: 5 Time step: 4 - - - Solving: Stress period: 5 Time step: 5 - - - Solving: Stress period: 5 Time step: 6 - - - Solving: Stress period: 5 Time step: 7 - - - Solving: Stress period: 5 Time step: 8 - - - Solving: Stress period: 5 Time step: 9 - - - Solving: Stress period: 5 Time step: 10 - - - Solving: Stress period: 6 Time step: 1 - -1 - STRESS PERIOD NO. 6, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 6 Time step: 2 - - - Solving: Stress period: 6 Time step: 3 - - - Solving: Stress period: 6 Time step: 4 - - - Solving: Stress period: 6 Time step: 5 - - - Solving: Stress period: 6 Time step: 6 - - - Solving: Stress period: 6 Time step: 7 - - - Solving: Stress period: 6 Time step: 8 - - - Solving: Stress period: 6 Time step: 9 - - - Solving: Stress period: 6 Time step: 10 - - - Solving: Stress period: 6 Time step: 11 - - - Solving: Stress period: 6 Time step: 12 - - - Solving: Stress period: 6 Time step: 13 - - - Solving: Stress period: 6 Time step: 14 - - - Solving: Stress period: 7 Time step: 1 - -1 - STRESS PERIOD NO. 7, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 7 Time step: 2 - - - Solving: Stress period: 7 Time step: 3 - - - Solving: Stress period: 7 Time step: 4 - - - Solving: Stress period: 7 Time step: 5 - - - Solving: Stress period: 7 Time step: 6 - - - Solving: Stress period: 7 Time step: 7 - - - Solving: Stress period: 7 Time step: 8 - - - Solving: Stress period: 7 Time step: 9 - - - Solving: Stress period: 7 Time step: 10 - - - Solving: Stress period: 8 Time step: 1 - -1 - STRESS PERIOD NO. 8, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 8 Time step: 2 - - - Solving: Stress period: 8 Time step: 3 - - - Solving: Stress period: 8 Time step: 4 - - - Solving: Stress period: 8 Time step: 5 - - - Solving: Stress period: 8 Time step: 6 - - - Solving: Stress period: 8 Time step: 7 - - - Solving: Stress period: 8 Time step: 8 - - - Solving: Stress period: 8 Time step: 9 - - - Solving: Stress period: 8 Time step: 10 - - - Solving: Stress period: 8 Time step: 11 - - - Solving: Stress period: 8 Time step: 12 - - - Solving: Stress period: 8 Time step: 13 - - - Solving: Stress period: 8 Time step: 14 - - - Solving: Stress period: 9 Time step: 1 - -1 - STRESS PERIOD NO. 9, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 9 Time step: 2 - - - Solving: Stress period: 9 Time step: 3 - - - Solving: Stress period: 9 Time step: 4 - - - Solving: Stress period: 9 Time step: 5 - - - Solving: Stress period: 9 Time step: 6 - - - Solving: Stress period: 9 Time step: 7 - - - Solving: Stress period: 9 Time step: 8 - - - Solving: Stress period: 9 Time step: 9 - - - Solving: Stress period: 9 Time step: 10 - - - Solving: Stress period: 10 Time step: 1 - -1 - STRESS PERIOD NO. 10, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 10 Time step: 2 - - - Solving: Stress period: 10 Time step: 3 - - - Solving: Stress period: 10 Time step: 4 - - - Solving: Stress period: 10 Time step: 5 - - - Solving: Stress period: 10 Time step: 6 - - - Solving: Stress period: 10 Time step: 7 - - - Solving: Stress period: 10 Time step: 8 - - - Solving: Stress period: 10 Time step: 9 - - - Solving: Stress period: 10 Time step: 10 - - - Solving: Stress period: 10 Time step: 11 - - - Solving: Stress period: 10 Time step: 12 - - - Solving: Stress period: 10 Time step: 13 - - - Solving: Stress period: 10 Time step: 14 - - - Solving: Stress period: 11 Time step: 1 - -1 - STRESS PERIOD NO. 11, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 11 Time step: 2 - - - Solving: Stress period: 11 Time step: 3 - - - Solving: Stress period: 11 Time step: 4 - - - Solving: Stress period: 11 Time step: 5 - - - Solving: Stress period: 11 Time step: 6 - - - Solving: Stress period: 11 Time step: 7 - - - Solving: Stress period: 11 Time step: 8 - - - Solving: Stress period: 11 Time step: 9 - - - Solving: Stress period: 11 Time step: 10 - - - Solving: Stress period: 12 Time step: 1 - -1 - STRESS PERIOD NO. 12, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 12 Time step: 2 - - - Solving: Stress period: 12 Time step: 3 - - - Solving: Stress period: 12 Time step: 4 - - - Solving: Stress period: 12 Time step: 5 - - - Solving: Stress period: 12 Time step: 6 - - - Solving: Stress period: 12 Time step: 7 - - - Solving: Stress period: 12 Time step: 8 - - - Solving: Stress period: 12 Time step: 9 - - - Solving: Stress period: 12 Time step: 10 - - - Solving: Stress period: 12 Time step: 11 - - - Solving: Stress period: 12 Time step: 12 - - - Solving: Stress period: 12 Time step: 13 - - - Solving: Stress period: 12 Time step: 14 - - - Solving: Stress period: 13 Time step: 1 - -1 - STRESS PERIOD NO. 13, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 13 Time step: 2 - - - Solving: Stress period: 13 Time step: 3 - - - Solving: Stress period: 13 Time step: 4 - - - Solving: Stress period: 13 Time step: 5 - - - Solving: Stress period: 13 Time step: 6 - - - Solving: Stress period: 13 Time step: 7 - - - Solving: Stress period: 13 Time step: 8 - - - Solving: Stress period: 13 Time step: 9 - - - Solving: Stress period: 13 Time step: 10 - - - Solving: Stress period: 14 Time step: 1 - -1 - STRESS PERIOD NO. 14, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 14 Time step: 2 - - - Solving: Stress period: 14 Time step: 3 - - - Solving: Stress period: 14 Time step: 4 - - - Solving: Stress period: 14 Time step: 5 - - - Solving: Stress period: 14 Time step: 6 - - - Solving: Stress period: 14 Time step: 7 - - - Solving: Stress period: 14 Time step: 8 - - - Solving: Stress period: 14 Time step: 9 - - - Solving: Stress period: 14 Time step: 10 - - - Solving: Stress period: 14 Time step: 11 - - - Solving: Stress period: 14 Time step: 12 - - - Solving: Stress period: 14 Time step: 13 - - - Solving: Stress period: 14 Time step: 14 - - - Solving: Stress period: 15 Time step: 1 - -1 - STRESS PERIOD NO. 15, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 15 Time step: 2 - - - Solving: Stress period: 15 Time step: 3 - - - Solving: Stress period: 15 Time step: 4 - - - Solving: Stress period: 15 Time step: 5 - - - Solving: Stress period: 15 Time step: 6 - - - Solving: Stress period: 15 Time step: 7 - - - Solving: Stress period: 15 Time step: 8 - - - Solving: Stress period: 15 Time step: 9 - - - Solving: Stress period: 15 Time step: 10 - - - Solving: Stress period: 16 Time step: 1 - -1 - STRESS PERIOD NO. 16, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 16 Time step: 2 - - - Solving: Stress period: 16 Time step: 3 - - - Solving: Stress period: 16 Time step: 4 - - - Solving: Stress period: 16 Time step: 5 - - - Solving: Stress period: 16 Time step: 6 - - - Solving: Stress period: 16 Time step: 7 - - - Solving: Stress period: 16 Time step: 8 - - - Solving: Stress period: 16 Time step: 9 - - - Solving: Stress period: 16 Time step: 10 - - - Solving: Stress period: 16 Time step: 11 - - - Solving: Stress period: 16 Time step: 12 - - - Solving: Stress period: 16 Time step: 13 - - - Solving: Stress period: 16 Time step: 14 - - - Solving: Stress period: 17 Time step: 1 - -1 - STRESS PERIOD NO. 17, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 17 Time step: 2 - - - Solving: Stress period: 17 Time step: 3 - - - Solving: Stress period: 17 Time step: 4 - - - Solving: Stress period: 17 Time step: 5 - - - Solving: Stress period: 17 Time step: 6 - - - Solving: Stress period: 17 Time step: 7 - - - Solving: Stress period: 17 Time step: 8 - - - Solving: Stress period: 17 Time step: 9 - - - Solving: Stress period: 17 Time step: 10 - - - Solving: Stress period: 18 Time step: 1 - -1 - STRESS PERIOD NO. 18, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 18 Time step: 2 - - - Solving: Stress period: 18 Time step: 3 - - - Solving: Stress period: 18 Time step: 4 - - - Solving: Stress period: 18 Time step: 5 - - - Solving: Stress period: 18 Time step: 6 - - - Solving: Stress period: 18 Time step: 7 - - - Solving: Stress period: 18 Time step: 8 - - - Solving: Stress period: 18 Time step: 9 - - - Solving: Stress period: 18 Time step: 10 - - - Solving: Stress period: 18 Time step: 11 - - - Solving: Stress period: 18 Time step: 12 - - - Solving: Stress period: 18 Time step: 13 - - - Solving: Stress period: 18 Time step: 14 - - - Solving: Stress period: 19 Time step: 1 - -1 - STRESS PERIOD NO. 19, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 19 Time step: 2 - - - Solving: Stress period: 19 Time step: 3 - - - Solving: Stress period: 19 Time step: 4 - - - Solving: Stress period: 19 Time step: 5 - - - Solving: Stress period: 19 Time step: 6 - - - Solving: Stress period: 19 Time step: 7 - - - Solving: Stress period: 19 Time step: 8 - - - Solving: Stress period: 19 Time step: 9 - - - Solving: Stress period: 19 Time step: 10 - - - Solving: Stress period: 20 Time step: 1 - -1 - STRESS PERIOD NO. 20, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 20 Time step: 2 - - - Solving: Stress period: 20 Time step: 3 - - - Solving: Stress period: 20 Time step: 4 - - - Solving: Stress period: 20 Time step: 5 - - - Solving: Stress period: 20 Time step: 6 - - - Solving: Stress period: 20 Time step: 7 - - - Solving: Stress period: 20 Time step: 8 - - - Solving: Stress period: 20 Time step: 9 - - - Solving: Stress period: 20 Time step: 10 - - - Solving: Stress period: 20 Time step: 11 - - - Solving: Stress period: 20 Time step: 12 - - - Solving: Stress period: 20 Time step: 13 - - - Solving: Stress period: 20 Time step: 14 - - - Solving: Stress period: 21 Time step: 1 - -1 - STRESS PERIOD NO. 21, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 21 Time step: 2 - - - Solving: Stress period: 21 Time step: 3 - - - Solving: Stress period: 21 Time step: 4 - - - Solving: Stress period: 21 Time step: 5 - - - Solving: Stress period: 21 Time step: 6 - - - Solving: Stress period: 21 Time step: 7 - - - Solving: Stress period: 21 Time step: 8 - - - Solving: Stress period: 21 Time step: 9 - - - Solving: Stress period: 21 Time step: 10 - - - Solving: Stress period: 22 Time step: 1 - -1 - STRESS PERIOD NO. 22, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 22 Time step: 2 - - - Solving: Stress period: 22 Time step: 3 - - - Solving: Stress period: 22 Time step: 4 - - - Solving: Stress period: 22 Time step: 5 - - - Solving: Stress period: 22 Time step: 6 - - - Solving: Stress period: 22 Time step: 7 - - - Solving: Stress period: 22 Time step: 8 - - - Solving: Stress period: 22 Time step: 9 - - - Solving: Stress period: 22 Time step: 10 - - - Solving: Stress period: 22 Time step: 11 - - - Solving: Stress period: 22 Time step: 12 - - - Solving: Stress period: 22 Time step: 13 - - - Solving: Stress period: 22 Time step: 14 - - - Solving: Stress period: 23 Time step: 1 - -1 - STRESS PERIOD NO. 23, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 23 Time step: 2 - - - Solving: Stress period: 23 Time step: 3 - - - Solving: Stress period: 23 Time step: 4 - - - Solving: Stress period: 23 Time step: 5 - - - Solving: Stress period: 23 Time step: 6 - - - Solving: Stress period: 23 Time step: 7 - - - Solving: Stress period: 23 Time step: 8 - - - Solving: Stress period: 23 Time step: 9 - - - Solving: Stress period: 23 Time step: 10 - - - Solving: Stress period: 24 Time step: 1 - -1 - STRESS PERIOD NO. 24, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 24 Time step: 2 - - - Solving: Stress period: 24 Time step: 3 - - - Solving: Stress period: 24 Time step: 4 - - - Solving: Stress period: 24 Time step: 5 - - - Solving: Stress period: 24 Time step: 6 - - - Solving: Stress period: 24 Time step: 7 - - - Solving: Stress period: 24 Time step: 8 - - - Solving: Stress period: 24 Time step: 9 - - - Solving: Stress period: 24 Time step: 10 - - - Solving: Stress period: 24 Time step: 11 - - - Solving: Stress period: 24 Time step: 12 - - - Solving: Stress period: 24 Time step: 13 - - - Solving: Stress period: 24 Time step: 14 - - - Solving: Stress period: 25 Time step: 1 - -1 - STRESS PERIOD NO. 25, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 25 Time step: 2 - - - Solving: Stress period: 25 Time step: 3 - - - Solving: Stress period: 25 Time step: 4 - - - Solving: Stress period: 25 Time step: 5 - - - Solving: Stress period: 25 Time step: 6 - - - Solving: Stress period: 25 Time step: 7 - - - Solving: Stress period: 25 Time step: 8 - - - Solving: Stress period: 25 Time step: 9 - - - Solving: Stress period: 25 Time step: 10 - - - Solving: Stress period: 26 Time step: 1 - -1 - STRESS PERIOD NO. 26, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 26 Time step: 2 - - - Solving: Stress period: 26 Time step: 3 - - - Solving: Stress period: 26 Time step: 4 - - - Solving: Stress period: 26 Time step: 5 - - - Solving: Stress period: 26 Time step: 6 - - - Solving: Stress period: 26 Time step: 7 - - - Solving: Stress period: 26 Time step: 8 - - - Solving: Stress period: 26 Time step: 9 - - - Solving: Stress period: 26 Time step: 10 - - - Solving: Stress period: 26 Time step: 11 - - - Solving: Stress period: 26 Time step: 12 - - - Solving: Stress period: 26 Time step: 13 - - - Solving: Stress period: 26 Time step: 14 - - - Solving: Stress period: 27 Time step: 1 - -1 - STRESS PERIOD NO. 27, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 27 Time step: 2 - - - Solving: Stress period: 27 Time step: 3 - - - Solving: Stress period: 27 Time step: 4 - - - Solving: Stress period: 27 Time step: 5 - - - Solving: Stress period: 27 Time step: 6 - - - Solving: Stress period: 27 Time step: 7 - - - Solving: Stress period: 27 Time step: 8 - - - Solving: Stress period: 27 Time step: 9 - - - Solving: Stress period: 27 Time step: 10 - - - Solving: Stress period: 28 Time step: 1 - -1 - STRESS PERIOD NO. 28, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 28 Time step: 2 - - - Solving: Stress period: 28 Time step: 3 - - - Solving: Stress period: 28 Time step: 4 - - - Solving: Stress period: 28 Time step: 5 - - - Solving: Stress period: 28 Time step: 6 - - - Solving: Stress period: 28 Time step: 7 - - - Solving: Stress period: 28 Time step: 8 - - - Solving: Stress period: 28 Time step: 9 - - - Solving: Stress period: 28 Time step: 10 - - - Solving: Stress period: 28 Time step: 11 - - - Solving: Stress period: 28 Time step: 12 - - - Solving: Stress period: 28 Time step: 13 - - - Solving: Stress period: 28 Time step: 14 - - - Solving: Stress period: 29 Time step: 1 - -1 - STRESS PERIOD NO. 29, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 29 Time step: 2 - - - Solving: Stress period: 29 Time step: 3 - - - Solving: Stress period: 29 Time step: 4 - - - Solving: Stress period: 29 Time step: 5 - - - Solving: Stress period: 29 Time step: 6 - - - Solving: Stress period: 29 Time step: 7 - - - Solving: Stress period: 29 Time step: 8 - - - Solving: Stress period: 29 Time step: 9 - - - Solving: Stress period: 29 Time step: 10 - - - Solving: Stress period: 30 Time step: 1 - -1 - STRESS PERIOD NO. 30, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 30 Time step: 2 - - - Solving: Stress period: 30 Time step: 3 - - - Solving: Stress period: 30 Time step: 4 - - - Solving: Stress period: 30 Time step: 5 - - - Solving: Stress period: 30 Time step: 6 - - - Solving: Stress period: 30 Time step: 7 - - - Solving: Stress period: 30 Time step: 8 - - - Solving: Stress period: 30 Time step: 9 - - - Solving: Stress period: 30 Time step: 10 - - - Solving: Stress period: 30 Time step: 11 - - - Solving: Stress period: 30 Time step: 12 - - - Solving: Stress period: 30 Time step: 13 - - - Solving: Stress period: 30 Time step: 14 - - - Solving: Stress period: 31 Time step: 1 - -1 - STRESS PERIOD NO. 31, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 31 Time step: 2 - - - Solving: Stress period: 31 Time step: 3 - - - Solving: Stress period: 31 Time step: 4 - - - Solving: Stress period: 31 Time step: 5 - - - Solving: Stress period: 31 Time step: 6 - - - Solving: Stress period: 31 Time step: 7 - - - Solving: Stress period: 31 Time step: 8 - - - Solving: Stress period: 31 Time step: 9 - - - Solving: Stress period: 31 Time step: 10 - - - Solving: Stress period: 32 Time step: 1 - -1 - STRESS PERIOD NO. 32, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 32 Time step: 2 - - - Solving: Stress period: 32 Time step: 3 - - - Solving: Stress period: 32 Time step: 4 - - - Solving: Stress period: 32 Time step: 5 - - - Solving: Stress period: 32 Time step: 6 - - - Solving: Stress period: 32 Time step: 7 - - - Solving: Stress period: 32 Time step: 8 - - - Solving: Stress period: 32 Time step: 9 - - - Solving: Stress period: 32 Time step: 10 - - - Solving: Stress period: 32 Time step: 11 - - - Solving: Stress period: 32 Time step: 12 - - - Solving: Stress period: 32 Time step: 13 - - - Solving: Stress period: 32 Time step: 14 - - - Solving: Stress period: 33 Time step: 1 - -1 - STRESS PERIOD NO. 33, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 33 Time step: 2 - - - Solving: Stress period: 33 Time step: 3 - - - Solving: Stress period: 33 Time step: 4 - - - Solving: Stress period: 33 Time step: 5 - - - Solving: Stress period: 33 Time step: 6 - - - Solving: Stress period: 33 Time step: 7 - - - Solving: Stress period: 33 Time step: 8 - - - Solving: Stress period: 33 Time step: 9 - - - Solving: Stress period: 33 Time step: 10 - - - Solving: Stress period: 34 Time step: 1 - -1 - STRESS PERIOD NO. 34, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 34 Time step: 2 - - - Solving: Stress period: 34 Time step: 3 - - - Solving: Stress period: 34 Time step: 4 - - - Solving: Stress period: 34 Time step: 5 - - - Solving: Stress period: 34 Time step: 6 - - - Solving: Stress period: 34 Time step: 7 - - - Solving: Stress period: 34 Time step: 8 - - - Solving: Stress period: 34 Time step: 9 - - - Solving: Stress period: 34 Time step: 10 - - - Solving: Stress period: 34 Time step: 11 - - - Solving: Stress period: 34 Time step: 12 - - - Solving: Stress period: 34 Time step: 13 - - - Solving: Stress period: 34 Time step: 14 - - - Solving: Stress period: 35 Time step: 1 - -1 - STRESS PERIOD NO. 35, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 35 Time step: 2 - - - Solving: Stress period: 35 Time step: 3 - - - Solving: Stress period: 35 Time step: 4 - - - Solving: Stress period: 35 Time step: 5 - - - Solving: Stress period: 35 Time step: 6 - - - Solving: Stress period: 35 Time step: 7 - - - Solving: Stress period: 35 Time step: 8 - - - Solving: Stress period: 35 Time step: 9 - - - Solving: Stress period: 35 Time step: 10 - - - Solving: Stress period: 36 Time step: 1 - -1 - STRESS PERIOD NO. 36, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 36 Time step: 2 - - - Solving: Stress period: 36 Time step: 3 - - - Solving: Stress period: 36 Time step: 4 - - - Solving: Stress period: 36 Time step: 5 - - - Solving: Stress period: 36 Time step: 6 - - - Solving: Stress period: 36 Time step: 7 - - - Solving: Stress period: 36 Time step: 8 - - - Solving: Stress period: 36 Time step: 9 - - - Solving: Stress period: 36 Time step: 10 - - - Solving: Stress period: 36 Time step: 11 - - - Solving: Stress period: 36 Time step: 12 - - - Solving: Stress period: 36 Time step: 13 - - - Solving: Stress period: 36 Time step: 14 - - - Solving: Stress period: 37 Time step: 1 - -1 - STRESS PERIOD NO. 37, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 37 Time step: 2 - - - Solving: Stress period: 37 Time step: 3 - - - Solving: Stress period: 37 Time step: 4 - - - Solving: Stress period: 37 Time step: 5 - - - Solving: Stress period: 37 Time step: 6 - - - Solving: Stress period: 37 Time step: 7 - - - Solving: Stress period: 37 Time step: 8 - - - Solving: Stress period: 37 Time step: 9 - - - Solving: Stress period: 37 Time step: 10 - - - Solving: Stress period: 38 Time step: 1 - -1 - STRESS PERIOD NO. 38, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 38 Time step: 2 - - - Solving: Stress period: 38 Time step: 3 - - - Solving: Stress period: 38 Time step: 4 - - - Solving: Stress period: 38 Time step: 5 - - - Solving: Stress period: 38 Time step: 6 - - - Solving: Stress period: 38 Time step: 7 - - - Solving: Stress period: 38 Time step: 8 - - - Solving: Stress period: 38 Time step: 9 - - - Solving: Stress period: 38 Time step: 10 - - - Solving: Stress period: 38 Time step: 11 - - - Solving: Stress period: 38 Time step: 12 - - - Solving: Stress period: 38 Time step: 13 - - - Solving: Stress period: 38 Time step: 14 - - - Solving: Stress period: 39 Time step: 1 - -1 - STRESS PERIOD NO. 39, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 39 Time step: 2 - - - Solving: Stress period: 39 Time step: 3 - - - Solving: Stress period: 39 Time step: 4 - - - Solving: Stress period: 39 Time step: 5 - - - Solving: Stress period: 39 Time step: 6 - - - Solving: Stress period: 39 Time step: 7 - - - Solving: Stress period: 39 Time step: 8 - - - Solving: Stress period: 39 Time step: 9 - - - Solving: Stress period: 39 Time step: 10 - - - Solving: Stress period: 40 Time step: 1 - -1 - STRESS PERIOD NO. 40, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 40 Time step: 2 - - - Solving: Stress period: 40 Time step: 3 - - - Solving: Stress period: 40 Time step: 4 - - - Solving: Stress period: 40 Time step: 5 - - - Solving: Stress period: 40 Time step: 6 - - - Solving: Stress period: 40 Time step: 7 - - - Solving: Stress period: 40 Time step: 8 - - - Solving: Stress period: 40 Time step: 9 - - - Solving: Stress period: 40 Time step: 10 - - - Solving: Stress period: 40 Time step: 11 - - - Solving: Stress period: 40 Time step: 12 - - - Solving: Stress period: 40 Time step: 13 - - - Solving: Stress period: 40 Time step: 14 - - - Solving: Stress period: 41 Time step: 1 - -1 - STRESS PERIOD NO. 41, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 41 Time step: 2 - - - Solving: Stress period: 41 Time step: 3 - - - Solving: Stress period: 41 Time step: 4 - - - Solving: Stress period: 41 Time step: 5 - - - Solving: Stress period: 41 Time step: 6 - - - Solving: Stress period: 41 Time step: 7 - - - Solving: Stress period: 41 Time step: 8 - - - Solving: Stress period: 41 Time step: 9 - - - Solving: Stress period: 41 Time step: 10 - - - Solving: Stress period: 42 Time step: 1 - -1 - STRESS PERIOD NO. 42, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 42 Time step: 2 - - - Solving: Stress period: 42 Time step: 3 - - - Solving: Stress period: 42 Time step: 4 - - - Solving: Stress period: 42 Time step: 5 - - - Solving: Stress period: 42 Time step: 6 - - - Solving: Stress period: 42 Time step: 7 - - - Solving: Stress period: 42 Time step: 8 - - - Solving: Stress period: 42 Time step: 9 - - - Solving: Stress period: 42 Time step: 10 - - - Solving: Stress period: 42 Time step: 11 - - - Solving: Stress period: 42 Time step: 12 - - - Solving: Stress period: 42 Time step: 13 - - - Solving: Stress period: 42 Time step: 14 - - - Solving: Stress period: 43 Time step: 1 - -1 - STRESS PERIOD NO. 43, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 43 Time step: 2 - - - Solving: Stress period: 43 Time step: 3 - - - Solving: Stress period: 43 Time step: 4 - - - Solving: Stress period: 43 Time step: 5 - - - Solving: Stress period: 43 Time step: 6 - - - Solving: Stress period: 43 Time step: 7 - - - Solving: Stress period: 43 Time step: 8 - - - Solving: Stress period: 43 Time step: 9 - - - Solving: Stress period: 43 Time step: 10 - - - Solving: Stress period: 44 Time step: 1 - -1 - STRESS PERIOD NO. 44, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 44 Time step: 2 - - - Solving: Stress period: 44 Time step: 3 - - - Solving: Stress period: 44 Time step: 4 - - - Solving: Stress period: 44 Time step: 5 - - - Solving: Stress period: 44 Time step: 6 - - - Solving: Stress period: 44 Time step: 7 - - - Solving: Stress period: 44 Time step: 8 - - - Solving: Stress period: 44 Time step: 9 - - - Solving: Stress period: 44 Time step: 10 - - - Solving: Stress period: 44 Time step: 11 - - - Solving: Stress period: 44 Time step: 12 - - - Solving: Stress period: 44 Time step: 13 - - - Solving: Stress period: 44 Time step: 14 - - - Solving: Stress period: 45 Time step: 1 - -1 - STRESS PERIOD NO. 45, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 45 Time step: 2 - - - Solving: Stress period: 45 Time step: 3 - - - Solving: Stress period: 45 Time step: 4 - - - Solving: Stress period: 45 Time step: 5 - - - Solving: Stress period: 45 Time step: 6 - - - Solving: Stress period: 45 Time step: 7 - - - Solving: Stress period: 45 Time step: 8 - - - Solving: Stress period: 45 Time step: 9 - - - Solving: Stress period: 45 Time step: 10 - - - Solving: Stress period: 46 Time step: 1 - -1 - STRESS PERIOD NO. 46, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 46 Time step: 2 - - - Solving: Stress period: 46 Time step: 3 - - - Solving: Stress period: 46 Time step: 4 - - - Solving: Stress period: 46 Time step: 5 - - - Solving: Stress period: 46 Time step: 6 - - - Solving: Stress period: 46 Time step: 7 - - - Solving: Stress period: 46 Time step: 8 - - - Solving: Stress period: 46 Time step: 9 - - - Solving: Stress period: 46 Time step: 10 - - - Solving: Stress period: 46 Time step: 11 - - - Solving: Stress period: 46 Time step: 12 - - - Solving: Stress period: 46 Time step: 13 - - - Solving: Stress period: 46 Time step: 14 - - - Solving: Stress period: 47 Time step: 1 - -1 - STRESS PERIOD NO. 47, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 47 Time step: 2 - - - Solving: Stress period: 47 Time step: 3 - - - Solving: Stress period: 47 Time step: 4 - - - Solving: Stress period: 47 Time step: 5 - - - Solving: Stress period: 47 Time step: 6 - - - Solving: Stress period: 47 Time step: 7 - - - Solving: Stress period: 47 Time step: 8 - - - Solving: Stress period: 47 Time step: 9 - - - Solving: Stress period: 47 Time step: 10 - - - Solving: Stress period: 48 Time step: 1 - -1 - STRESS PERIOD NO. 48, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 48 Time step: 2 - - - Solving: Stress period: 48 Time step: 3 - - - Solving: Stress period: 48 Time step: 4 - - - Solving: Stress period: 48 Time step: 5 - - - Solving: Stress period: 48 Time step: 6 - - - Solving: Stress period: 48 Time step: 7 - - - Solving: Stress period: 48 Time step: 8 - - - Solving: Stress period: 48 Time step: 9 - - - Solving: Stress period: 48 Time step: 10 - - - Solving: Stress period: 48 Time step: 11 - - - Solving: Stress period: 48 Time step: 12 - - - Solving: Stress period: 48 Time step: 13 - - - Solving: Stress period: 48 Time step: 14 - - - Solving: Stress period: 49 Time step: 1 - -1 - STRESS PERIOD NO. 49, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 49 Time step: 2 - - - Solving: Stress period: 49 Time step: 3 - - - Solving: Stress period: 49 Time step: 4 - - - Solving: Stress period: 49 Time step: 5 - - - Solving: Stress period: 49 Time step: 6 - - - Solving: Stress period: 49 Time step: 7 - - - Solving: Stress period: 49 Time step: 8 - - - Solving: Stress period: 49 Time step: 9 - - - Solving: Stress period: 49 Time step: 10 - - - Solving: Stress period: 50 Time step: 1 - -1 - STRESS PERIOD NO. 50, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 50 Time step: 2 - - - Solving: Stress period: 50 Time step: 3 - - - Solving: Stress period: 50 Time step: 4 - - - Solving: Stress period: 50 Time step: 5 - - - Solving: Stress period: 50 Time step: 6 - - - Solving: Stress period: 50 Time step: 7 - - - Solving: Stress period: 50 Time step: 8 - - - Solving: Stress period: 50 Time step: 9 - - - Solving: Stress period: 50 Time step: 10 - - - Solving: Stress period: 50 Time step: 11 - - - Solving: Stress period: 50 Time step: 12 - - - Solving: Stress period: 50 Time step: 13 - - - Solving: Stress period: 50 Time step: 14 - - - Solving: Stress period: 51 Time step: 1 - -1 - STRESS PERIOD NO. 51, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 51 Time step: 2 - - - Solving: Stress period: 51 Time step: 3 - - - Solving: Stress period: 51 Time step: 4 - - - Solving: Stress period: 51 Time step: 5 - - - Solving: Stress period: 51 Time step: 6 - - - Solving: Stress period: 51 Time step: 7 - - - Solving: Stress period: 51 Time step: 8 - - - Solving: Stress period: 51 Time step: 9 - - - Solving: Stress period: 51 Time step: 10 - - - Solving: Stress period: 52 Time step: 1 - -1 - STRESS PERIOD NO. 52, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 52 Time step: 2 - - - Solving: Stress period: 52 Time step: 3 - - - Solving: Stress period: 52 Time step: 4 - - - Solving: Stress period: 52 Time step: 5 - - - Solving: Stress period: 52 Time step: 6 - - - Solving: Stress period: 52 Time step: 7 - - - Solving: Stress period: 52 Time step: 8 - - - Solving: Stress period: 52 Time step: 9 - - - Solving: Stress period: 52 Time step: 10 - - - Solving: Stress period: 52 Time step: 11 - - - Solving: Stress period: 52 Time step: 12 - - - Solving: Stress period: 52 Time step: 13 - - - Solving: Stress period: 52 Time step: 14 - - - Solving: Stress period: 53 Time step: 1 - -1 - STRESS PERIOD NO. 53, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 53 Time step: 2 - - - Solving: Stress period: 53 Time step: 3 - - - Solving: Stress period: 53 Time step: 4 - - - Solving: Stress period: 53 Time step: 5 - - - Solving: Stress period: 53 Time step: 6 - - - Solving: Stress period: 53 Time step: 7 - - - Solving: Stress period: 53 Time step: 8 - - - Solving: Stress period: 53 Time step: 9 - - - Solving: Stress period: 53 Time step: 10 - - - Solving: Stress period: 54 Time step: 1 - -1 - STRESS PERIOD NO. 54, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 54 Time step: 2 - - - Solving: Stress period: 54 Time step: 3 - - - Solving: Stress period: 54 Time step: 4 - - - Solving: Stress period: 54 Time step: 5 - - - Solving: Stress period: 54 Time step: 6 - - - Solving: Stress period: 54 Time step: 7 - - - Solving: Stress period: 54 Time step: 8 - - - Solving: Stress period: 54 Time step: 9 - - - Solving: Stress period: 54 Time step: 10 - - - Solving: Stress period: 54 Time step: 11 - - - Solving: Stress period: 54 Time step: 12 - - - Solving: Stress period: 54 Time step: 13 - - - Solving: Stress period: 54 Time step: 14 - - - Solving: Stress period: 55 Time step: 1 - -1 - STRESS PERIOD NO. 55, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 55 Time step: 2 - - - Solving: Stress period: 55 Time step: 3 - - - Solving: Stress period: 55 Time step: 4 - - - Solving: Stress period: 55 Time step: 5 - - - Solving: Stress period: 55 Time step: 6 - - - Solving: Stress period: 55 Time step: 7 - - - Solving: Stress period: 55 Time step: 8 - - - Solving: Stress period: 55 Time step: 9 - - - Solving: Stress period: 55 Time step: 10 - - - Solving: Stress period: 56 Time step: 1 - -1 - STRESS PERIOD NO. 56, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 56 Time step: 2 - - - Solving: Stress period: 56 Time step: 3 - - - Solving: Stress period: 56 Time step: 4 - - - Solving: Stress period: 56 Time step: 5 - - - Solving: Stress period: 56 Time step: 6 - - - Solving: Stress period: 56 Time step: 7 - - - Solving: Stress period: 56 Time step: 8 - - - Solving: Stress period: 56 Time step: 9 - - - Solving: Stress period: 56 Time step: 10 - - - Solving: Stress period: 56 Time step: 11 - - - Solving: Stress period: 56 Time step: 12 - - - Solving: Stress period: 56 Time step: 13 - - - Solving: Stress period: 56 Time step: 14 - - - Solving: Stress period: 57 Time step: 1 - -1 - STRESS PERIOD NO. 57, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 57 Time step: 2 - - - Solving: Stress period: 57 Time step: 3 - - - Solving: Stress period: 57 Time step: 4 - - - Solving: Stress period: 57 Time step: 5 - - - Solving: Stress period: 57 Time step: 6 - - - Solving: Stress period: 57 Time step: 7 - - - Solving: Stress period: 57 Time step: 8 - - - Solving: Stress period: 57 Time step: 9 - - - Solving: Stress period: 57 Time step: 10 - - - Solving: Stress period: 58 Time step: 1 - -1 - STRESS PERIOD NO. 58, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 58 Time step: 2 - - - Solving: Stress period: 58 Time step: 3 - - - Solving: Stress period: 58 Time step: 4 - - - Solving: Stress period: 58 Time step: 5 - - - Solving: Stress period: 58 Time step: 6 - - - Solving: Stress period: 58 Time step: 7 - - - Solving: Stress period: 58 Time step: 8 - - - Solving: Stress period: 58 Time step: 9 - - - Solving: Stress period: 58 Time step: 10 - - - Solving: Stress period: 58 Time step: 11 - - - Solving: Stress period: 58 Time step: 12 - - - Solving: Stress period: 58 Time step: 13 - - - Solving: Stress period: 58 Time step: 14 - - - Solving: Stress period: 59 Time step: 1 - -1 - STRESS PERIOD NO. 59, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 59 Time step: 2 - - - Solving: Stress period: 59 Time step: 3 - - - Solving: Stress period: 59 Time step: 4 - - - Solving: Stress period: 59 Time step: 5 - - - Solving: Stress period: 59 Time step: 6 - - - Solving: Stress period: 59 Time step: 7 - - - Solving: Stress period: 59 Time step: 8 - - - Solving: Stress period: 59 Time step: 9 - - - Solving: Stress period: 59 Time step: 10 - - - Solving: Stress period: 60 Time step: 1 - -1 - STRESS PERIOD NO. 60, LENGTH = 0.1840860E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 14 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 60 Time step: 2 - - - Solving: Stress period: 60 Time step: 3 - - - Solving: Stress period: 60 Time step: 4 - - - Solving: Stress period: 60 Time step: 5 - - - Solving: Stress period: 60 Time step: 6 - - - Solving: Stress period: 60 Time step: 7 - - - Solving: Stress period: 60 Time step: 8 - - - Solving: Stress period: 60 Time step: 9 - - - Solving: Stress period: 60 Time step: 10 - - - Solving: Stress period: 60 Time step: 11 - - - Solving: Stress period: 60 Time step: 12 - - - Solving: Stress period: 60 Time step: 13 - - - Solving: Stress period: 60 Time step: 14 - - - Solving: Stress period: 61 Time step: 1 - -1 - STRESS PERIOD NO. 61, LENGTH = 0.1314900E+08 - ------------------------------------------ - NUMBER OF TIME STEPS = 10 - MULTIPLIER FOR DELT = 1.000 - INITIAL TIME STEP SIZE = 1314900. - - Solving: Stress period: 61 Time step: 2 - - - Solving: Stress period: 61 Time step: 3 - - - Solving: Stress period: 61 Time step: 4 - - - Solving: Stress period: 61 Time step: 5 - - - Solving: Stress period: 61 Time step: 6 - - - Solving: Stress period: 61 Time step: 7 - - - Solving: Stress period: 61 Time step: 8 - - - Solving: Stress period: 61 Time step: 9 - - - Solving: Stress period: 61 Time step: 10 - - - - Solution SLN_1 summary - ---------------------------------------------------------------------- - Total formulate time: 1.303031000000611 seconds - Total solution time: .7496689999999060 seconds - - - - Solution SLN_2 summary - ---------------------------------------------------------------------- - Total formulate time: 53.76477000000007 seconds - Total solution time: 129.3818160000000 seconds - - - MEMORY MANAGER TOTAL STORAGE BY DATA TYPE, IN MEGABYTES - ------------------------------- - ALLOCATED - DATA TYPE MEMORY - ------------------------------- - Character 1.71450000E-02 - Logical 6.80000000E-05 - Integer 13.460984 - Real 24.887384 - ------------------------------- - Total 38.365581 - Virtual 0.0000000 - ------------------------------- - - - - Run end date and time (yyyy/mm/dd hh:mm:ss): 2025/11/13 20:34:57 - Elapsed run time: 4 Minutes, 43.903 Seconds - Normal termination of simulation.