Fixed sign error

Author: danieljvickers

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
@@ -378,8 +378,8 @@ contains
         xy_local = matmul(inverse_rotation, xy_local)
 
         normal_vector = xy_local
-        normal_vector(2) = normal_vector(2)*(ellipse_coeffs(1)/ellipse_coeffs(2)) &
-                      & **2._wp  ! get the normal direction via the coordinate transformation method
+        ! get the normal direction via the coordinate transformation method
+        normal_vector(2) = normal_vector(2)*(ellipse_coeffs(1)/ellipse_coeffs(2))**2._wp
         normal_vector = normal_vector/sqrt(dot_product(normal_vector, normal_vector))  ! normalize the vector
         gp%levelset_norm = matmul(rotation, normal_vector)  ! save after rotating the vector to the global frame