(dh/dt)*Ss = Kxx * (d^2h/dx^2) + Kyy * (d^2h/dy^2) + q
# set the region accordingly g.region res=25 res3=25 t=100 b=0 n=1000 s=0 w=0 e=1000 #now create the input raster maps for confined and unconfined aquifers r.mapcalc "phead=if(row() == 1 , 50, 40)" r.mapcalc "status=if(row() == 1 , 2, 1)" r.mapcalc "well=if(row() == 20 && col() == 20 , -0.001, 0)" r.mapcalc "hydcond=0.00025" r.mapcalc "recharge=0" r.mapcalc "top_conf=20.0" r.mapcalc "top_unconf=70.0" r.mapcalc "bottom=0.0" r.mapcalc "null=0.0" r.mapcalc "poros=0.15" r.mapcalc "syield=0.0001" #confined groundwater flow with cg solver and sparse matrix r.gwflow --o -s solver=cg top=top_conf bottom=bottom phead=phead status=status \ hc_x=hydcond hc_y=hydcond q=well s=syield r=recharge output=gwresult_conf \ dt=8640000 type=confined velocity=gwresult_conf_velocity #unconfined groundwater flow with cg solver and sparse matrix r.gwflow --o -s solver=cg top=top_unconf bottom=bottom phead=phead \ status=status hc_x=hydcond hc_y=hydcond q=well s=poros r=recharge \ output=gwresult_unconf dt=8640000 type=unconfined velocity=gwresult_unconf_velocity # The data can be visulaized with paraview when exported with r.out.vtk r.out.vtk -p in=gwresult_conf,status vector=gwresult_conf_velocity_x,gwresult_conf_velocity_y,null out=/tmp/gwdata_conf2d.vtk r.out.vtk -p elevation=gwresult_unconf in=gwresult_unconf,status vector=gwresult_unconf_velocity_x,gwresult_unconf_velocity_y,null out=/tmp/gwdata_unconf2d.vtk #now load the data into paraview paraview --data=/tmp/gwdata_conf2d.vtk & paraview --data=/tmp/gwdata_unconf2d.vtk &
Last changed: $Date$