#---------------------------------------------------------------------------------
~Simulation Title Card
#---------------------------------------------------------------------------------
1,
STOMP Example Problem W-2,
Mart Oostrom,
PNNL,
June 03,
15:00,
4,
Two aquifer problem (Segol 1994, p. 423-432),
Two aquifers (5 m thick) separated by an aquitard (10 m thick),
Pumping well is screened in the lower aquifer only,
No drawdown allowed on other end,


#---------------------------------------------------------------------------------
~Solution Control Card
#---------------------------------------------------------------------------------
Normal,
Water,
1,
0,s,100,yr,36,s,50,d,1.414,24,1.e-6,
10000,
0,


#---------------------------------------------------------------------------------
~Grid Card
#---------------------------------------------------------------------------------
Cylindrical,
13,1,20,
1.0,m,1.195,m,2.365,m,3.955,m,7.285,m,12.715,m,22.885,m,40.315,m,72.085,m,127.915,m,
504.485,m,1495.515,m,4828.485,m,10000,m,
0.,deg,18.,deg,
0,m,20@1.0,m,


#---------------------------------------------------------------------------------
~Rock/Soil Zonation Card
#---------------------------------------------------------------------------------
3,
Aquifer1,1,13,1,1,16,20,
Aquitard,1,13,1,1,6,15,
Aquifer2,1,13,1,1,1,5,


#---------------------------------------------------------------------------------
~Mechanical Properties Card
#---------------------------------------------------------------------------------
Aquifer1,2650,kg/m^3,0.30,0.30,1.e-04,1/m,Millington and Quirk,
Aquitard,2650,kg/m^3,0.43,0.43,8.e-04,1/m,Millington and Quirk,
Aquifer2,2650,kg/m^3,0.30,0.30,1.e-04,1/m,Millington and Quirk,


#---------------------------------------------------------------------------------
~Hydraulic Properties Card
#---------------------------------------------------------------------------------
Aquifer1,2.e-05,hc m/s,,,2.e-05,hc m/s,
Aquitard,1.e-08,hc m/s,,,1.e-08,hc m/s,
Aquifer2,2.e-05,hc m/s,,,2.e-05,hc m/s,


#---------------------------------------------------------------------------------
~Saturation Function Card
#---------------------------------------------------------------------------------
Aquifer1,van Genuchten,0.133,1/cm,1.88,0.268,,
Aquitard,van Genuchten,0.133,1/cm,1.88,0.268,,
Aquifer2,van Genuchten,0.133,1/cm,1.88,0.268,,


#---------------------------------------------------------------------------------
~Aqueous Relative Permeability Card
#---------------------------------------------------------------------------------
Aquifer1,Mualem,,
Aquitard,Mualem,,
Aquifer2,Mualem,,


#---------------------------------------------------------------------------------
~Initial Conditions Card
#---------------------------------------------------------------------------------
Gas Pressure,Aqueous Pressure, 
1,
Aqueous Pressure,684036.8,Pa,,,,,-9793.5,1/m,1,13,1,1,1,20,

#---------------------------------------------------------------------------------
~Boundary Conditions Card
#---------------------------------------------------------------------------------
2,
west,neumann,
1,1,1,1,1,5,1,
0,d,-5.02929620e-05,m/s,
east,hydraulic gradient,
13,13,1,1,1,20,1,
0,d,684036.8,Pa,


#---------------------------------------------------------------------------------
~Output Options Card
#---------------------------------------------------------------------------------
4,
3,1,3,
7,1,3,
3,1,18,
7,1,10,
1,1,yr,m,deg,6,6,6,
3,
aqueous saturation,,
aqueous pressure,Pa,
aqueous hydraulic head,m,
5,
0.5,d,
1,d,
1,yr,
10,yr,
50,yr,
8,
no restart,,
aqueous hydraulic head,m,
aqueous pressure,Pa,
aqueous saturation,,
XNC aqueous volumetric flux,m/hr,
ZNC aqueous volumetric flux,m/hr,
X aqueous volumetric flux,m/hr,
Z aqueous volumetric flux,m/hr,


#---------------------------------------------------------------------------------
~Surface Flux Card
#---------------------------------------------------------------------------------
2,
aqueous volumetric flux,m^3/day,m^3,west,1,1,1,1,1,5,
aqueous volumetric flux,m^3/day,m^3,east,13,13,1,1,1,20,