hi Soren,<br><br>my constan head boundary conditions on the edges are causing water tables to build up everywhere else. so i estimated a constant flux 0f 0.0019 m3/s that i want to apply for edge cells. this is what i am trying in order to impose a constant flux on the edges<br>
<br>r.mapcalc "sink.init=if(row()==1 || row()==444 ||col()==1 ||col==477,-0.0019,null())<br><br>#sink.init is the constant flux on edges<br># my r.gwflow script runs at monthly time step- ihave a loop $month-here's the r.gwflow snippet<br>
<br>r.gwflow --o -s solver=cg top=top bottom=bottom status=bc2 hc_x=k.1 hc_y=k.1 s=s.1 type=unconfined dt=2592000 error=0.05 phead=sim.$((prevmonth)) r=gwnaturalms.$((month)) q=sink.init output=sim.$((month))<br><br>my boundary condition map (bc2) has contant head in strea pixels. all else are calculated.<br>
<br>do you think the above implementation will correctly impose the constant flux on the edges? it seems though that since this is not a boundary condition, how can it ensure a constant flux at edges?<br><br>Thanks for any help you can provide,<br>
Vishal<br><div class="gmail_quote">On Mon, Jul 16, 2012 at 1:46 PM, Vishal Mehta <span dir="ltr"><<a href="mailto:vishalm1975@gmail.com" target="_blank">vishalm1975@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
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<div class="HOEnZb"><div class="h5"><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:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
sorry for the delay.<br>
<br>
2012/7/11 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com" target="_blank">vishalm1975@gmail.com</a>>:<br>
<div>> 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><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><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></div></div><div class="HOEnZb"><div class="h5">-- <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" target="_blank">www.sei-us.org</a><br>

</div></div></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>