g.region rows=80, cols=140, res3=1, t=100, b=0, n=800,s=0, w=0, e=1400 --overwrite

r.mapcalc expression="fracture = 0" --overwrite

r.mapcalc expression="phead = 65" --overwrite
r.mapcalc expression="top = 77" --overwrite
r.mapcalc expression="null = 0.0" --overwrite
r.mapcalc expression="status = if(col() == 1 || col() == 140 || row()==1 || row() == 80, 2, 1)" --overwrite
r.mapcalc expression="recharge= 0.2/31536000" --overwrite

r.mapcalc expression="hydcondx = if(fracture==0,0.00001,0.001)" --overwrite
r.mapcalc expression="hydcondy = if(fracture==0,0.00001,0.01)" --overwrite
r.mapcalc expression="poros = if(fracture==0,0.25,0.02)" --overwrite

r.mapcalc expression="q= if((col()== 30 && row() == 20) || (col()== 110 && row()==20) || (col()==80 && row()==60), -0.01,0)" --overwrite

# Compute a steady state groundwater flow
r.gwflow top=top bottom=null phead=phead status=status hc_x=hydcondx hc_y=hydcondy s=poros q=q output=gwresult dt=864000000 type="unconfined" vx=velx vy=vely --overwrite