@@ -81,8 +81,8 @@ contains
8181
8282 x_centroid = patch_ib(ib_patch_id)%x_centroid
8383 y_centroid = patch_ib(ib_patch_id)%y_centroid
84- inverse_rotation(:, :) = patch_ib(patch_id )%rotation_matrix_inverse(:, :)
85- rotation(:, :) = patch_ib(patch_id )%rotation_matrix(:, :)
84+ inverse_rotation(:, :) = patch_ib(ib_patch_id )%rotation_matrix_inverse(:, :)
85+ rotation(:, :) = patch_ib(ib_patch_id )%rotation_matrix(:, :)
8686
8787 do i = 0 , m
8888 do j = 0 , n
@@ -182,7 +182,7 @@ contains
182182
183183 if (xyz_local(2 ) >= y_centroid) then
184184 do k = 1 , Np
185- dist_vec(1 ) = xzy_local (1 ) - airfoil_grid_u(k)%x
185+ dist_vec(1 ) = xyz_local (1 ) - airfoil_grid_u(k)%x
186186 dist_vec(2 ) = xyz_local(2 ) - airfoil_grid_u(k)%y
187187 dist_vec(3 ) = 0
188188 dist_surf = sqrt (sum (dist_vec** 2 ))
@@ -364,7 +364,7 @@ contains
364364 do j = 0 , n
365365 do k = 0 , p
366366
367- xyz_local = [x_cc(i) - x_centroid, cart_y - y_centroid, cart_z - z_centroid] ! get coordinate frame centered on IB
367+ xyz_local = [x_cc(i) - x_centroid, y_cc(j) - y_centroid, z_cc(k) - z_centroid] ! get coordinate frame centered on IB
368368 xyz_local = matmul (inverse_rotation, xyz_local) ! rotate the frame into the IB' s coordiantes
369369
370370 if ((xyz_local(1) > Left .and. xyz_local(1) < Right) .or. &
@@ -498,7 +498,7 @@ contains
498498 do i = 0, m
499499 do j = 0, n
500500 do k = 0, p
501- xyz_local = [x_cc(i) - x_centroid, cart_y - y_centroid, cart_z - z_centroid] ! get coordinate frame centered on IB
501+ xyz_local = [x_cc(i) - x_centroid, y_cc(j) - y_centroid, z_cc(k) - z_centroid] ! get coordinate frame centered on IB
502502 xyz_local = matmul(inverse_rotation, xyz_local) ! rotate the frame into the IB' s coordiantes
503503
504504 ! get distance to flat edge of cylinder
0 commit comments