Updating solute transport model parameters

Author: ConcettaDAmato

src/main/java/it/geoframe/blogspot/closureequation/closureequation/Parameters.java

@@ -35,7 +35,7 @@ public static Parameters getInstance(double referenceTemperatureSWRC, double bet
 		return uniqueInstance;
 	}
 
-	public static Parameters getInstance(double molecularDiffusion,double longitudinalDispersivity, double referenceTemperatureSWRC, double beta0,
+	public static Parameters getInstance(double[] molecularDiffusion,double[] longitudinalDispersivity, double referenceTemperatureSWRC, double beta0,
 			double[] thetaS, double[] thetaR, double[] par1, double[] par2, double[] par3, double[] par4, double[] par5, double[] kappaSaturation,
 			double[] alphaSpecificStorage, double[] betaSpecificStorage) {
 		if (uniqueInstance == null) {
@@ -72,8 +72,8 @@ public static Parameters getInstance(double molecularDiffusion,double longitudin
 	public double[] alphaSpecificStorage;
 	public double[] betaSpecificStorage;
 
-	public double molecularDiffusion;
-	public double longitudinalDispersivity;
+	public double[] molecularDiffusion;
+	public double[] longitudinalDispersivity;
 
 
 	private Parameters(double waterDensity, double iceDensity, double specificThermalCapacityWater,
@@ -139,12 +139,12 @@ private Parameters(double referenceTemperatureSWRC, double beta0, double[] theta
 
 	}
 
-	private Parameters(double molecularDiffusion,double longitudinalDispersivity, double referenceTemperatureSWRC, double beta0,
+	private Parameters(double[] molecularDiffusion,double[] longitudinalDispersivity, double referenceTemperatureSWRC, double beta0,
 			double[] thetaS, double[] thetaR, double[] par1, double[] par2, double[] par3, double[] par4, double[] par5, double[] kappaSaturation,
 			double[] alphaSpecificStorage, double[] betaSpecificStorage) {
 
-		this.molecularDiffusion = molecularDiffusion;
-		this.longitudinalDispersivity = longitudinalDispersivity;
+		this.molecularDiffusion = molecularDiffusion.clone();
+		this.longitudinalDispersivity = longitudinalDispersivity.clone();
 		this.referenceTemperatureSWRC = referenceTemperatureSWRC;
 		this.beta0 = beta0;
 		this.thetaS = thetaS.clone();