Jump to content
GMS, SMS, and WMS User Forum

Bruce Campbell

  • Content Count

  • Joined

  • Last visited

Community Reputation

-2 Poor

About Bruce Campbell

  • Rank
    Senior Member

Recent Profile Visitors

The recent visitors block is disabled and is not being shown to other users.

  1. Bruce Campbell

    Transient spatially distributed recharge

    Hi Sylvain -- I've had good results using a 2-D grid. I) create a 2-D grid from the 3-D grid 2) do a linear (or other type) interpolation of your raster recharge data to the to 2-D grid. These will stack up for your 120 rasters under the 2-D grid 3) Go to MODFLOW -> Optional Packages -> Recharge and the use the 2D Dataset - > Array to put the 120 datasets into the model at the appropriate stress period
  2. Bruce Campbell

    Export native modflow files with Pilot Points

    Hi Michal -- I do all of my pest work outside of GMS. I set up the pilot points and other variables in GMS, save and run the files - up to the point where the *.PST file is written. I stop the execution in GMS and do all of the PEST work in DOS prompt windows. The *.PST file usually requires extensive edits but when I get in all back together, PEST works fine.
  3. Bruce Campbell

    PEST Observation Groups

    Hi -- Wanted to see if an attribute could be added to the water-level observation coverage indicating a PEST observation group. Of course. this would have to travel through the PEST setup process and end up in the *.PST file, somehow. In conversations with John Doherty, he recommends putting the water-level observations into groups - by aquifer / model layer / other, to help with the PEST process. I've tried to think of way to edit the *.PST file but haven't been successful.
  4. Bruce Campbell

    Export native modflow files with Pilot Points

    Hi Michal -- I've had this same problem. What I've done is interpolate the pilot point values to a 2-D grid (hopefully, the same as the model 3-D grid) and then import the interpolated values into the corresponding model layer input. For example, horizontal hydraulic conductivity pilot points for the models surficial layer are interpolated to the 2-D grid then in your MODFLOW package, use the "2D Dataset -> Layer" command to place the interpolated horizontal hydraulic conductivity values into layer 1 of the model. Do this for all of your pilot point sets. Remove all other PEST variables, set up the GMS files for a forward run and try the export to native text again. From what I've seen, GMS does not export files with PEST inputs.
  5. Bruce Campbell

    Computational time of the transient model

    Hi Petra -- Sounds like a good compromise to reduce the number of model cells. Groundwater modeling always seems to be a series of compromises... Here's the NWT solver settings that Rich passed along. I went back and looked at his description of the settings and he wrote that these settings would be a good starting point for large NWT models, 1.0E-1 1500.0 100 1.0e-8 2 1 1 SPECIFIED CONTINUE 0.90 0.00001 0.000 0.10 0 50 3.0 0.9 linmeth=1 (GMRES); =2 (xmd) 1 0 3 7 1 0.0 1 5.0E-3 1.0e-3 100 XMD I've had good results with these settings and was able to reduce my runtimes. The first 2 (the closure criteria) are probably the most model-specific.
  6. Bruce Campbell

    Computational time of the transient model

    Hi Petra -- I've been struggling with long run times for a couple of large (up to about 8 million active cell) models, also. I'm doing PEST parameter estimation calibration runs across my local network and see differences in run times depending on the computer capabilities. On older machines (4-5 years old) the run times are usually around 10 hours. On a brand new machine (Dell Latitude - lots of RAM and a fast processor) the run times are about 4 hours. I'm using MODFLOW-NWT and have a unconfined surficial layer and have gotten very good results. I suspect the difference in run times you see from HUF as compared to NWT is the ability of NWT to handle dry and flooded cells better than HUF. If you don't already have one, suggest you get a new fast machine with at least 32 GB of RAM. Also, I've been back and forth with Rich Niswonger, one of the authors of the NWT code on my large model solver settings. He's run the model(s) and suggested a set of NWT solver settings for my models. Would be glad to post them here, if you would like to have them.
  7. Bruce Campbell

    PEST Observation Groups

    Hi -- Wanted to see if there was a way to designate PEST observation groups directly in GMS. I have tried editing the *.pst file after it's created by GMS but with about 92,000 head observations (from multiple map coverages representing multiple model layers) it quickly becomes problematic.
  8. Bruce Campbell

    SFR Package - Arc interpolation control

    Hi Sean -- I did the SFR2 package construction all in GMS. The segments have to be in order - from upstream to downstream and the GMS MAP module programming does a good job of this. The trick I found was having a good DEM and matching flowlines to work with. Even with all the prep work we did, there we still some segments that had to fixed by hand.
  9. Bruce Campbell

    SFR Package - Arc interpolation control

    I'm running with about 8,000 segments in the larger of the 2 model I'm working on. I don't think the type of topography is that important, just get the closed depressions out and make the flowlines match the DEM
  10. Bruce Campbell

    SFR Package - Arc interpolation control

    Hi Sean -- What we did recently to get the SFR2 package running was mainly a GIS exercise. We simplified and put all of the streams in a "blanket" type surficial layer that covers the entire active model area. Then we put a Digital Elevation Model (DEM) together of the land surface altitudes and and made it match the model grid - grid-wise. Next, use GIS magic to get the sinks (closed topographic) areas out of the DEM. Use more GIS magic to create flowlines from the DEM. We didn't have any luck with the National Hydrography Dataset - the flowlines don't match the topography and you get some arcs with the downstream ends higher than the upstream ends. When all of the GIS work was done (it's fairly high-level work) we had a DEM and flowlines that world match and work in GMS and were able to get functioning SFR2 packages. There are several advantages to the SFR2 package over the Drain or River packages - the primary one being the removal of the "infinite" term in the equation that lets the Drain and River packages take or add as much water as the head differences and conductivity allows. The SFR2 budget is limited to the groundwater that has been discharged to the stream and any surface water volumes that are directly applied.
  11. Bruce Campbell


    FYI - A new version of MODFLOW was posted on the USGS page yesterday: https://water.usgs.gov/ogw/modflow/6beta.html
  12. Bruce Campbell

    PEST Variables - Drain Conductance

    Hi -- I'm using PEST to estimate MODFLOW Drain Package conductance terms. I have a number of stream gages where I have analysed the data and calculated the groundwater baseflow component and I'm using these data as observations in PEST. In several of the basins, I have 2 (or more) gages and have a question about the PEST variable. In the stream arcs, I'm setting up a negative number in the attribute table for the conductance term (-201, -202, etc). What is the best procedure to use for having the 2 gages in the same basin. What I've done in the past is have 2 separate variables for the 2 gages. However, it seems like that cuts the simulated flow term at the lower gage. Does anyone have a better procedure?
  13. By any chance is there a Linux-compiled version of the MODFLOW-NWT_h5_64 executable available? I wanted to use a cloud environment to calibrate a GMS created model but the cloud operating system is Linux.
  14. Bruce Campbell

    Coordinate Significant Figures

    Hi Michael -- Thanks for the quick response. How about the name of a well or observation point - I have some long names from our database. As much as 30-50 characters with punctuation such as parentheses, pound signs, and periods. Could these long names be a problem? I can replace them with simple consecutive numbers.
  15. Bruce Campbell

    Coordinate Significant Figures

    Hi -- How many significant figures does GMS use when assigning coordinate locations for groundwater level observations? I've got a little over 400,000 groundwater level observations in a GMS 10.1.4 model and having trouble using PEST outside of GMS. I think what's happening is GMS is creating duplicate locations for the wells. I have 8 significant figures on the X and Y coordinates and there are no duplicates. However, if GMS is only using the first 2 or 3 of these significant figures in the coordinates, I'll have lots of duplicates.