WEBVTT 00:00:00.000 --> 00:00:05.400 The letters "USGS" and the words "science for a changing world appear in green then fade to black. 00:00:05.400 --> 00:00:10.333 Sometimes when you are working with MODFLOW, 00:00:10.333 --> 00:00:17.999 what you want to do is to get a time series of the flow into or out of 00:00:18.000 --> 00:00:21.233 a certain feature through time. 00:00:21.233 --> 00:00:29.633 I would like to show how to do that using ModelMuse, MODFLOW, GW_Chart, and ZONEBUDGET. 00:00:29.633 --> 00:00:32.566 I'll start off with ModelMuse. 00:00:32.566 --> 00:00:34.766 Here's a simple model. 00:00:34.766 --> 00:00:37.332 I have an object here that defines a river. 00:00:37.333 --> 00:00:43.866 Let me briefly show you exactly where the river is. 00:00:43.866 --> 00:00:50.966 Boundary Conditions, River, River Stage. 00:00:50.966 --> 00:01:01.899 And these cells in the top layer of the model all define a river. 00:01:01.900 --> 00:01:04.266 Let's see the grid. 00:01:04.266 --> 00:01:14.432 Let's suppose I wanted to see the flux out through these four cells over a period of time. 00:01:14.433 --> 00:01:15.266 Here is how I would do it. 00:01:15.266 --> 00:01:22.532 First go Model|MODFLOW Packages and Programs. 00:01:22.533 --> 00:01:29.099 Then Post procecessors, ZONEBUDGET. 00:01:29.100 --> 00:01:32.600 ZONEBUDGET allows you to define composite zones but we don't need one of those here. 00:01:32.600 --> 00:01:35.633 All we need is an ordinary zone. 00:01:35.633 --> 00:01:41.533 It exports a couple of other output files besides the regular listing file. 00:01:41.533 --> 00:01:45.966 Those are useful so I'll leave those checked. 00:01:45.966 --> 00:01:50.232 No composite zones so just click OK. 00:01:50.233 --> 00:01:54.066 Now I want to define a ZONEBUDGET zone. 00:01:54.066 --> 00:01:57.032 I want to create a polygon. 00:01:57.033 --> 00:02:00.599 I want to surround these cells. 00:02:00.600 --> 00:02:03.633 I'm going to get a few extras but that's OK. 00:02:03.633 --> 00:02:13.333 This is a polygon that extends from the top of the model to the bottom of the upper aquifer. 00:02:13.333 --> 00:02:22.199 When I go to "Data Sets", I find under ZoneBudget, there is a zone. 00:02:22.200 --> 00:02:29.700 I'll make that zone number 1. 00:02:29.700 --> 00:02:32.533 Now I'm ready to go. 00:02:32.533 --> 00:02:43.366 First I will run MODFLOW. 00:02:43.366 --> 00:02:49.799 We are exporting the input files for MODFLOW. 00:02:49.800 --> 00:02:58.533 This is a model with 5 stress periods and varying recharge and daily time steps: 00:02:58.533 --> 00:03:03.833 365 time steps. 00:03:03.833 --> 00:03:08.633 It's finished exporting the MODFLOW input files and is starting to run the model. 00:03:08.633 --> 00:03:12.733 We'll go back to this ModelMonitor window. 00:03:12.733 --> 00:03:30.999 It's looking at the data generated by MODFLOW: percent discrepancy vs. time. 00:03:31.000 --> 00:03:37.066 If something went wrong with the model, you might want to kill the model early. 00:03:37.066 --> 00:03:41.299 This would allow you to figure that out. 00:03:41.300 --> 00:03:43.566 It looks like it is almost done. 00:03:43.566 --> 00:03:45.799 There. It's done. 00:03:45.800 --> 00:03:48.000 We'll close this. 00:03:48.000 --> 00:03:50.766 It's going to open the listing file. 00:03:50.766 --> 00:03:56.899 I don't want to look at that this time so I'll just close that as well. 00:03:56.900 --> 00:04:00.366 I'll press the Enter key to close that. 00:04:00.366 --> 00:04:04.332 Now I'm ready to run ZONEBUDGET. 00:04:04.333 --> 00:04:14.699 File|Export| ZONEBUDGET. 00:04:14.700 --> 00:04:18.733 Click Save. 00:04:18.733 --> 00:04:26.566 Now it's running ZONEBUDGET. 00:04:26.566 --> 00:04:30.266 That's all done. 00:04:30.266 --> 00:04:35.832 Here is the listing file. I don't want to look at that either. 00:04:35.833 --> 00:04:48.433 In ModelMuse, one thing I could do, If I wanted to, would be to import model results. 00:04:48.433 --> 00:04:51.866 I can look at the budget file. 00:04:51.866 --> 00:04:59.432 I'll change this to binary flow files. 00:04:59.433 --> 00:05:03.766 It has the right one preselected. 00:05:03.766 --> 00:05:05.832 I can click Open. 00:05:05.833 --> 00:05:16.533 It's going to take a while to read a huge file so I'm going to pause briefly. 00:05:16.533 --> 00:05:18.666 Now I'm back. 00:05:18.666 --> 00:05:30.332 I could go down here and look at the river leakage. 00:05:30.333 --> 00:05:34.866 OK, here it is. 00:05:34.866 --> 00:05:38.566 You can see that there are different values for these cells. 00:05:38.566 --> 00:05:40.766 All the rest are zeros. 00:05:40.766 --> 00:05:44.032 Let's filter those out. 00:05:44.033 --> 00:05:46.166 Filter. 00:05:46.166 --> 00:05:48.966 The value to ignore is zero. 00:05:48.966 --> 00:05:53.899 Apply. 00:05:53.900 --> 00:06:04.633 I could look at those individual values here. 00:06:04.633 --> 00:06:08.899 They are showing here what the value is as I move the cursor over it. 00:06:08.900 --> 00:06:12.866 You can see the values changing in these different cells. 00:06:12.866 --> 00:06:16.966 You could just copy those values to the clipboard and add them up. 00:06:16.966 --> 00:06:22.099 That would be a feasible way of doing it but rather time-consuming 00:06:22.100 --> 00:06:23.500 and we don't need to do that. 00:06:23.500 --> 00:06:27.433 We've run ZONEBUDGET; that has all the information we really need. 00:06:27.433 --> 00:06:36.899 So let's look at how you could use GW_Chart to look at that time series. 00:06:36.900 --> 00:06:42.900 Select "Chart type|Water budgets. 00:06:42.900 --> 00:06:46.700 ZONEBUGET is preselected. 00:06:46.700 --> 00:06:54.000 Select Open Zonebudget file. 00:06:54.000 --> 00:07:02.766 I need to find the listing file generated by ZONEBUDGET. 00:07:02.766 --> 00:07:06.199 It's going to read that. 00:07:06.200 --> 00:07:18.733 Let's select the discretization file as well. 00:07:18.733 --> 00:07:25.299 Here it is. 00:07:25.300 --> 00:07:28.866 Click OK. 00:07:28.866 --> 00:07:31.532 Now it's reading the listing file and 00:07:31.533 --> 00:07:42.833 extracting all the data for all the different budget terms from it. 00:07:42.833 --> 00:07:53.299 Let's hide everything except the one we are interested in 00:07:53.300 --> 00:07:55.766 Let's hide this one too. 00:07:55.766 --> 00:07:58.132 Let's look at river leakage. 00:07:58.133 --> 00:08:02.633 That's a plot of river leakage through time. 00:08:02.633 --> 00:08:08.499 We could save this to a file. 00:08:08.500 --> 00:08:21.733 That file is going to have all our data in it: SampleModel.txt. 00:08:21.733 --> 00:08:28.566 It saved that file. Now, I'm going to open that file with a spreadsheet. 00:08:28.566 --> 00:08:30.632 I'll use Excel this time. 00:08:30.633 --> 00:08:34.266 There are other spreadsheets; you could do the same thing with those. 00:08:34.266 --> 00:08:42.166 Select File|Open. 00:08:42.166 --> 00:08:48.466 It's going to be a text file. 00:08:48.466 --> 00:08:56.566 I will need to go to the correct directory. 00:08:56.566 --> 00:09:02.399 Find my model, my output file. 00:09:02.400 --> 00:09:04.500 Here it is. 00:09:04.500 --> 00:09:06.566 Open. 00:09:06.566 --> 00:09:13.599 It is a delimited file with tabs as the delimiter. 00:09:13.600 --> 00:09:19.700 Next. Tabs. There we go. 00:09:19.700 --> 00:09:24.866 Finish. 00:09:24.866 --> 00:09:26.566 Here are all the budget terms. 00:09:26.566 --> 00:09:28.232 Here's the one we are interested in. 00:09:28.233 --> 00:09:35.933 It shows the data all nicely arranged into a column. 00:09:35.933 --> 00:09:37.266 Here is the stress period. 00:09:37.266 --> 00:09:40.366 Here is the time step and here is the elapsed time. 00:09:40.366 --> 00:09:44.932 The letters "USGS" and the words "science for a changing world appear in green.