Fixed airfoil shape, but it is not rotating yet

Author: danieljvickers

src/common/m_compute_levelset.fpp

@@ -90,8 +90,8 @@ contains
                 xy_local = matmul(inverse_rotation, xy_local) ! rotate the frame into the IB's coordiantes
 
                 if (xy_local(2) >= 0._wp) then
+                    ! finds the location on the airfoil grid with the minimum distance (closest)
                     do k = 1, Np
-                        ! TODO :: I think that witht he IBM implementation that airfoil_grid_u(k).y = 0._wp for all indices. Delete?
                         dist_vec(1) = xy_local(1) - airfoil_grid_u(k)%x
                         dist_vec(2) = xy_local(2) - airfoil_grid_u(k)%y
                         dist_vec(3) = 0
@@ -111,6 +111,7 @@ contains
                     dist_vec(3) = 0
                     dist = global_dist
                 else
+                    ! TODO :: This looks identical to the above code but using the lower array. Refactor this.
                     do k = 1, Np
                         dist_vec(1) = xy_local(1) - airfoil_grid_l(k)%x
                         dist_vec(2) = xy_local(2) - airfoil_grid_l(k)%y

src/common/m_ib_patches.fpp

@@ -209,6 +209,7 @@ contains
         allocate (airfoil_grid_u(1:Np)) 
         allocate (airfoil_grid_l(1:Np))
 
+        ! TODO :: The below instantiations are already handles by the loop below
         airfoil_grid_u(1)%x = 0._wp
         airfoil_grid_u(1)%y = 0._wp
 
@@ -218,6 +219,7 @@ contains
         eta = 1._wp
 
         do i = 1, Np1 + Np2 - 1
+          ! TODO :: This allcoated the upper and lower airfoil arrays, and does not need to be performed each time the IB markers are updated. Place this as a separate subroutine.
             if (i <= Np1) then
                 xc = i*(pa*ca_in/Np1)
                 xa = xc/ca_in
@@ -274,19 +276,19 @@ contains
                         yc = (ma/(1 - pa)**2)*(1 - 2*pa + 2*pa*xa - xa**2)
                         dycdxc = (2*ma/(1 - pa)**2)*(pa - xa)
                     end if
-                    if (xy_local(1) >= 0._wp) then
+                    if (xy_local(2) >= 0._wp) then
                         k = 1
                         do while (airfoil_grid_u(k)%x < xy_local(1) .and. k <= Np)
                             k = k + 1
                         end do
                         if (f_approx_equal(airfoil_grid_u(k)%x, xy_local(1))) then
-                            if (xy_local(1) <= airfoil_grid_u(k)%y) then
+                            if (xy_local(2) <= airfoil_grid_u(k)%y) then
                                 !!IB
                                 ib_markers_sf(i, j, 0) = patch_id
                             end if
                         else
                             f = (airfoil_grid_u(k)%x - xy_local(1))/(airfoil_grid_u(k)%x - airfoil_grid_u(k - 1)%x)
-                            if (xy_local(1) <= ((1._wp - f)*airfoil_grid_u(k)%y + f*airfoil_grid_u(k - 1)%y)) then
+                            if (xy_local(2) <= ((1._wp - f)*airfoil_grid_u(k)%y + f*airfoil_grid_u(k - 1)%y)) then
                                 !!IB
                                 ib_markers_sf(i, j, 0) = patch_id
                             end if
@@ -297,14 +299,14 @@ contains
                             k = k + 1
                         end do
                         if (f_approx_equal(airfoil_grid_l(k)%x, xy_local(1))) then
-                            if (xy_local(1) >= airfoil_grid_l(k)%y) then
+                            if (xy_local(2) >= airfoil_grid_l(k)%y) then
                                 !!IB
                                 ib_markers_sf(i, j, 0) = patch_id
                             end if
                         else
                             f = (airfoil_grid_l(k)%x - xy_local(1))/(airfoil_grid_l(k)%x - airfoil_grid_l(k - 1)%x)
 
-                            if (xy_local(1) >= ((1._wp - f)*airfoil_grid_l(k)%y + f*airfoil_grid_l(k - 1)%y)) then
+                            if (xy_local(2) >= ((1._wp - f)*airfoil_grid_l(k)%y + f*airfoil_grid_l(k - 1)%y)) then
                                 !!IB
                                 ib_markers_sf(i, j, 0) = patch_id
                             end if

src/simulation/m_ibm.fpp

@@ -396,6 +396,7 @@ contains
         integer :: index
 
         do q = 1, num_gps
+          print *, q, num_gps
             gp = ghost_points_in(q)
             i = gp%loc(1)
             j = gp%loc(2)
@@ -445,7 +446,6 @@ contains
                                .or. temp_loc > s_cc(index + 1)))
                         index = index + dir
                         if (index < -buff_size .or. index > bound) then
-                            print *, q, index, bound, buff_size
                             print *, "temp_loc=", temp_loc, " s_cc(index)=", s_cc(index), " s_cc(index+1)=", s_cc(index + 1)
                             print *, "Increase buff_size further in m_helper_basic (currently set to a minimum of 10)"
                             error stop "Increase buff_size"

src/simulation/m_start_up.fpp

@@ -581,7 +581,6 @@ contains
                 if (patch_ib(i)%c > 0) then
                     Np = int((patch_ib(i)%p*patch_ib(i)%c/dx(0))*20) + int(((patch_ib(i)%c - patch_ib(i)%p*patch_ib(i)%c)/dx(0))*20) + 1
                     allocate (MPI_IO_airfoil_IB_DATA%var(1:2*Np))
-                    print *, "HERE Np", Np
                 end if
             end do
         end if
@@ -943,8 +942,6 @@ contains
             do j = 1, num_ibs
                 if (patch_ib(j)%c > 0) then
 
-                    print *, "HERE Np", Np
-
                     allocate (airfoil_grid_u(1:Np))
                     allocate (airfoil_grid_l(1:Np))
 

src/simulation/p_main.fpp

@@ -53,7 +53,6 @@ program p_main
     call nvtxStartRange("INIT-GPU-VARS")
     call s_initialize_gpu_vars()
     call nvtxEndRange
-    print *, patch_ib(1)%moving_ibm, patch_ib(1)%vel(1), patch_ib(1)%vel(2), patch_ib(1)%vel(3)
 
     ! Setting the time-step iterator to the first time-step
     if (cfl_dt) then