Jump to content
GMS, SMS, and WMS User Forum

All Activity

This stream auto-updates     

  1. Last week
  2. Chris Smemoe

    GSSHA - Long-term with GAR infiltration & Penman ET

    Joel, I answered a question in another post. The answer to this question is that it's hard to answer without seeing your model. You can make sure you have good starting values for all the evapotranspiration parameters in your model. I'd recommend the GSSHA course or contacting tech support with your model as I described in your other post. Let me know if you have any other questions. Chris
  3. Joel, I don't have an immediate answer to your question off the top of my head. I'd recommend contacting tech support to see if they can help you with your specific questions. Their contact information is here: https://www.aquaveo.com/technical-support If you could provide a copy of your model to tech support and explain the specific problem you're having, that would be best. With specific engineering-related questions, consulting is available at an additional charge. There are 2 good resources for answering GSSHA-specific questions. The first is the GSSHA wiki, located at https://gsshawiki.com Also notice on the GSSHA wiki web site that a training course is being offered next month in Alabama. This 2-day course is offered free of charge if there are still seats available. You would need to provide your own time, lodging, and transportation to the course. I think this is a great opportunity to discuss any problems you have with the GSSHA developers and to get to know them. Contact information for the course is on gsshawiki.com. This would be a good use of your time and you would get to know other people using GSSHA for their projects. Let me know if you have any other questions, Chris
  4. Michal

    RT3D

    You are welcome, I am happy it worked
  5. Kanchana

    RT3D

    Many Thanks Michal, It worked fine with smaller time steps and after I increase the grid resolution.
  6. I'm having trouble understanding the results of my long-term simulation. When I look in the summary file, there is a significant volume of "Total downward flux from unsaturated zone" (which is equal to "volume of groundwater recharge"). However, the infiltrated depth according to the grid output file says that the maximum depth of infiltration is only 0.4 meters (which is less than my Soil Moisture Depth of 1 meter). I thought that the Soil Moisture Depth parameter for long-term Green & Amt with Redistribution defined the depth of the unsaturated zone. But that doesn't jibe with the aforementioned results. If maximum infiltrated depth is 0.4 meters, but the unsaturated zone is 1 meter deep, shouldn't the "Total downward flux from unsaturated zone" be zero? Can someone help clear this up? Have I misunderstood what all of these parameters are referring to? Thanks! -Joel
  7. Earlier
  8. Michal

    RT3D

    Hi Kanchana, its realy hard to tell what may be causing the troubles. I guess it will be something in the reaction module related to the reaction rate. Try to run the simulation for lets say a minute with time steps in the order of seconds or fractions of a second. Also increasing grid resolution may improve the solution.
  9. Kanchana

    RT3D

    Hi Michal, Yes it runs OK only as conservative transport.
  10. cansudemir

    ERROR: transient modflow forward run

    Hi, I have been using Modflow-usg and made successful runs. However, there are still many bugs in USG since it's still in development. Therefore I wanted to switch to Modflow with the same data set. I have 20336 cells with cell size of 150m in x and y direction, HK and starting head values were assigned to the model from 2-D scatter data- Interpolate to Modflow. When all the conceptual model data is converted to modflow, all the things look fine (as they were before in the modflow-usg, everything is same). Simulation check is not giving any errors. However, when I run the model I continuously receive the same error and MODFLOW terminates with error without doing nothing. I have tried all the solvers. I have decreased the number of cells. Instead of starting heads I also tried checking steady state in stress period dialog. The error is: ERROR READING ARRAY. Error opening data set: Pilot Points/Index1/start vals. Error opening data set: Pilot Points/Index1/weights. MODFLOW terminated with error! Error reading ......hed. File formatted incorrectly. Other output encountered: ---Premature end of file file_name.ccf Could you please help me solving this problem. Thanks! Cansu
  11. Michal

    RT3D

    I guess it can. For instance if the reaction module causes abrupt changes in concentration of some specie. Does it run ok with MT3DMS? I mean without the reactions only as conservative transport?
  12. Kanchana

    RT3D

    Hi Michal, Many Thanks for the comment. I tried increasing Max.Trans.steps as you mentioned, but still the same error pops up. Can this be a problem related to the edited source code in user defined modules?
  13. I'm attempting to model a small basin (~0.2 sq mi) using a long-term simulation. I can get the model to run without errors or warnings, but some of the model results are counter-intuitive. So I'd like to verify if I've set up my model correctly and throw out a couple concerns/questions that some of you might be able to help me out with. Some quick background about my model setup: I'm using a single rain gauge located very near the basin boundary, defined using a .gag file with hourly data for the duration of November 2009. I've compiled WES format hydrometeorological data file that also gives hourly data for November 2009. My HMET file includes no radiation data, so I'm relying on the GSSHA calculated radiation from lat/lon/time&date. Active modules are as follows: Overland flow computed by ADE with Interception; Evapotranspiration via Penman method; Infiltration via G & A with Redistribution (Soil moisture depth = 1.0 m, Top layer depth = 1.0 m); Channel routing via Diffusive wave. My primary concern with the results is the lack of evapotranspiration calculated by the simulation. Out of ~100,000 m^3 of precipitation during the simulation period, only ~100 m^3 are evaporated. (Roughly 20% of the precip is intercepted and the other 79.9% infiltrates.) 0.1% evapotranspiration seems alarmingly low, even for a very wet month with consistent cloud cover in the Pacific Northwest. I've tried manipulating all of the parameters that feed into the Penman ET calcs (including cloud cover and wind speeds in my HMET file), and I can only increase ET by about 50 m^3. Has anyone else had such a problem? Or can anyone think of a reason that ET would be so seemingly underrepresented? Secondarily, I'm not sure if I've understood the GAR layer depth parameters. I've read the Downer paper which explains the GAR (http://www.dtic.mil/get-tr-doc/pdf?AD=ADA473824) but still am unsure about the physical meaning of "Soil moisture depth" and "Top layer depth". Can someone clear this up? Ultimately, my basin of interest has approx. 1 m of silty loam resting atop an impermeable siltstone bedrock. A perched water table develops above the siltstone during the rainy season and it would be nice to recreate that with the model, but until now I've been unsuccessful with any of the Infiltration options. Could I model such a configuration with Richards Eq. or multi-layer G&A? Thanks in advance to anyone who reads! Hope it's all clear
  14. Michal

    RT3D

    Hi Kanchana, on the stress period dialog, try to increase the Max. trans. steps variable. It is the maximum number of transport steps allowed during each time step of a flow solution. Give it a high value like 1e6 and see what happens. While performing transport simulation with MT3D based model, time steps are further divided into transport steps. When you leave 0.0 as transport step size, the model calculates transport steps size itself according to certain criteria based on the flow solution. These stability criteria relates among others to the Courant number, you can find more about these in the MT3DMS manual (Zheng, Wang, 1999).
  15. Kanchana

    RT3D

    Hi, I am trying to develop a RT3D user defined module. Following the instructions given in the tutorials, I developed a dll file and, batchrxn ran successfully. After that, when running the new RT3D simulation, it gives an error message "Number of Transport steps exceeds specified maximum (MXSTRN) = 1000", and the program terminated. But I defined only 1 stress period (75 days long) with only 75 time steps (Images attached). Could anyone please help me with this? Many Thanks Subhani
  16. Chris Smemoe

    Format of GSSHA.dep file

    David, The .dep file just has water depth solution values for each of the cells for each output time step. I don't know that the file format is documented anywhere, but maybe I can try to track something down if you need me to. The important values you would use if you were reading this file are the number of cells in the model, which is in the "ND" card at the top of the file, the name of the dataset (NAME card), whether each cell is active/inactive (these values do not change for each time step), and the solution values for each cell (these values DO change for each time step). It's just a text file, so you can open the file in a text editor and look at it to see what I'm talking about. The active/inactive and the solution values are for each cell (based on the cell ID) and you cannot determine the actual location with just this file. The values are listed in the .dep file in the order of the cell ID's. If you have the location of the center of the upper left cell of your grid and you know the number of rows and columns in your grid, you can determine the location of the center of any other cell from the cell ID. The cells are numbered starting from the upper left corner of the grid. Let me know if you have any other questions. Chris
  17. hi everybody, How do link GMS and Matlab software? can you give me any instruction for this, please? This case is very much emergency for me, It will be very helpful if you can share your experiences or suggestions. It will be very much appreciated. All the best,
  18. hi everybody, How do link GMS and Matlab software? can you give me any instruction for this, please? This case is very much emergency for me, It will be very helpful if you can share your experiences or suggestions. It will be very much appreciated.All the best,
  19. disgreen

    Format of GSSHA.dep file

    Hello, I'm assuming that *.dep file the is the data file of simulated water depths for a given GSSHA simulation. I'd like to be able to import this file into Python or R to automate evaluation of simulated water depths at various locations within a model domain. Also, what is the general structure of the .dep file, and how is this linked to the model domain? Is the linkage through cell index or through the coordinates of the cell centroids? Thanks, David
  20. glarsen

    SMS 12.3 has been Released!

    SMS 12.3 has been Updated! The 64-bit full installation can be downloaded here. The Tutorials can be found here. Bug fixes can be found here.
  21. cansudemir

    Observation Wells in Transient Model

    Hello, I am having a serious problem with the observations in my transient model. The observation wells have measurement data from 2009, 2010, 2011, 2012 and 2013. Each well has different measurement dates. Some wells have measurement of all years but some of them have less. When I run the model and plotted computed vs observed graph, in the view values part I see observations of different months and years. For example, in the plot of 16.03.2013 there are some observation values of 2009 (a well whose last observation is 2009). What is the problem? How can I erase these wells from the times that they don't belong to? Need help ASAP. Thank you!! Cansu
  22. disgreen

    Converting GIS Point Features Into Gaging Stations

    Hi Chris, Sorry for the repeated posts. I wrote a small Python script to transform the point feature data into the gage file format you gave. It worked quite well. Thank you.
  23. disgreen

    Converting GIS Point Features Into Gaging Stations

    Hi Chris, An additional question: the gage file format looks to be space delimited. Is this correct? Thanks again, David
  24. disgreen

    Converting GIS Point Features Into Gaging Stations

    Hi Chris, As always, thank you for looking into this. I'll try your suggestion and let you know if it works. Perhaps this would be a good feature to add to WMS. Also, is there a way in WMS to do a bulk extraction of time-series from multiple gages? For example, one (ideal) export structure would be gages organized by ID (as columns) with a common time-stamp. Thanks again, David
  25. mina.88

    RCH output

    okay Michael, Thanks for your quick respond but this ASCII.dat file is just one time step and i want all of them together ( 168 Time steps ). Also i need cell properties along with the values. sincerely Mina
  26. Chris Smemoe

    Converting GIS Point Features Into Gaging Stations

    David, I looked at the documentation for gages in WMS and I'm afraid it's not very helpful. There are 2 types of gages in WMS: Rain gages and monitoring gages. You can convert a shapefile to rain gages in a "Rain Gage" coverage but there's no way to directly convert a shapefile to monitoring gages. However, there is a monitoring gage file format that you can import into WMS in the Gages dialog. Here's a sample file: GAGE BEGGAG NAM Gage1 LOC 1.000000 1.000000 0.000000 DIR 0.000000 1.000000 0.000000 COL 0 0 255 BEGGAG NAM Gage2 LOC 10.000000 10.000000 0.000000 DIR 0.000000 1.000000 0.000000 COL 0 0 255 BEGGAG NAM Gage3 LOC 20.000000 20.000000 0.000000 DIR 0.000000 1.000000 0.000000 COL 0 0 255 The "BEGGAG" line tells WMS to create a new gage. The "NAM" line is the name of the gage and the "LOC" line is the X Y Z location of the gage. You could just export the XYZ points from a shapefile to a text file and modify the file to include the necessary lines to import as a gage file. The only thing you might want to change in this file (besides the location and name) is the "COL" card. This is the Red, Green, and Blue color for the gage. The "DIR" isn't really used in WMS. To import this file, just go to the 2D grid module and select Data | Gages from the menus. It's also accessible from the "Select Monitoring Gage" tool. Select import from the Gages dialog and read your file. You could also just create a bunch of new gages in the Gages dialog and copy/paste your XYZ values into the dialog spreadsheet. Let me know if I can help with anything else. Chris
  1. Load more activity
×