~Simulation Title Card 
1, 
STOMP Example Problem CO2E-1, 
Mark White, 
Pacific Northwest Laboratory, 
01 June 2011, 
09:37 PDT, 
8, 
2.1 Definition of benchmark problem 1: CO2 plume evolution and 
leakage through an abandoned well
2.1.1 Formulated by A. Ebigbo, J. Nordbotten and H. Class
Problem description CO2 is injected into an aquifer; spreads within the 
aquifer and, upon reaching a leaky well, rises up to a shallower 
aquifer. A quantification of the leakage rate which depends on the 
pressure build- up in the aquifer due to injection and on the plume 
evolution is the goal of the simulation.

~Solution Control Card 
Normal,
STOMP-CO2e Isothermal w/ Invariant Fluid Density and Viscosity,
1, 
0,day,1000,day,1,s,10,day,1.25,16,1.e-06,0.001,s,0.2,
10000, 
Variable Aqueous Diffusion, 
Variable Gas Diffusion, 
1045,kg/m^3,2.535e-4,Pa s,479,kg/m^3,3.950e-5,Pa s,
0, 

~Grid Card 
Cartesian,
89,1,44,
-700.000,m,-675.000,m,-650.000,m,-625.000,m,
-600.000,m,-575.000,m,-550.000,m,-525.000,m,
-500.000,m,-475.000,m,-450.000,m,-425.000,m,-400.000,m,-375.000,m,
-350.000,m,-325.000,m,-300.000,m,-275.000,m,-255.000,m,-235.000,m,
-215.000,m,-200.000,m,-185.000,m,-170.000,m,-160.000,m,-150.000,m,
-140.000,m,-130.000,m,-122.500,m,-115.000,m,-107.500,m,-102.500,m,
-97.500,m,-92.500,m,-85.000,m,-77.500,m,-70.000,m,-60.000,m,
-50.000,m,-40.000,m,-30.000,m,-20.000,m,-12.500,m,-8.000,m,
-5.000,m,-2.800,m,-1.500,m,-0.800,m,-0.400,m,-0.133,m,
0.133,m,0.400,m,0.800,m,1.500,m,2.800,m,5.000,m,
8.000,m,12.500,m,20.000,m,30.000,m,40.000,m,50.000,m,
60.000,m,70.000,m,85.000,m,100.000,m,125.000,m,150.000,m,
175.000,m,200.000,m,225.000,m,250.000,m,275.000,m,300.000,m,
325.000,m,350.000,m,375.000,m,400.000,m,425.000,m,450.000,m,
475.000,m,500.000,m,525.000,m,550.000,m,575.000,m,600.000,m,
625.000,m,650.000,m,675.000,m,700.000,m,
0.000,m,1.000,m,
0.0,m,4.41,m,8.32,m,11.80,m,14.90,m,17.65,m,20.09,m,22.27,m,24.20,m,
25.92,m,27.44,m,28.80,m,30.0,m,20@5.000,m,12@2.500,m,

~Inactive Nodes Card
2,
1,49,1,1,13,32,
51,89,1,1,13,32,

~Rock/Soil Zonation Card 
2,
aquifer,1,89,1,1,1,44,
leaky well,50,50,1,1,1,44,

~Mechanical Properties Card 
aquifer,2650,kg/m^3,0.15,0.15,Compressibility,1.e-9,1/psi,,,constant,1.0,1.0,
leaky well,2650,kg/m^3,0.15,0.15,Compressibility,1.e-9,1/psi,,,constant,1.0,1.0,

~Hydraulic Properties Card
aquifer,2.e-14,m^2,2.e-14,m^2,2.e-14,m^2,
leaky well,1.e-12,m^2,1.e-12,m^2,1.e-12,m^2,

~Saturation Function Card
aquifer,Brooks and Corey,0.1,m,2.0,,,
leaky well,Brooks and Corey,0.1,m,2.0,,,

~Aqueous Relative Permeability Card 
aquifer,tabular,2,
0.0,0.0,
1.0,1.0,
leaky well,tabular,2,
0.0,0.0,
1.0,1.0,

~Gas Relative Permeability Card 
aquifer,tabular,2,
0.0,0.0,
1.0,1.0,
leaky well,tabular,2,
0.0,0.0,
1.0,1.0,

~Salt Transport Card 
aquifer,0.0,ft,0.0,ft,
leaky well,0.0,ft,0.0,ft,

~Initial Conditions Card 
Gas Pressure,Aqueous Pressure, 
3, 
Gas Pressure,30.81955,MPa,,,,,-0.01025,1/m,1,89,1,1,1,44,
Aqueous Pressure,30.81955,MPa,,,,,-0.01025,1/m,1,89,1,1,1,44,
Temperature,34,C,,,,,,,1,89,1,1,1,44,

~Boundary Conditions Card 
2, 
West,Aqu. Initial Condition,Gas Initial Condition,Aqu. Mass Frac.,
1,1,1,1,1,44,1, 
0,s,,,,,,,0.0,,,, 
East,Aqu. Initial Condition,Gas Initial Condition,Aqu. Mass Frac.,
89,89,1,1,1,44,1, 
0,s,,,,,,,0.0,,,, 

~Coupled Well Card
1,
CO2 Injection Well,Water Relative Saturation,1.0,0.5,1.0,0.383184,MMT,
1,
-100.0,m,0.5,m,30.0,m,-100.0,m,0.5,m,0.0,m,0.15,m,0.0,screened,
1,
0.0,hr,0.025,kg/s,45,MPa,0.0,
#0.0,hr,4.435,kg/s,45,MPa,0.0,

~Output Options Card 
5, 
32,1,1,
32,1,12,
50,1,12,
50,1,22,
50,1,33,
1,1,day,m,6,6,6, 
16,
Phase Condition,,
Gas Saturation,,
Gas Relative Permeability,,
#Trapped Gas Saturation,,
Integrated CO2 Mass,kg, 
Integrated CO2 Aqueous Mass,kg, 
Integrated CO2 Gas Mass,kg, 
Integrated CO2 Trapped-Gas Mass,kg, 
Salt Aqueous Mass Fraction,, 
CO2 Aqueous Mass Fraction,, 
Gas Density,kg/m^3,
Aqueous Density,kg/m^3,
Gas Viscosity,Pa s,
Aqueous Viscosity,Pa s,
Gas Pressure,MPa, 
Aqueous Pressure,MPa, 
#Diffusive Porosity,,
Coupled-Well Press,1,MPa,
#Coupled-Well CO2 Mass Rate,1,kg/s,
#Coupled-Well CO2 Mass Integral,1,kg,
8, 
0.1,day,
0.5,day,
1.0,day,
5.0,day,
10.0,day,
50.0,day,
100.0,day,
500.0,day,
13, 
Rock/Soil Type,,
Gas Saturation,, 
Trapped Gas Saturation,,
Salt Saturation,, 
CO2 Aqueous Concentration,gm/cm^3, 
Salt Aqueous Mass Fraction,, 
CO2 Aqueous Mass Fraction,, 
Gas Pressure,MPa, 
Aqueous Pressure,MPa, 
Diffusive Porosity,, 
Gas Density,kg/m^3,
Aqueous Density,kg/m^3,
Phase Condition,,

~Surface Flux Card
1, 
CO2 Mass Flux,kg/s,kg,top,50,50,1,1,22,22,