diff --git a/src/simulation/m_compute_levelset.fpp b/src/simulation/m_compute_levelset.fpp
index 9e7519790c..1f88456399 100644
--- a/src/simulation/m_compute_levelset.fpp
+++ b/src/simulation/m_compute_levelset.fpp
@@ -166,9 +166,9 @@ contains
                     end if
                 end if
             end do
-            dist_vec(1) = xy_local(1) - airfoil_grid_l(global_id)%x
-            dist_vec(2) = xy_local(2) - airfoil_grid_l(global_id)%y
-            dist_vec(3) = 0
+            dist_vec(1) = airfoil_grid_l(global_id)%x - xy_local(1)
+            dist_vec(2) = airfoil_grid_l(global_id)%y - xy_local(2)
+            dist_vec(3) = 0._wp
             dist = global_dist
         end if
 
@@ -230,14 +230,14 @@ contains
                     end if
                 end if
             end do
-            dist_vec(1) = xyz_local(1) - airfoil_grid_u(global_id)%x
-            dist_vec(2) = xyz_local(2) - airfoil_grid_u(global_id)%y
+            dist_vec(1) = airfoil_grid_u(global_id)%x - xyz_local(1)
+            dist_vec(2) = airfoil_grid_u(global_id)%y - xyz_local(2)
             dist_vec(3) = 0._wp
             dist_surf = global_dist
         else
             do k = 1, Np
-                dist_vec(1) = xyz_local(1) - airfoil_grid_l(k)%x
-                dist_vec(2) = xyz_local(2) - airfoil_grid_l(k)%y
+                dist_vec(1) = airfoil_grid_l(k)%x - xyz_local(1)
+                dist_vec(2) = airfoil_grid_l(k)%y - xyz_local(2)
                 dist_vec(3) = 0
                 dist_surf = sqrt(sum(dist_vec**2))
                 if (k == 1) then
@@ -250,8 +250,8 @@ contains
                     end if
                 end if
             end do
-            dist_vec(1) = xyz_local(1) - airfoil_grid_l(global_id)%x
-            dist_vec(2) = xyz_local(2) - airfoil_grid_l(global_id)%y
+            dist_vec(1) = airfoil_grid_l(global_id)%x - xyz_local(1)
+            dist_vec(2) = airfoil_grid_l(global_id)%y - xyz_local(2)
             dist_vec(3) = 0._wp
             dist_surf = global_dist
         end if