--- a/+rv/+diffops/constructSymmetricD2.m	Wed Aug 07 13:28:21 2019 +0200
+++ b/+rv/+diffops/constructSymmetricD2.m	Wed Aug 07 15:23:42 2019 +0200
@@ -20,9 +20,10 @@
                 M =  Dm - Hi*B;
                 D2 = @(Viscosity) M*Viscosity*Dp;
             else
-                % TODO: Fix closure for D2Variable
                 % TODO: Fix Viscosity not being vector
-                D2 = @(Viscosity)ops{1}.D2(diag(Viscosity));
+                d1_r = ops{1}.d1_r';
+                d1_l = ops{1}.d1_l';
+                D2 = @(Viscosity)ops{1}.D2(diag(Viscosity)) + Hi*(Viscosity(1,1)*e_l*d1_l - e_r*Viscosity(end,end)*d1_r);
             end
         case 2
             % TODO: