Aquifer diffusivity - dormant bug

[MattBrst]

In the calibration of the 1D models the aquifer diffusivity is relevant.

According to Kruseman and De Ridder this should be calculated as $aquifer diffusivity = kD / S$ with dimension $L^2/T$. So S should be dimensionless.

My AI assistant found a dormant bug in Timflow in line where $Saq$ is used with dimension $1/L$:

timflow/timflow/transient/aquifer.py

Line 87 in 882bc67

# self.D = self.T / self.Saq

For future reference:
$aquifer diffusivity=kD/S = (kaq * Haq) / (Saq * Haq) = kaq / Saq$

[dbrakenhoff]

I believe that line can be removed. The aquifer diffusivity is computed in AquiferData.initialize() as self.D = self.T / self.Scoefaq where Scoefaq is the dimensionless storage coefficient.

[mbakker7]

I agree with @dbrakenhoff. The diffusivity is indeed computed correctly as self.D = self.T / self.Scoefaq. The reported line is already commented out and can be removed.