Jump to content
GMS, SMS, and WMS User Forum


  • Posts

  • Joined

  • Last visited

Recent Profile Visitors

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

yacob's Achievements


Newbie (1/14)



  1. Thanks for the explanation. I knew that sometimes, after a certain about of iterations, clean dam could not clear all depressions and I would indeed have to manually adjust. However for this case the Cleandam program was reporting that all dams were cleared and I have no black dots on my grid. Only after running GSSHA did I notice erroneous ponding and after viewing blue flow direction arrows did I discover that some depressions were still there. I just was wondering if that was a bug.
  2. I am having a problem with my WMS 10.1.17 where when I run the CleanDam tool and accept changes, but the resultant grid still has depressions. It isn't readily apparent until I run a GSSHA simulation and look at the ponding. When I turn on digital dams in display properties I get no display of digital dams, however when I look at the flow directions with my depth grid, I can plainly see all four arrows pointing inward and a pond. I turn on elevations display values which confirms the flow directions honor the elevation file. I have tried running CleanDam externally from WMS as well with no luck.
  3. The current version of gssha.exe (official release 5.1) does not allow you to use spatially variable HMET data in your long term simulations.
  4. Hi Chris, Thanks for your help. I figured out how the precip data is spatially recorded in the GSSHA .gag file and wrote a program to print RADAR data for all precipitation events. The way wms reads spatial precip is as follows Suppose arcinfo ascii file below with indicies 1 to 20. WMS reads the precip file from 1 to 20, and loads only those values in the file that fall within your watershed. (The best way to visualize this is to load the ascii file and a polygon of your watershed into ArcGIS) WMS then stores them in the .gag file as RADAR YYYY MM DD HH MN 1 2 3 ... 20. i.e. if your watershed polygon only covers 9,14,3,8,13,18,7,12 then your RADAR line will read RADAR YYYY MM DD HH MN 3 7 8 9 12 13 14 18 5 10 15 20 4 9 14 19 3 8 13 18 2 7 12 17 1 6 11 16 If WMS is giving you a hard time with this, just recreate the same process with your own code. As long as it is in the right format, GSSHA will read it. Also, if you import a nexrad .asc file and save the project, it will write out your coordinates in the .gag file. So keep this part of the file and paste your precip data underneath. This is just the way I did it; it may not be applicable to everyone.
  5. Hi, I am attempting to use observed precipitation data that I have aquired from TRMM (Tropical Rainfall Measuring Mission) in my GSSHA project. The data is spatially and temporally varied and I have converted it to the arcinfo grid .asc format. I am using the WMS Nexrad Radar grid conversion tool to preprocess my precipitation into the .gag file. The data is near real time, 3 hourly global precip at 0.25 degrees (-180 to 180 lon and -60 to 60 lat). When read in a 48 hour rainfall (which composes of 16 asc files) I get weird values in my .gag file. It appears that values for a given gridded gage seem to double which does NOT match the .asc file...at all. I get output like this in my .gag file 0.0 0.0 0.749 0.0 0.749 0.0 0.0 1.499 2.997 0.0 0.749 5.994 12.738 0.0 0.0 26.045 As you can see it starts to get exponentially out of control. I have tried doing Judy's Branch Tutorial (using NEXRAD raw data) and found that the .gag file looks normal after adding radar data. The only difference between the TRMM and the NEXRAD is that TRMM is global and I am building a model that is located outside the US. Can WMS help me build a precip file using the NEXRAD radar option with data that is not NEXRAD but is still in the correct arcinfo grid format or do I have to build the .gag file myself? Thanks.
  6. I think I have found a solution. Apparently it was a parameter issue. I was digging around on the gsshawiki and under primer -> capabilities -> infiltration I found that for green ampt inf "The initial moisture content must be less than or equal to the porosity and should be greater than the residual water content." I changed this in my mapping tables and was able to get a stable run. Some more testing is in order but I think I'm on the right track. I'll post if have anymore problems with this issue. Thank you all for your help and feedback. Here is the link. Happy Modeling! http://gsshawiki.com/gssha/Capabilities:Infiltration#Green_.26_Ampt_Infiltration
  7. @ Cameron If I use a text editor on the soils.idx file I find that all gridcells inside of my watershed are mapped to the correct uscs soil classification number. All cells outside of my watershed are zeros. The frustrating thing is it all looks normal. For some reason WMS doesnt show these grid cells that are mapped zero the lie outside of my water shed (which is where my model is blowing up). I have to look at the .idx file to see them all. In the mapping tables I only set up parameters for the possitive uscs classification integer. @ Chris Unfortunatley for me, I cannot disclose my model to the WMS tech support because it is classified. I'm am currently waiting for 8.3 to be approved for update and I hope that will help. I will also look at the tutorial.
  8. On second thought, how long did you run it? You may have not given it enough time for the water to runoff. Try to increase the Total time and see if that helps.
  9. I am new to gssha but it seems to me that your .sum file is saying that all the water is either being stored in the surface/channels or has been infiltrated. Load your surface depth (.dep file) solution to your grid, contour it, and look for any depressions. With that amount of volume remaining on the surface it should be evident where all the water is pooling up. When you find it you can edit your elevations to move the water. Also did you clean digital dams before you ran the model? Hope this helps.
  10. Hi. I am attempting to run a gssha model with infiltration (green ampt with soil moisture redistribution or "GAR"), however either one of two things are happening. My model crashes if I use Ksat and initial soil moisture that are too low. Conversely if I increase these parameters to the point where the model will run, all the water infiltrates and I get no runoff. The crash error states "h NaN in main()- point 3, grid cell xxxxxx". Something interesting to note here is that the gridcell it says has crashed is actually outside of my watershed. How I set up the index map: I map shapefile to "soil type" coverage. In the "GIS to Feature objects wizard, I map USCS soil classification attribute in my shapefile to USCS soil type. Then I create an index map for gssha. The input coverage is my soils coverage and the coverage attribute is ID. I then set up the mapping tables and run the model. The model runs fine without infiltration and I can get overland flow (ADE-PC) and stream flow routing(Diffusive Wave). I seems I have tried every conceivable combinations of infiltration parameter values and soil moistures. (within typical ranges of the wiki table) If I create a single soil type for the entire watershed it still fails. I have redrawn my shapefiles in case they may have been corrupted. I have lowered the timestep down to 1 sec. The highest initial moisture values it will take is .29 which, too me, is pretty low so no wonder its all infiltrating. I'm not really sure what is causing cells outside of my grid to blow up. If it helps my gridcell resolution is 500m x 500m and I'm using WMS 8.1. Thanks in advance for the help.
  • Create New...