Thanks Soren,<br>for your script.<br><br>i'm using r.gwflow in v6.5svn, which does not have the budget parameter option. do you know how i can get the 7 version to work in my 6.5 version?<br><br>also, i'd like to know how all th rest if it is set up - the river_bed parameters etc. do you have complete documentation of it? i'm trying to model the coupled water supply-extraction and groundwater system of a city, so more complex real world examples of the use of r.gwflow would be really, really useful. have you used it in a ocomplex setting perhaps for your dissertation?<br>
<br>thanks again,<br>Vishal<br><br><div class="gmail_quote">On Wed, Jul 18, 2012 at 2:47 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 Vishal,<br>
<br>
2012/7/18 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>>:<br>
<div class="im">> hi Soren,<br>
><br>
> my constan head boundary conditions on the edges are causing water tables to<br>
> build up everywhere else. so i estimated a constant flux 0f 0.0019 m3/s that<br>
> i want to apply for edge cells. this is what i am trying in order to impose<br>
> a constant flux on the edges<br>
><br>
> r.mapcalc "sink.init=if(row()==1 || row()==444 ||col()==1<br>
> ||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<br>
> 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<br>
> s=s.1 type=unconfined dt=2592000 error=0.05 phead=sim.$((prevmonth))<br>
> r=gwnaturalms.$((month)) q=sink.init output=sim.$((month))<br>
><br>
> my boundary condition map (bc2) has contant head in strea pixels. all else<br>
> are calculated.<br>
><br>
> do you think the above implementation will correctly impose the constant<br>
> flux on the edges? it seems though that since this is not a boundary<br>
> condition, how can it ensure a constant flux at edges?<br>
<br>
</div>Your script looks reasonable, except that the error term is much to large.<br>
You may need to use a smaller number like 10⁻7,<br>
otherwise your results my be wrong.<br>
<br>
I have attached a small r.gwflow example to show you how to estimate<br>
the constant flux at<br>
a western boundary using the budget computation. Since the flow will<br>
be specified<br>
as a source term, it will be constant the whole computational time.<br>
<br>
Maybe you can apply this method to estimate the boundary flux in your area?<br>
<br>
JFYI i am using r.gwflow of GRASS 7.<br>
<br>
Best regards<br>
Soeren<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> Thanks for any help you can provide,<br>
> Vishal<br>
><br>
> On Mon, Jul 16, 2012 at 1:46 PM, Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>> wrote:<br>
>><br>
>> Thanks Soren,<br>
>> From you response, can you please tell me how to do the following two<br>
>> tasks (which i dont find in the online manual for r.gwflow): the remainder<br>
>> 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<br>
>> now- how can i get the budget raster maps you mention?<br>
>><br>
>> thanks again,<br>
>> Vishal<br>
>><br>
>><br>
>> On Sat, Jul 14, 2012 at 12:05 PM, Sören Gebbert<br>
>> <<a href="mailto:soerengebbert@googlemail.com">soerengebbert@googlemail.com</a>> wrote:<br>
>>><br>
>>> Hi,<br>
>>> sorry for the delay.<br>
>>><br>
>>> 2012/7/11 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>>:<br>
>>> > Thanks Soren,<br>
>>> > That explains some of the results i'm getting, with water piling up<br>
>>> > 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<br>
>>> > constant<br>
>>> > head? If constant flux, should that be in m/s units?<br>
>>><br>
>>> 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<br>
>>> r.mapcalc.<br>
>>><br>
>>> ><br>
>>> > My problem though is that i dont know what a constant flux or head at<br>
>>> > the<br>
>>> > edges should be set to. For now the only bc i have put in there<br>
>>> > deliberately<br>
>>> > (beyond the default you mention) is that i have set constant head in<br>
>>> > stream<br>
>>> > pixels. I'll have to let flow through at the edges but i have no idea<br>
>>> > what<br>
>>><br>
>>> You can use the river boundary condition to specify the flux in stream<br>
>>> pixel.<br>
>>><br>
>>> > that flow should be. Are there some ways of setting the edge conditions<br>
>>> > such<br>
>>> > that the gw evolution in the central areas of interest are not highly<br>
>>> > influenced?<br>
>>><br>
>>> 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>
>><br>
>><br>
>><br>
>><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" target="_blank">www.sei-us.org</a><br>
><br>
><br>
><br>
><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" 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>