Jump to content
GMS, SMS, and WMS User Forum


  • Content Count

  • Joined

  • Last visited

Community Reputation

0 Neutral

About Michal

  • Rank
    Senior Member

Profile Information

  • Location
    Czech Republic

Recent Profile Visitors

1,505 profile views
  1. Here is the relationship from the manual:
  2. Hi, I think what you are trying to achieve could be done only aproximately, because you only have to assume what the head in the model will be. Anyway, for a given cell Qb = cb (hb - h), where Qb is the flow into the cell from the source cb is the conductance between external source and the cell hb is the head assigned to the external source and h is the model calculated head in the reference cell You problably have measured hb, so I would calculate cb this way: cb = Qb / (hb-h) But as I said, you can only assume what the value of h will be. Maybe calibrating the model against the measured flow could be a better practice.
  3. Hi, it is actually quite simple. Each stress period is divided into one or more time steps. The time steps within a stress period does not have to be the same length, but can increase exponentialy in length in order to capture system response to change in stress more accurately. So by specifying a multiplier f, each following step is f multiple of the previous step. So if you have stress period of lets say 10 days and you want to solve it in 3 time steps with 4.0 as multiplier, modflow would calculate the following step lengths: {0.47619, 1.90476, 7.61905} It is not so hard to figure the realtion as for timesteps 0 to n: step_length_i=(perlen*f^i)/Sum[f^i, {i, 0, n - 1}]
  4. There is a trick that could be used with MODFLOW. One could think of the Z coordinate as Y and model it as single layer model. But this is off the table for SEAWAT due to gravity vector orientation.
  5. Hi, I think you could use the conceptual model approach. You need to think of your cross section model as a normal layered aquifer model, that has only a single row. Therefore everything has to be treated layer-wise. You need to follow the MODFLOW logic, there is no special case for cross section models. And you obviously need to add Y coordinate to your elevation dataset.
  6. Hi, well the message desribes the issue. There is only one point in the scatter point dataset you interpolate from. I would double check that the scatter points do not have this error.
  7. Hi Jonathan, I was dealing with similar issue some time ago. There is not a satisfying way how to do it in GMS. The solution depends on the complexity of the task, but it involves grid level approach. You can isolate cells by material, select boundary cells by intersecting with boundary Arc (from Top view, Orto off) and adding a Sink/Source term. In some situtions it may be convenient to map BC form conceptual model to all layers and delete from grid where appropriate. This works if the boundaries do not overlap. Or you can modify the Package outside GMS with some script.
  8. Hi Rachel, The problem is probably that the concentration observation file is not generated after your model finishes. Therefore PEST complains that it cannot find it. The reasons for this could be many. I suggest you to study how model equivalent concentration observations are generated with VM. Alternatively, you could do it manualy with MOD2OBS program from PEST utility suit for instance. When not sure what to do, read PEST manual.
  9. Hi, I think that what you describe could be achieved by general head boundary. Without specifying head at the boundary modflow cannot calculate groundwater gradient between the cell and the external region. So I think the problem is unsolvable the way you try to aproach it. I suggest using GHB package. If you have no information about heads outside the model domain, you could for instace specify GHB head the same as some initial model result head and together with lower conductance the head in the boundary condition cells will vary based on pumping/injection in the model domain. For higher conductance values the boundary will resemble constant head boundary while for lower values it will act more as no flow BC. The key is to find the right value which could be achieved through parameter estimation process.
  10. Hello, I think that MTD3MS was released long before MODFLOW-NWT. MODFLOW-NWT does not change IBOUND value to 0 if cell goes dry and this is probably the cause of problems you are facing. I think MT3D-USGS has improved handling of such situation: MT3D-USGS also includes the capability to route a solute through dry cells that may occur in the Newton-Raphson formulation of MODFLOW (that is, MODFLOW-NWT).
  11. Hi, I think it is not directly posible through MODFLOW, but maybe you could create the second drawdown dataset in postprocessing using GMS data calculator. Just compute difference in head between the first timestep of the first simulation and the transient head dataset of the second sim. That should be the drawdown you are looking for.
  12. You are welcome. I suggest you to read the MT3DMS v5.3 supplemental manual here: http://gmsdocs.aquaveo.com/mt3dms_v5_supplemental.pdf There is a section about the new TOB package. If you need data for each cell of the boundary condition, then you would probably have to prepare the input files by hand and include each BC cell as a mass-flux object/group for which the budget data are going to be saved.
  13. Hi Lalith, In the MT3D output file, there is cumulative mass budget associated with each BC type already listed. If you need more control then there is TOB package for MT3D. https://www.xmswiki.com/wiki/GMS:Transport_Observation_Package
  14. Michal

    Flow Budget

    Hi bokuhata, I assume you have maped a general head boundary (GHB) at the coastline. When you select an arc representing this boundary you see, in the tooltip area of GMS, a flow across this boundary to/from an external reservoir. On the other hand, the Zone flow in Zonebudget is the flow between the two zones in the model. From this you see it has to be different.
  15. Michal


    Hi Aly, I think there is no straight forward way to do this in GMS. I would suggest the following procedure: 1. Create the profile as a regular 3d grid with high enough resolution to carry the material information (1 row, many columns, many layers) 2. Export your grid into a text file with the following format: Cell_id X Y Z Y coordinate should be constant for all cells. 3. Swap Y<->Z 4. Import into GIS 5. Digitize the material polygons in GIS in X-Y space 6. Select the cell centers according to the material polygon and asign material ID into a new column in the atribute table (spatial join) 7. Export the atribute table into a text file 8. Copy the material column into a new file and format the new file as GMS 3d dataset format (see GMS wiki for details) 9. Import the new material dataset to the existing 3d grid 10. Select cells in GMS according to the material dataset and specify material properties as you like
  • Create New...