diff --git a/Meshes/multi_drop_gauss.geo b/Meshes/multi_drop_gauss.geo
index 80ddb2a021615c135566ae1dcf467b740f461642..a0c98d93ce2357ca5524bacfa84f60ed2e2b952b 100644
--- a/Meshes/multi_drop_gauss.geo
+++ b/Meshes/multi_drop_gauss.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
 Physical Volume("volumedomain") = {1};
 
 Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4.5)*(x-4.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-3.5)*(x-3.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4.5)*(y-4.5)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-3.5)*(y-3.5)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[8] = Min;
 Field[8].FieldsList = {3,4,5,6,7};
diff --git a/Meshes/multi_drop_gauss_far.geo b/Meshes/multi_drop_gauss_far.geo
index e9a25ee8befec81647b818e78a0e284a3fff822a..d9e4f86a251613e9c125125332e61a5535393fef 100644
--- a/Meshes/multi_drop_gauss_far.geo
+++ b/Meshes/multi_drop_gauss_far.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
 Physical Volume("volumedomain") = {1};
 
 Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-5.5)*(x-5.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-2.5)*(x-2.5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-5.5)*(y-5.5)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-2.5)*(y-2.5)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[8] = Min;
 Field[8].FieldsList = {3,4,5,6,7};
diff --git a/Meshes/multi_drop_gauss_med.geo b/Meshes/multi_drop_gauss_med.geo
index f189c78c532e8efca5f82231f275ed48b65c1a33..5a28df73f6f6ce2c9a48022d5c5155af805ef106 100644
--- a/Meshes/multi_drop_gauss_med.geo
+++ b/Meshes/multi_drop_gauss_med.geo
@@ -5,19 +5,19 @@ Physical Surface("surfacedomain") = {1};
 Physical Volume("volumedomain") = {1};
 
 Field[3] = MathEval;
-Field[3].F = ".5 - (.5-0.01)*exp(-(1./0.015) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[3].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[4] = MathEval;
-Field[4].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[4].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-5)*(x-5)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[5] = MathEval;
-Field[5].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.5)*(z-0.5))-0.25)))";
+Field[5].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-3)*(x-3)+(y-4)*(y-4)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[6] = MathEval;
-Field[6].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.5)*(z-0.5))-0.25)))";
+Field[6].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-5)*(y-5)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[7] = MathEval;
-Field[7].F = ".5 - (.5-0.01)*exp(-(1./0.025) *((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.5)*(z-0.5))-0.25))*((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.5)*(z-0.5))-0.25)))";
+Field[7].F = ".25 - .24*exp(-(1./0.25) *((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.25)*(z-0.25))-0.25))*((sqrt((x-4)*(x-4)+(y-3)*(y-3)+(z-0.25)*(z-0.25))-0.25)))";
 
 Field[8] = Min;
 Field[8].FieldsList = {3,4,5,6,7};