Repository: parflow/parflow (language: python)

Issue:
Deprecated tests from local Makefile, using CTest to run tests



Discussion / hints:

Fix patch:
diff --git a/pfsimulator/parflow_lib/nl_function_eval.c b/pfsimulator/parflow_lib/nl_function_eval.c
index bf86c8c8b..8016095d0 100644
--- a/pfsimulator/parflow_lib/nl_function_eval.c
+++ b/pfsimulator/parflow_lib/nl_function_eval.c
@@ -1068,7 +1068,7 @@ void NlFunctionEval(Vector *     pressure, /* Current pressure values */
                                               rpp[ip - sz_p] * dp[ip - sz_p], rpp[ip] * dp[ip])
                                      / viscosity;
 
-                             sep = dz * z_mult_dat[ip] / 2.0;
+                             sep = dz * z_mult_dat[ip];
 
                              lower_cond = value / sep - 0.25 * dp[ip] * gravity;
                              upper_cond = pp[ip] / sep + 0.25 * dp[ip] * gravity;
@@ -1111,7 +1111,7 @@ void NlFunctionEval(Vector *     pressure, /* Current pressure values */
                                               rpp[ip] * dp[ip], rpp[ip + sz_p] * dp[ip + sz_p])
                                      / viscosity;
 
-                             sep = dz * z_mult_dat[ip] / 2.0;
+                             sep = dz * z_mult_dat[ip];
 
                              lower_cond = (pp[ip] / sep) - 0.25 * dp[ip] * gravity * z_dir_g;
                              upper_cond = (value / sep) + 0.25 * dp[ip] * gravity * z_dir_g;


