Perirhizal watersolutes

.gitignore

@@ -64,3 +64,4 @@ tutorial/jupyter/.ipynb_checkpoints/
 # images
 *.png
 .vscode/settings.json
+

src/functional/van_genuchten.py

@@ -154,9 +154,10 @@ def matric_potential_mfp(mfp, sp):
     call create_mfp_lookup first, once for each soil parameter @param sp"""
 
 
-def create_mfp_lookup(sp, wilting_point = -16000, n = 16001):
+def create_mfp_lookup(sp, wilting_point = -16000, n = 16001, verbose = True):
     """ initializes the look up tables for soil parameter to use fast_mfp, and fast_imfp """
-    print("initializing look up tables")
+    if verbose:
+        print("initializing mfp look up tables")
     global fast_mfp
     global fast_imfp
 
@@ -176,7 +177,8 @@ def create_mfp_lookup(sp, wilting_point = -16000, n = 16001):
         imfp[i] = h_[i]
     fast_imfp[sp] = interpolate.interp1d(mfp, imfp, bounds_error = False, fill_value = (imfp[0], imfp[-1]))  #
 
-    print("done")
+    if verbose:
+        print("done")
 
 # fast_specific_moisture_storage = {}
 # fast_water_content = {}

tutorial/chapter7_coupled/example7_3_lookup.py

@@ -9,4 +9,4 @@
 filename = "hydrus_loam"
 sp = vg.Parameters(hydrus_loam)  # |\label{l73l:soil_end}|
 vg.create_mfp_lookup(sp)  # |\label{l73l:mfp}|
-peri.create_lookup_mpi("results/" + filename, sp)  # |\label{l73l:lookup}|
+peri.create_lookup("results/" + filename, sp)  # |\label{l73l:lookup}|

tutorial/chapter7_coupled/results/.gitignore

@@ -1,2 +1,4 @@
 *.vtk
 *.vtu
+*.npy
+*.npz