<div dir="ltr"><div><div>Looks like changing the projection system of the data seems to be working.</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 25 June 2015 at 13:33, vishal tiwari <span dir="ltr"><<a href="mailto:hi.vishal123@gmail.com" target="_blank">hi.vishal123@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear Ma'am, <br><br></div>Thank you for your response. Over the past couple of days I tried the r.sim.water over various number of data-sets and I was able to get the results for elevation data which have resolution of size greater than 10m.<br>Although the results for very high resolution data-set were not getting generated even though the dimensions of the grid was not too big.<br><br></div>So I created two kind of data-sets each having the following meta data:<br><br>First (Not working) =><br><br> | Type of Map: raster Number of Categories: 255 |<br> | Data Type: FCELL |<br> | Rows: 49 |<br> | Columns: 131 |<br> | Total Cells: 6419 |<br> | Projection: UTM (zone 13) |<br> | N: 52.33036335 S: 52.32309305 Res: 0.00014837 |<br> | E: 13.04243697 W: 13.03048443 Res: 0.00009124 |<br> | Range of data: min = 43 max = 103<br><br></div>Second (working) =><br><br> | Type of Map: raster Number of Categories: 255 |<br> | Data Type: FCELL |<br> | Rows: 53 |<br> | Columns: 140 |<br> | Total Cells: 7420 |<br> | Projection: UTM (zone 13) |<br> | N: 52.33058068 S: -477.66941932 Res: 10 |<br> | E: 1413.03011844 W: 13.03011844 Res: 10 |<br> | Range of data: min = 43 max = 103 |<br> | |<br><div><br></div><div>First (original) and second both are the same data-set. Second is generated by changing the resolution size to 10 (in degress, thus changing the S and E coordinates, which are wrong). But the original data-set have the pixel size of 10m X 10m (but shown in degrees).<br><br>My question being does the format of the datasets projection have to be in a specific projection format for r.sim.water to work?<br><br></div><div>Also can you please help me figure out how are the "Time step" value and "Number of iterations" in the simulation are being calculated, by pointing me to some references?<br></div><div><br></div><div>Also apologies for the late relay.<br></div><div><br></div><div>Thank you and kind regards,<br></div><div>Vishal Tiwari<br></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On 13 June 2015 at 16:38, Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><div>On Wed, Jun 10, 2015 at 11:21 AM, vishal tiwari <span dir="ltr"><<a href="mailto:hi.vishal123@gmail.com" target="_blank">hi.vishal123@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi all,<br><br>I am trying to run the <code>r.sim.water</code> in GRASS. I have calculated the partial derivatives dx and dy rater map using <code>r.slope.aspect</code>.
I am using this elevation map [1] to get the time-series depth map. But
I think the output file is not being saved in the "outputDepth" rater map
(says failed to open, on seeing the metadata).<br><br>> r.sim.water -t elevin=testElevationLive@sampleMapSet dxin=dx@sampleMapSet dyin=dy@sampleMapSet rain_val=150 manin_val=0.015 depth=outputDepth nwalk=50000 outiter=20<br><br></div><div>STDOUT:<br><pre><code>Min elevation = 43.00 m
Max elevation = 103.00 m
Mean Source Rate (rainf. excess or sediment) = 0.000030 m/s or kg/m2s
Mean flow velocity = 21.486384 m/s
Mean Mannings = 0.015121
Number of iterations = 647482826 cells
Time step = 0.00 s
</code></pre>
<p>and on the command console:</p>
<pre><code>Command 'd.rast map=outputDepth@sampleMapSet' failed
Details: Raster map <outputDepth@sampleMapSet> not found
<br></code></pre>[1] <a href="https://www.dropbox.com/s/ooxcj0y6v2kzeml/testElevationLive.tif?dl=0" rel="nofollow" target="_blank">https://www.dropbox.com/s/ooxcj0y6v2kzeml/testElevationLive.tif?dl=0</a><br><br></div><div>Thanks and Regards,</div></div></blockquote><div><br><br></div></div></div><div>For the list archive record, the offlist answer by Helena Mitasova:<span><br><br>On Fri, Jun 12, 2015 at 10:52 PM, Helena Mitasova wrote:<br>> It is
possible that there is no output - with time step smaller than 0.001 it
looks like you are running it at very high resolution<br>> on a huge DEM for some extreme storm event (150mm/hour - is that what you are trying to model?)<br>> or you just don’t have your region set properly. To test, try to increase the grid cell size and use 50mm/hr storm.<br>><br>> Helena<br>><br></span></div><span><font color="#888888"><br><br></font></span></div><span><font color="#888888"><div class="gmail_quote">Markus<br></div></font></span></div></div>
</blockquote></div><br><br clear="all"><br></div></div><span class="">-- <br><div><div dir="ltr"><div><div dir="ltr">Vishal Tiwari<div>Computer Science and Engineering</div><div>Lab for Spatial Informatics</div><div>IIIT-Hyderabad</div></div></div></div></div>
</span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">Vishal Tiwari<div>Computer Science and Engineering</div><div>Lab for Spatial Informatics</div><div>IIIT-Hyderabad</div></div></div></div></div>
</div>