Test patch:
diff --git a/examples/correct_output/test_XPlusYPlusZ.out.press.00001.pfb b/examples/correct_output/test_XPlusYPlusZ.out.press.00001.pfb
index eaa94146d..ae8a4ed16 100644
Binary files a/examples/correct_output/test_XPlusYPlusZ.out.press.00001.pfb and b/examples/correct_output/test_XPlusYPlusZ.out.press.00001.pfb differ
diff --git a/examples/correct_output/test_XYZTPlus1.out.press.00100.pfb b/examples/correct_output/test_XYZTPlus1.out.press.00100.pfb
index 1854a39c4..71732152f 100644
Binary files a/examples/correct_output/test_XYZTPlus1.out.press.00100.pfb and b/examples/correct_output/test_XYZTPlus1.out.press.00100.pfb differ
diff --git a/examples/correct_output/test_XYZTPlus1.out.satur.00100.pfb b/examples/correct_output/test_XYZTPlus1.out.satur.00100.pfb
index 1854a39c4..71732152f 100644
Binary files a/examples/correct_output/test_XYZTPlus1.out.satur.00100.pfb and b/examples/correct_output/test_XYZTPlus1.out.satur.00100.pfb differ
diff --git a/test/correct_output/Dirichlet.out.perm_x.pfb b/test/correct_output/Dirichlet.out.perm_x.pfb
new file mode 100644
index 000000000..68d15e4ae
Binary files /dev/null and b/test/correct_output/Dirichlet.out.perm_x.pfb differ
diff --git a/test/correct_output/Dirichlet.out.perm_y.pfb b/test/correct_output/Dirichlet.out.perm_y.pfb
new file mode 100644
index 000000000..68d15e4ae
Binary files /dev/null and b/test/correct_output/Dirichlet.out.perm_y.pfb differ
diff --git a/test/correct_output/Dirichlet.out.perm_z.pfb b/test/correct_output/Dirichlet.out.perm_z.pfb
new file mode 100644
index 000000000..68d15e4ae
Binary files /dev/null and b/test/correct_output/Dirichlet.out.perm_z.pfb differ
diff --git a/test/correct_output/Dirichlet.out.press.00000.pfb b/test/correct_output/Dirichlet.out.press.00000.pfb
new file mode 100644
index 000000000..62672d3d3
Binary files /dev/null and b/test/correct_output/Dirichlet.out.press.00000.pfb differ
diff --git a/test/correct_output/Dirichlet.out.press.00001.pfb b/test/correct_output/Dirichlet.out.press.00001.pfb
new file mode 100644
index 000000000..62672d3d3
Binary files /dev/null and b/test/correct_output/Dirichlet.out.press.00001.pfb differ
diff --git a/test/correct_output/Dirichlet.out.satur.00000.pfb b/test/correct_output/Dirichlet.out.satur.00000.pfb
new file mode 100644
index 000000000..1b6c58f21
Binary files /dev/null and b/test/correct_output/Dirichlet.out.satur.00000.pfb differ
diff --git a/test/correct_output/Dirichlet.out.satur.00001.pfb b/test/correct_output/Dirichlet.out.satur.00001.pfb
new file mode 100644
index 000000000..1b6c58f21
Binary files /dev/null and b/test/correct_output/Dirichlet.out.satur.00001.pfb differ
diff --git a/test/correct_output/Dirichlet.out.velx.00000.pfb b/test/correct_output/Dirichlet.out.velx.00000.pfb
new file mode 100644
index 000000000..e927805fe
Binary files /dev/null and b/test/correct_output/Dirichlet.out.velx.00000.pfb differ
diff --git a/test/correct_output/Dirichlet.out.velx.00001.pfb b/test/correct_output/Dirichlet.out.velx.00001.pfb
new file mode 100644
index 000000000..e927805fe
Binary files /dev/null and b/test/correct_output/Dirichlet.out.velx.00001.pfb differ
diff --git a/test/correct_output/Dirichlet.out.vely.00000.pfb b/test/correct_output/Dirichlet.out.vely.00000.pfb
new file mode 100644
index 000000000..5afe73ed4
Binary files /dev/null and b/test/correct_output/Dirichlet.out.vely.00000.pfb differ
diff --git a/test/correct_output/Dirichlet.out.vely.00001.pfb b/test/correct_output/Dirichlet.out.vely.00001.pfb
new file mode 100644
index 000000000..5afe73ed4
Binary files /dev/null and b/test/correct_output/Dirichlet.out.vely.00001.pfb differ
diff --git a/test/correct_output/Dirichlet.out.velz.00000.pfb b/test/correct_output/Dirichlet.out.velz.00000.pfb
new file mode 100644
index 000000000..89921b549
Binary files /dev/null and b/test/correct_output/Dirichlet.out.velz.00000.pfb differ
diff --git a/test/correct_output/Dirichlet.out.velz.00001.pfb b/test/correct_output/Dirichlet.out.velz.00001.pfb
new file mode 100644
index 000000000..3787ccce8
Binary files /dev/null and b/test/correct_output/Dirichlet.out.velz.00001.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00001.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00001.pfb
index 6b5b53ba9..b79eea430 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00001.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00001.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00002.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00002.pfb
index a4ea99199..2537c787f 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00002.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00002.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00003.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00003.pfb
index 8c77680d6..20f6ccbec 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00003.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00003.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00004.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00004.pfb
index 26d26d445..9c617a9c4 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00004.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00004.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00005.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00005.pfb
index 8a1d22f79..c22fd68a4 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00005.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00005.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00006.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00006.pfb
index 982e2b0bb..9740da082 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00006.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00006.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00007.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00007.pfb
index dc70c495f..0604364b3 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00007.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00007.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00008.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00008.pfb
index 27223d23d..a3057927a 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00008.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00008.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00009.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00009.pfb
index 0ea505a63..cfe3d3a9f 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00009.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00009.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00010.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00010.pfb
index 0a6fd9831..f46d51cf0 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00010.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00010.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00011.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00011.pfb
index 7ef6cdc2a..2fcda3708 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00011.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00011.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00012.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00012.pfb
index 32389818e..e79d17cd8 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00012.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00012.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00013.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00013.pfb
index cc86f7419..5e8bc4720 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00013.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00013.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00014.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00014.pfb
index 4d1cd81d8..cd169ac7e 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00014.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00014.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00015.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00015.pfb
index 482c907a0..bc8025055 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00015.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00015.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00016.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00016.pfb
index c8c574be6..5bbdb8c8b 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00016.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00016.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00017.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00017.pfb
index d79ddc31b..8b987257c 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00017.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00017.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00018.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00018.pfb
index b0dbd71e8..586f6cd29 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00018.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00018.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00019.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00019.pfb
index b300236da..d1025da10 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00019.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00019.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00020.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00020.pfb
index bf9ee95db..f74bd3e7b 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00020.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00020.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00021.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00021.pfb
index 45bc43092..9a94a76a9 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00021.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00021.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00022.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00022.pfb
index 041a31d83..f5f5ff8e6 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00022.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00022.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00023.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00023.pfb
index 0c7c778b2..e7c2d8b78 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00023.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00023.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.press.00024.pfb b/test/correct_output/clm_output/pfclm_sc.out.press.00024.pfb
index 3180254e3..1939fba49 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.press.00024.pfb and b/test/correct_output/clm_output/pfclm_sc.out.press.00024.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00001.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00001.pfb
index 861f17a63..a47d431ab 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00001.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00001.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00002.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00002.pfb
index f2b319f75..a28fe68cf 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00002.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00002.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00003.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00003.pfb
index 4dcc944fb..07c2ca5f9 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00003.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00003.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00004.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00004.pfb
index d06caaeac..1c0f15156 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00004.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00004.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00005.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00005.pfb
index 4362c5e91..97d52df58 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00005.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00005.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00006.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00006.pfb
index 7c33802df..4f1d35f58 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00006.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00006.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00007.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00007.pfb
index 2635401c1..ab9b09d5e 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00007.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00007.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00008.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00008.pfb
index 8069c194d..33c25625a 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00008.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00008.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00009.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00009.pfb
index afa0eb31b..7d5bdde46 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00009.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00009.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00010.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00010.pfb
index 2e44cdca6..7eaec0274 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00010.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00010.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00011.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00011.pfb
index 7d09143a9..b79b743ab 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00011.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00011.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00012.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00012.pfb
index a33a72409..2ea544cbb 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00012.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00012.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00013.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00013.pfb
index e2152fd21..24aa0a978 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00013.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00013.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00014.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00014.pfb
index 4cf560f6d..66197cef1 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00014.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00014.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00015.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00015.pfb
index 89f337eeb..d0b575e71 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00015.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00015.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00016.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00016.pfb
index b9eafd527..d37f2bbdf 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00016.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00016.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00017.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00017.pfb
index 13906fcae..55f913440 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00017.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00017.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00018.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00018.pfb
index d35beb8fa..e2707e613 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00018.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00018.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00019.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00019.pfb
index d7b7d72af..6c7723276 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00019.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00019.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00020.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00020.pfb
index 8d8eee71a..16fce5590 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00020.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00020.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00021.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00021.pfb
index b2e5da18a..b9825e00e 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00021.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00021.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00022.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00022.pfb
index ecb4f1a93..14d26165e 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00022.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00022.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00023.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00023.pfb
index d13dec842..758615285 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00023.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00023.pfb differ
diff --git a/test/correct_output/clm_output/pfclm_sc.out.satur.00024.pfb b/test/correct_output/clm_output/pfclm_sc.out.satur.00024.pfb
index feda12551..c716c1a04 100644
Binary files a/test/correct_output/clm_output/pfclm_sc.out.satur.00024.pfb and b/test/correct_output/clm_output/pfclm_sc.out.satur.00024.pfb differ
diff --git a/test/tcl/CMakeLists.txt b/test/tcl/CMakeLists.txt
index 34463e60b..c03c090b5 100644
--- a/test/tcl/CMakeLists.txt
+++ b/test/tcl/CMakeLists.txt
@@ -2,6 +2,7 @@
 include(ParflowTest)
 
 set(TESTS
+  Dirichlet.tcl
   default_single.tcl
   default_richards_wells.tcl
   forsyth2.tcl
diff --git a/test/tcl/Dirichlet.tcl b/test/tcl/Dirichlet.tcl
new file mode 100644
index 000000000..fcc7a8072
--- /dev/null
+++ b/test/tcl/Dirichlet.tcl
@@ -0,0 +1,339 @@
+#
+# Import the ParFlow TCL package
+#
+lappend auto_path $env(PARFLOW_DIR)/bin 
+package require parflow
+namespace import Parflow::*
+
+
+#-----------------------------------------------------------------------------
+# File input version number
+#-----------------------------------------------------------------------------
+pfset FileVersion 4
+
+#-----------------------------------------------------------------------------
+# Process Topology
+#-----------------------------------------------------------------------------
+
+pfset Process.Topology.P        [lindex $argv 0]
+pfset Process.Topology.Q        [lindex $argv 1]
+pfset Process.Topology.R        [lindex $argv 2]
+
+#-----------------------------------------------------------------------------
+# Computational Grid
+#-----------------------------------------------------------------------------
+pfset ComputationalGrid.Lower.X                 0.0
+pfset ComputationalGrid.Lower.Y                 0.0
+pfset ComputationalGrid.Lower.Z                 0.0
+
+pfset ComputationalGrid.DX	                1.
+pfset ComputationalGrid.DY                      1.
+pfset ComputationalGrid.DZ	                0.1
+
+pfset ComputationalGrid.NX                      1
+pfset ComputationalGrid.NY                      1
+pfset ComputationalGrid.NZ                     10 
+
+#-----------------------------------------------------------------------------
+# The Names of the GeomInputs
+#-----------------------------------------------------------------------------
+pfset GeomInput.Names "domain_input"
+
+
+#-----------------------------------------------------------------------------
+# Domain Geometry Input
+#-----------------------------------------------------------------------------
+pfset GeomInput.domain_input.InputType            Box
+pfset GeomInput.domain_input.GeomName             domain
+
+#-----------------------------------------------------------------------------
+# Domain Geometry
+#-----------------------------------------------------------------------------
+pfset Geom.domain.Lower.X                        0.0 
+pfset Geom.domain.Lower.Y                        0.0
+pfset Geom.domain.Lower.Z                        0.0
+
+pfset Geom.domain.Upper.X                        1.0
+pfset Geom.domain.Upper.Y                        1.0
+pfset Geom.domain.Upper.Z                        1.0
+
+pfset Geom.domain.Patches "left right front back bottom top"
+
+#-----------------------------------------------------------------------------
+# Perm
+#-----------------------------------------------------------------------------
+pfset Geom.Perm.Names "domain"
+
+pfset Geom.domain.Perm.Type     Constant
+pfset Geom.domain.Perm.Value    4.0
+
+pfset Perm.TensorType               TensorByGeom
+
+pfset Geom.Perm.TensorByGeom.Names  "domain"
+
+pfset Geom.domain.Perm.TensorValX  1.0
+pfset Geom.domain.Perm.TensorValY  1.0
+pfset Geom.domain.Perm.TensorValZ  1.0
+
+#-----------------------------------------------------------------------------
+# Specific Storage
+#-----------------------------------------------------------------------------
+# specific storage does not figure into the impes (fully sat) case but we still
+# need a key for it
+
+pfset SpecificStorage.Type            Constant
+pfset SpecificStorage.GeomNames       ""
+pfset Geom.domain.SpecificStorage.Value 1.0e-4
+
+#-----------------------------------------------------------------------------
+# Phases
+#-----------------------------------------------------------------------------
+
+pfset Phase.Names "water"
+
+pfset Phase.water.Density.Type	Constant
+pfset Phase.water.Density.Value	1.0
+
+pfset Phase.water.Viscosity.Type	Constant
+pfset Phase.water.Viscosity.Value	1.0
+
+#-----------------------------------------------------------------------------
+# Contaminants
+#-----------------------------------------------------------------------------
+pfset Contaminants.Names			""
+
+#-----------------------------------------------------------------------------
+# Retardation
+#-----------------------------------------------------------------------------
+pfset Geom.Retardation.GeomNames           ""
+
+#-----------------------------------------------------------------------------
+# Gravity
+#-----------------------------------------------------------------------------
+
+pfset Gravity				1.0
+
+#-----------------------------------------------------------------------------
+# Setup timing info
+#-----------------------------------------------------------------------------
+
+pfset TimingInfo.BaseUnit		1.0
+pfset TimingInfo.StartCount		0
+pfset TimingInfo.StartTime		0.0
+pfset TimingInfo.StopTime               1.0
+pfset TimingInfo.DumpInterval	       -1
+pfset TimeStep.Type                     Constant
+pfset TimeStep.Value                    1.
+
+#-----------------------------------------------------------------------------
+# Porosity
+#-----------------------------------------------------------------------------
+
+pfset Geom.Porosity.GeomNames          domain
+
+pfset Geom.domain.Porosity.Type    Constant
+pfset Geom.domain.Porosity.Value   0.5
+
+#-----------------------------------------------------------------------------
+# Domain
+#-----------------------------------------------------------------------------
+pfset Domain.GeomName domain
+
+#-----------------------------------------------------------------------------
+# Relative Permeability
+#-----------------------------------------------------------------------------
+
+pfset Phase.RelPerm.Type               VanGenuchten
+pfset Phase.RelPerm.GeomNames          domain
+pfset Geom.domain.RelPerm.Alpha        1.
+pfset Geom.domain.RelPerm.N            2.
+
+#---------------------------------------------------------
+# Saturation
+#---------------------------------------------------------
+
+pfset Phase.Saturation.Type            VanGenuchten
+pfset Phase.Saturation.GeomNames       domain
+pfset Geom.domain.Saturation.Alpha     1.
+pfset Geom.domain.Saturation.N         2.
+pfset Geom.domain.Saturation.SRes      0.2
+pfset Geom.domain.Saturation.SSat      1.
+
+#-----------------------------------------------------------------------------
+# Wells
+#-----------------------------------------------------------------------------
+pfset Wells.Names                           ""
+
+
+#-----------------------------------------------------------------------------
+# Time Cycles
+#-----------------------------------------------------------------------------
+pfset Cycle.Names constant
+pfset Cycle.constant.Names		"alltime"
+pfset Cycle.constant.alltime.Length	 1
+pfset Cycle.constant.Repeat		-1
+
+#-----------------------------------------------------------------------------
+# Boundary Conditions: Pressure
+#-----------------------------------------------------------------------------
+pfset BCPressure.PatchNames "left right front back bottom top"
+
+pfset Patch.left.BCPressure.Type			DirEquilRefPatch
+pfset Patch.left.BCPressure.Cycle			"constant"
+pfset Patch.left.BCPressure.RefGeom			domain
+pfset Patch.left.BCPressure.RefPatch			bottom
+pfset Patch.left.BCPressure.alltime.Value		0.5
+
+pfset Patch.right.BCPressure.Type			DirEquilRefPatch
+pfset Patch.right.BCPressure.Cycle			"constant"
+pfset Patch.right.BCPressure.RefGeom			domain
+pfset Patch.right.BCPressure.RefPatch			bottom
+pfset Patch.right.BCPressure.alltime.Value		0.5
+
+pfset Patch.front.BCPressure.Type			DirEquilRefPatch
+pfset Patch.front.BCPressure.Cycle			"constant"
+pfset Patch.front.BCPressure.RefGeom			domain
+pfset Patch.front.BCPressure.RefPatch			bottom
+pfset Patch.front.BCPressure.alltime.Value		0.5
+
+pfset Patch.back.BCPressure.Type			DirEquilRefPatch
+pfset Patch.back.BCPressure.Cycle			"constant"
+pfset Patch.back.BCPressure.RefGeom			domain
+pfset Patch.back.BCPressure.RefPatch			bottom
+pfset Patch.back.BCPressure.alltime.Value		0.5
+
+pfset Patch.bottom.BCPressure.Type			DirEquilRefPatch
+pfset Patch.bottom.BCPressure.Cycle			"constant"
+pfset Patch.bottom.BCPressure.RefGeom			domain
+pfset Patch.bottom.BCPressure.RefPatch			bottom
+pfset Patch.bottom.BCPressure.alltime.Value		0.5
+
+pfset Patch.top.BCPressure.Type				DirEquilRefPatch
+pfset Patch.top.BCPressure.Cycle			"constant"
+pfset Patch.top.BCPressure.RefGeom			domain
+pfset Patch.top.BCPressure.RefPatch			top
+pfset Patch.top.BCPressure.alltime.Value		-0.5
+
+#---------------------------------------------------------
+# Topo slopes in x-direction
+#---------------------------------------------------------
+# topo slopes do not figure into the impes (fully sat) case but we still
+# need keys for them
+
+pfset TopoSlopesX.Type "Constant"
+pfset TopoSlopesX.GeomNames ""
+
+pfset TopoSlopesX.Geom.domain.Value 0.0
+
+#---------------------------------------------------------
+# Topo slopes in y-direction
+#---------------------------------------------------------
+
+pfset TopoSlopesY.Type "Constant"
+pfset TopoSlopesY.GeomNames ""
+
+pfset TopoSlopesY.Geom.domain.Value 0.0
+
+#---------------------------------------------------------
+# Mannings coefficient 
+#---------------------------------------------------------
+# mannings roughnesses do not figure into the impes (fully sat) case but we still
+# need a key for them
+
+pfset Mannings.Type "Constant"
+pfset Mannings.GeomNames ""
+pfset Mannings.Geom.domain.Value 0.
+
+#---------------------------------------------------------
+# Initial conditions: water pressure
+#---------------------------------------------------------
+
+pfset ICPressure.Type                                   HydroStaticPatch
+pfset ICPressure.GeomNames                              domain
+pfset Geom.domain.ICPressure.Value                      0.5
+pfset Geom.domain.ICPressure.RefGeom                    domain
+pfset Geom.domain.ICPressure.RefPatch                   bottom
+
+#-----------------------------------------------------------------------------
+# Phase sources:
+#-----------------------------------------------------------------------------
+
+pfset PhaseSources.water.Type                         Constant
+pfset PhaseSources.water.GeomNames                    domain
+pfset PhaseSources.water.Geom.domain.Value        0.0
+
+#-----------------------------------------------------------------------------
+# Exact solution specification for error calculations
+#-----------------------------------------------------------------------------
+
+pfset KnownSolution                                    NoKnownSolution
+
+#-----------------------------------------------------------------------------
+# Set solver parameters
+#-----------------------------------------------------------------------------
+pfset Solver                                             Richards
+pfset Solver.MaxIter                                     5
+
+pfset Solver.Nonlinear.MaxIter                           10
+pfset Solver.Nonlinear.ResidualTol                       1e-9
+pfset Solver.Nonlinear.EtaChoice                         EtaConstant
+pfset Solver.Nonlinear.EtaValue                          1e-5
+pfset Solver.Nonlinear.UseJacobian                       True
+pfset Solver.Nonlinear.DerivativeEpsilon                 1e-2
+
+pfset Solver.Linear.KrylovDimension                      10
+
+pfset Solver.Linear.Preconditioner                       PFMG
+pfset Solver.Linear.Preconditioner.MGSemi.MaxIter        1
+pfset Solver.Linear.Preconditioner.MGSemi.MaxLevels      100
+
+pfset Solver.PrintVelocities True
+
+
+#-----------------------------------------------------------------------------
+# Run and Unload the ParFlow output files
+#-----------------------------------------------------------------------------
+pfrun Dirichlet
+pfundist Dirichlet
+
+
+# Tests 
+#
+source pftest.tcl
+set passed 1
+
+if ![pftestFile Dirichlet.out.perm_x.pfb "Max difference in perm_x" $sig_digits] {
+    set passed 0
+}
+if ![pftestFile Dirichlet.out.perm_y.pfb "Max difference in perm_y" $sig_digits] {
+    set passed 0
+}
+if ![pftestFile Dirichlet.out.perm_z.pfb "Max difference in perm_z" $sig_digits] {
+    set passed 0
+}
+
+foreach i "00000 00001" {
+    if ![pftestFile Dirichlet.out.press.$i.pfb "Max difference in Pressure for timestep $i" $sig_digits] {
+    set passed 0
+    }
+    if ![pftestFile Dirichlet.out.satur.$i.pfb "Max difference in Saturation for timestep $i" $sig_digits] {
+    set passed 0
+    }
+#use abs value test to prevent machine precision effects
+   set abs_value 1e-12
+   if ![pftestFileWithAbs Dirichlet.out.velx.$i.pfb "Max difference in x-velocity for timestep $i" $sig_digits $abs_value] {
+   set passed 0
+   }
+   if ![pftestFileWithAbs Dirichlet.out.vely.$i.pfb "Max difference in y-velocity for timestep $i" $sig_digits $abs_value] {
+   set passed 0
+   }
+   if ![pftestFileWithAbs Dirichlet.out.velz.$i.pfb "Max difference in z-velocity for timestep $i" $sig_digits $abs_value] {
+   set passed 0
+   }
+}
+
+if $passed {
+   puts "Dirichlet : PASSED"
+} {
+   puts "Dirichlet : FAILED"
+}