--- a/+scheme/LaplaceCurvilinearVirtaMin.m	Mon Jun 10 13:27:29 2019 +0200
+++ b/+scheme/LaplaceCurvilinearVirtaMin.m	Mon Jun 10 14:39:14 2019 +0200
@@ -388,7 +388,7 @@
             bp = v.bp;
 
             lambda_mat = reshape(lambda_v, my, mx);
-            switch boundary
+            switch neighbour_boundary
             case 'w'
                 lambda_min = min(lambda_mat(:,1:bp), [], 2);
                 lambda_mat = repmat(lambda_min, 1, mx);
@@ -406,20 +406,20 @@
 
             switch boundary
             case {'w', 'e'}
-                b1_u = gamm_u*u.lambda_min_u(I_u)./u.a11(I_u).^2;
-                b2_u = h11_u *u.lambda(I_u)./u.a22(I_u).^2;
+                b1_u = gamm_u*lambda_min_u(I_u)./u.a11(I_u).^2;
+                b2_u = h11_u *lambda_u(I_u)./u.a22(I_u).^2;
             case {'s', 'n'}
-                b1_u = h11_u *u.lambda(I_u)./u.a11(I_u).^2;
-                b2_u = gamm_u*u.lambda_min(I_u)./u.a22(I_u).^2;
+                b1_u = h11_u *lambda_u(I_u)./u.a11(I_u).^2;
+                b2_u = gamm_u*lambda_min_u(I_u)./u.a22(I_u).^2;
             end
 
             switch neighbour_boundary
             case {'w', 'e'}
-                b1_v = gamm_v*v.lambda_min_v(I_v)./v.a11(I_v).^2;
-                b2_v = h11_v *v.lambda(I_v)./v.a22(I_v).^2;
+                b1_v = gamm_v*lambda_min_v(I_v)./v.a11(I_v).^2;
+                b2_v = h11_v *lambda_v(I_v)./v.a22(I_v).^2;
             case {'s', 'n'}
-                b1_v = h11_v *v.lambda(I_v)./v.a11(I_v).^2;
-                b2_v = gamm_v*v.lambda_min(I_v)./v.a22(I_v).^2;
+                b1_v = h11_v *lambda_v(I_v)./v.a11(I_v).^2;
+                b2_v = gamm_v*lambda_min_v(I_v)./v.a22(I_v).^2;
             end
 
             tau1 = -1./(4*b1_u) -1./(4*b1_v) -1./(4*b2_u) -1./(4*b2_v);