rhof 1e3[kg/m^3] "Fluid density"
poro_1 0.339 "Porosity, layer 1"
poro_2 0.399 "Porosity, layer 2"
thetar_1 0.001 "Residual saturation, layer 1"
thetar_2 0.001 "Residual saturation, layer 2"
Sp_1 0.339[1/m]/(g_const*rhof) "Storage coefficient, layer 1"
Sp_2 0.399[1/m]/(g_const*rhof) "Storage coefficient, layer 2"
Ks_1 0.454[m/d] "Saturated hydraulic conductivity, layer 1"
Ks_2 0.298[m/d] "Saturated hydraulic conductivity, layer 2"
alpha_1 1.39[1/m] "Van Genuchten parameter alpha, layer 1"
alpha_2 1.74[1/m] "Van Genuchten parameter alpha, layer 2"
n_1 1.6 "Van Genuchten parameter n, layer 1"
n_2 1.38 "Van Genuchten parameter n, layer 2"
Hp0 0.01[m] "Pressure head in the ring"
N0 0.01*Ks_1*rhof "Leak from the base"
c0 1[mol/m^3] "Solute concentration in the ring"
rhob 1300[kg/m^3] "Bulk density"
kp_a 1e-4[m^3/kg] "Partition coefficient, aldicarb"
kp_asx 0.5e-4[m^3/kg] "Partition coefficient, aldicarb sulfoxide"
kp_asn 2e-4[m^3/kg] "Partition coefficient, aldicarb sulfone"
Dl 3.74e-3[m^2/d] "Diffusion coefficient, liquid phase"
Dg 0.432[m^2/d] "Diffusion coefficient, gas phase"
alphar 0.005[m] "Longitudinal dispersivity"
alphaz 0.001[m] "Transverse dispersivity"
kg_a 1.33e-7 Volatilization
kg_asn 1.33e-3 Volatilization
d_s 0.05[m] "Thickness of the concentration boundary layer"