Fix ksat estimation by SWRC_PTF_Cosby1984_for_Campbell1974()

- now coefficient of the clay predictor is correct, see Table 4 of Cosby et al. 1984 - example SWRC ksat inputs are corrected - cf. commit 5653f7f "Added saturated hydraulic conductivity as parameter to swrcp" from June 10, 2022
DrylandEcology · Sep 6, 2024 · 2351511 · 2351511
1 parent 67cb40a
commit 2351511
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 9 deletions.
diff --git a/src/SW_Site.c b/src/SW_Site.c
@@ -499,7 +499,7 @@ void SWRC_PTF_Cosby1984_for_Campbell1974(
     swrcp[2] = -0.3 * sand + 15.7 * clay + 3.10;
     /* swrcp[3] = K_saturated: originally with units [inches / day]
        here re-formulated with units [cm / day] */
-    swrcp[3] = 2.54 * 24. * powe(10.0, 1.26 * sand - 6.4 * clay - 0.60);
+    swrcp[3] = 2.54 * 24. * powe(10.0, 1.26 * sand - 0.64 * clay - 0.60);
 }
 
 /**

diff --git a/tests/example/Input/swrc_params.in b/tests/example/Input/swrc_params.in
@@ -29,11 +29,11 @@
 
 
 # param1  param2   param3   param4     param5   param6
-18.6080  0.42703   5.3020   24.03047   0.0000   0.0000
-20.4644  0.43290   7.0500   14.94351   0.0000   0.0000
-22.8402  0.44013   9.4320   13.71177   0.0000   0.0000
-24.0381  0.44291  10.0690   13.43171   0.0000   0.0000
-24.2159  0.44359  10.3860   14.14351   0.0000   0.0000
-23.3507  0.44217  10.3830   40.63764   0.0000   0.0000
-12.3880  0.41370   7.3250   37.22899   0.0000   0.0000
-12.3880  0.41370   7.3250   37.22899   0.0000   0.0000
+18.6080  0.42703   5.3020   53.90697   0.0000   0.0000
+20.4644  0.43290   7.0500   37.41493   0.0000   0.0000
+22.8402  0.44013   9.4320   23.10175   0.0000   0.0000
+24.0381  0.44291  10.0690   19.96388   0.0000   0.0000
+24.2159  0.44359  10.3860   18.82977   0.0000   0.0000
+23.3507  0.44217  10.3830   19.38407   0.0000   0.0000
+12.3880  0.41370   7.3250   52.97183   0.0000   0.0000
+12.3880  0.41370   7.3250   52.97183   0.0000   0.0000