Hi Soren,<br><br>Thanks much for your example on how to use the river parameters. i'll see about switching to Grass 7.<br><br>If you should be interested in collaborating in this kind of complex application, do let me know. i'll send you a conference paper i'm writing soon. after that, if you find it interesting enough and have the time, we could perhaps develop it together and i can put you down as co-author if that works out..<br>
<br>thanks for your patience with my repeated emails.FYI  i have done distributed surface hydrological modeling 10 years ago in GRASS, and eventually i'd like to do combined surface-groundwater modeling in GRASS. This is the first time i've attempted groundwater modeling in GRASS.<br>
<br>Vishal<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2012/7/19 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>>:<br>
> Thanks Soren,<br>
<div class="im">> for your script.<br>
><br>
> i'm using r.gwflow in v6.5svn, which does not have the budget parameter<br>
> option. do you know how i can get the 7 version to work in my 6.5 version?<br>
<br>
</div>Unfortunately the grass7 version of r.gwflow does not run in grass6.5<br>
and i have no plans o back port the budget feature.<br>
So maybe you would like to switch to grass7 to compute the budget?<br>
<div class="im"><br>
> also, i'd like to know how all th rest if it is set up - the river_bed<br>
> parameters etc. do you have complete documentation of it? i'm trying to<br>
<br>
</div>You need three parameters to model the river:<br>
1.) The water table of the river in [m], mostly derived from a river<br>
vector map and converted into a raster map<br>
2.) The bed of the river in [m], mostly derived from the river water<br>
table: r.mapcalc "river_bed = river_table - 2"<br>
3.) The leakage coefficient of the river bed in [1/s]<br>
<br>
River table, river bed and leakage are occupying the same pixels.<br>
<br>
Simple and stupid example:<br>
<br>
# River with a water table of 20m and the bed at 18m with a<br>
# leakage coefficient of 0.0001 1/s<br>
<br>
# First convert the stream network into a raster map<br>
v.to.rast input=streams output=streams type=line<br>
<br>
# Then compute the water table of the river at 20m<br>
# and the height of the bed that is 2m lower than<br>
# the water table of the river<br>
r.mapcalc "river_head = if(isnull(streams), null(), 20)"<br>
r.mapcalc "river_bed = river_head - 2"<br>
<br>
# The river leakage coefficient can be computed from the hydraulic conductivity<br>
# of the river bed (0.0001m/s) and the thickness of the river bed (1m)<br>
r.mapcalc "river_leakage = if(isnull(streams), null(), 0.0001/1)"<br>
<div class="im"><br>
> model the coupled water supply-extraction and groundwater system of a city,<br>
> so more complex real world examples of the use of r.gwflow would be really,<br>
> really useful. have you used it in a ocomplex setting perhaps for your<br>
> dissertation?<br>
<br>
</div>The only documentation that is available is the r.gwflow manual page<br>
and my diploma thesis [1] (in German)<br>
that describes the mathematical details. I never utilized r.gwflow to<br>
compute large complex real world problems.<br>
<br>
Best regards<br>
Soeren<br>
<br>
[1] <a href="http://www.hydrogeologie.tu-berlin.de/fileadmin/fg66/_hydro/Diplomarbeiten/2007_Diplomarbeit_Soeren_Gebbert.pdf" target="_blank">http://www.hydrogeologie.tu-berlin.de/fileadmin/fg66/_hydro/Diplomarbeiten/2007_Diplomarbeit_Soeren_Gebbert.pdf</a><br>

<div class="HOEnZb"><div class="h5"><br>
><br>
> thanks again,<br>
> Vishal<br>
><br>
><br>
> On Wed, Jul 18, 2012 at 2:47 PM, Sören Gebbert<br>
> <<a href="mailto:soerengebbert@googlemail.com">soerengebbert@googlemail.com</a>> wrote:<br>
>><br>
>> Hi Vishal,<br>
>><br>
>> 2012/7/18 Vishal Mehta <<a href="mailto:vishalm1975@gmail.com">vishalm1975@gmail.com</a>>:<br>
>> > hi Soren,<br>
>> ><br>
>> > my constan head boundary conditions on the edges are causing water<br>
>> > tables to<br>
>> > build up everywhere else. so i estimated a constant flux 0f 0.0019 m3/s<br>
>> > that<br>
>> > i want to apply for edge cells. this is what i am trying in order to<br>
>> > 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<br>
>> > $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<br>
>> > 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<br>
>> > 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>
>> 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<br>
>> area?<br>
>><br>
>> JFYI i am using r.gwflow of GRASS 7.<br>
>><br>
>> Best regards<br>
>> Soeren<br>
>><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>><br>
>> > 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<br>
>> >> 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<br>
>> >>> > 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<br>
>> >>> > 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<br>
>> >>> > conditions<br>
>> >>> > such<br>
>> >>> > that the gw evolution in the central areas of interest are not<br>
>> >>> > 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>
><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>