Thanks Soren,<br>From you response, can you please tell me how to do the following two tasks (which i dont find in the online manual for r.gwflow): the remainder of your comments i have figured out.<br><br>- How can i compute the flux in [m^3/s] for each cell with r.mapcalc?<br>
and<br><br>- i have set the bc of the edges and stream cells at constant head for now- how can i get the budget raster maps you mention?<br><br>thanks again,<br>Vishal<br><br><div class="gmail_quote">On Sat, Jul 14, 2012 at 12:05 PM, Sören Gebbert <span dir="ltr"><<a href="mailto:soerengebbert@googlemail.com" target="_blank">soerengebbert@googlemail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
sorry for the delay.<br>
<br>
2012/7/11 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>>:<br>
<div class="im">> Thanks Soren,<br>
> That explains some of the results i'm getting, with water piling up above<br>
> the surface in the edges in low-lying areas.<br>
><br>
> Can you please tell me how i can change that to constant flux or constant<br>
> head? If constant flux, should that be in m/s units?<br>
<br>
</div>Constant flux can currently only be defined using sources/sinks with<br>
unit [m^3/s], that is option q.<br>
But i can add two new options (fn, fe) that defines the flux in<br>
northern or eastern direction using the unit [m/s]<br>
that will be multiplied internally with the northern or eastern face<br>
area of the cell?<br>
<br>
Otherwise you need to compute the flux in [m^3/s] for each cell with r.mapcalc.<br>
<div class="im"><br>
><br>
> My problem though is that i dont know what a constant flux or head at the<br>
> edges should be set to. For now the only bc i have put in there deliberately<br>
> (beyond the default you mention) is that i have set constant head in stream<br>
> pixels. I'll have to let flow through at the edges but i have no idea what<br>
<br>
</div>You can use the river boundary condition to specify the flux in stream pixel.<br>
<div class="im"><br>
> that flow should be. Are there some ways of setting the edge conditions such<br>
> that the gw evolution in the central areas of interest are not highly<br>
> influenced?<br>
<br>
</div>You can set the boundary of interest to constant head pixel and<br>
compute the flow throw the boundary pixel using the budget option.<br>
The resulting budget raster map shows the flow from active cell into<br>
sources, sinks and constant heads in [m^3/s].<br>
<br>
Best regards<br>
Soeren<br>
</blockquote></div><br><br clear="all"><br>-- <br>Vishal K. Mehta, PhD<br>Scientist<br>Stockholm Environment Institute - US<br>133 D St Suite F<br>Davis CA 95616<br><a href="http://www.sei-us.org">www.sei-us.org</a><br>