Spelling and formatting

Author: danieljvickers

docs/documentation/case.md

@@ -309,7 +309,7 @@ This is enabled by adding `'elliptic_smoothing': "T",` and `'elliptic_smoothing_
 | `model_threshold`    | Real    | Ray fraction inside the model patch above which the fraction is set to one.|
 | `moving_ibm`         | Integer | Sets the method used for IB movement. |
 | `vel(i)`             | Real    | Initial velocity of the moving IB in the i-th direction. |
-| `angular_vel(i)`     | Real    | Initial angular velcoity of the moving IB in the i-th direction. |
+| `angular_vel(i)`     | Real    | Initial angular velocity of the moving IB in the i-th direction. |
 
 These parameters should be prepended with `patch_ib(j)%` where $j$ is the patch index.
 
@@ -336,9 +336,9 @@ Additional details on this specification can be found in [The Naca Airfoil Serie
 
 - `moving_ibm` sets the method by which movement will be applied to the immersed boundary. Using 0 will result in no movement. Using 1 will result 1-way coupling where the boundary moves at a constant rate and applied forces to the fluid based upon it's own motion. In 1-way coupling, the fluid does not apply forces back onto the IB.
 
-- `vel(i)` is the initial linear velocity of the IB in the x, y, z direction for i=1, 2, 3. When `moving_ibm` equals 1, this velcoity is constant.
+- `vel(i)` is the initial linear velocity of the IB in the x, y, z direction for i=1, 2, 3. When `moving_ibm` equals 1, this velocity is constant.
 
-- `angular_vel(i)` is the initial angular velocity of the IB about the x, y, z axes for i=1, 2, 3 in radians per second. When `moving_ibm` equals 1, this angular velcoity is constant.
+- `angular_vel(i)` is the initial angular velocity of the IB about the x, y, z axes for i=1, 2, 3 in radians per second. When `moving_ibm` equals 1, this angular velocity is constant.
 
 ### 5. Fluid Material’s
 

src/common/m_mpi_common.fpp

@@ -486,10 +486,10 @@ contains
         ! Performing the reduction procedure
         if (loc(var_loc) == loc(var_glb)) then
             call MPI_Allreduce(MPI_IN_PLACE, var_glb, num_vectors*vector_length, &
-                              mpi_p, MPI_SUM, MPI_COMM_WORLD, ierr)
+                               mpi_p, MPI_SUM, MPI_COMM_WORLD, ierr)
         else
             call MPI_Allreduce(var_loc, var_glb, num_vectors*vector_length, &
-                              mpi_p, MPI_SUM, MPI_COMM_WORLD, ierr)
+                               mpi_p, MPI_SUM, MPI_COMM_WORLD, ierr)
         end if
 
 #endif

src/simulation/m_ibm.fpp

@@ -202,7 +202,7 @@ contains
                 do j = 0, m
                     do k = 0, n
                         do l = 0, p
-                            if (ib_markers%sf(j,k,l) == patch_id) then
+                            if (ib_markers%sf(j, k, l) == patch_id) then
                                 q_prim_vf(E_idx)%sf(j, k, l) = 1._wp
                             end if
                         end do
@@ -264,8 +264,8 @@ contains
                     if (patch_ib(patch_id)%moving_ibm == 0) then
                         q_prim_vf(E_idx)%sf(j, k, l) = pres_IP
                     else
-                      ! TODO :: improve for two fluid
-                      q_prim_vf(E_idx)%sf(j, k, l) = pres_IP / (1._wp - 2._wp*abs(levelset%sf(j,k,l,patch_id)*alpha_rho_IP(q)/pres_IP) * dot_product(patch_ib(patch_id)%force/patch_ib(patch_id)%mass, levelset_norm%sf(j,k,l,patch_id,:)))
+                        ! TODO :: improve for two fluid
+                        q_prim_vf(E_idx)%sf(j, k, l) = pres_IP/(1._wp - 2._wp*abs(levelset%sf(j, k, l, patch_id)*alpha_rho_IP(q)/pres_IP)*dot_product(patch_ib(patch_id)%force/patch_ib(patch_id)%mass, levelset_norm%sf(j, k, l, patch_id, :)))
                     end if
                 end do
 
@@ -1019,21 +1019,21 @@ contains
                         if (patch_ib(ib_idx)%moving_ibm == 2) then
                             if (num_dims == 3) then
                                 radial_vector = [x_cc(i), y_cc(j), z_cc(k)] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, patch_ib(ib_idx)%z_centroid] ! get the vector pointing to the grid cell
-                            else 
+                            else
                                 radial_vector = [x_cc(i), y_cc(j), 0._wp] - [patch_ib(ib_idx)%x_centroid, patch_ib(ib_idx)%y_centroid, 0._wp] ! get the vector pointing to the grid cell
                             end if
-                            pressure_divergence(1) = (pressure(i+1, j, k) - pressure(i-1, j, k)) / (2._wp * x_cc(i))
-                            pressure_divergence(2) = (pressure(i, j+1, k) - pressure(i, j-1, k)) / (2._wp * y_cc(j))
-                            cell_volume = x_cc(i) * y_cc(j)
+                            pressure_divergence(1) = (pressure(i + 1, j, k) - pressure(i - 1, j, k))/(2._wp*x_cc(i))
+                            pressure_divergence(2) = (pressure(i, j + 1, k) - pressure(i, j - 1, k))/(2._wp*y_cc(j))
+                            cell_volume = x_cc(i)*y_cc(j)
                             if (num_dims == 3) then
-                                pressure_divergence(3) = (pressure(i, j, k+1) - pressure(i, j, k-1)) / (2._wp * z_cc(k))
-                                cell_volume = cell_volume * z_cc(k)
+                                pressure_divergence(3) = (pressure(i, j, k + 1) - pressure(i, j, k - 1))/(2._wp*z_cc(k))
+                                cell_volume = cell_volume*z_cc(k)
                             else
                                 pressure_divergence(3) = 0._wp
                             end if
 
-                            forces(ib_idx, :) = forces(ib_idx, :) - (pressure_divergence * cell_volume)
-                            torques(ib_idx, :) = torques(ib_idx, :) + (cross_product(radial_vector, pressure_divergence) * cell_volume)
+                            forces(ib_idx, :) = forces(ib_idx, :) - (pressure_divergence*cell_volume)
+                            torques(ib_idx, :) = torques(ib_idx, :) + (cross_product(radial_vector, pressure_divergence)*cell_volume)
                         end if
                     end if
                 end do
@@ -1068,9 +1068,9 @@ contains
 
         real(wp), dimension(3), optional :: axis !< the axis about which we compute the moment. Only required in 3D.
         integer, intent(in) :: ib_marker
-        
-        real (wp) :: moment, distance_to_axis, cell_volume
-        real (wp), dimension(3) :: position, closest_point_along_axis, vector_to_axis
+
+        real(wp) :: moment, distance_to_axis, cell_volume
+        real(wp), dimension(3) :: position, closest_point_along_axis, vector_to_axis
         integer :: i, j, k, count
 
         if (p == 0) then
@@ -1080,26 +1080,25 @@ contains
             patch_ib(ib_marker)%moment = 1._wp
             return
         else
-            axis = axis / sqrt(sum(axis))
+            axis = axis/sqrt(sum(axis))
         end if
-        
+
         ! if the IB is in 2D or a 3D sphere, we can compute this exactly
         if (patch_ib(ib_marker)%geometry == 2) then ! circle
-            patch_ib(ib_marker)%moment = 0.5 * patch_ib(ib_marker)%mass * (patch_ib(ib_marker)%radius)**2
+            patch_ib(ib_marker)%moment = 0.5*patch_ib(ib_marker)%mass*(patch_ib(ib_marker)%radius)**2
         elseif (patch_ib(ib_marker)%geometry == 3) then ! rectangle
-            patch_ib(ib_marker)%moment = patch_ib(i)%mass * (patch_ib(ib_marker)%length_x**2 + patch_ib(ib_marker)%length_y**2) / 6._wp                
+            patch_ib(ib_marker)%moment = patch_ib(i)%mass*(patch_ib(ib_marker)%length_x**2 + patch_ib(ib_marker)%length_y**2)/6._wp
         elseif (patch_ib(ib_marker)%geometry == 8) then ! sphere
-            patch_ib(ib_marker)%moment = 0.4 * patch_ib(ib_marker)%mass * (patch_ib(ib_marker)%radius)**2
-          
+            patch_ib(ib_marker)%moment = 0.4*patch_ib(ib_marker)%mass*(patch_ib(ib_marker)%radius)**2
+
         else ! we do not have an analytic moment of inertia calculation and need to approximate it directly
             count = 0
             moment = 0._wp
             cell_volume = (x_cc(1) - x_cc(0))*(y_cc(1) - y_cc(0)) ! computed without grid stretching. Update in the loop to perform with stretching
             if (p /= 0) then
-              cell_volume = cell_volume * (z_cc(1) - z_cc(0))
+                cell_volume = cell_volume*(z_cc(1) - z_cc(0))
             end if
 
-            
             $:GPU_PARALLEL_LOOP(private='[position,closest_point_along_axis,vector_to_axis,distance_to_axis]', copy='[moment,count]', copyin='[ib_marker,cell_volume,axis]', collapse=3)
             do i = 0, m
                 do j = 0, j
@@ -1130,7 +1129,7 @@ contains
             $:END_GPU_PARALLEL_LOOP()
 
             ! write the final moment assuming the points are all uniform density
-            patch_ib(ib_marker)%moment = moment * patch_ib(ib_marker)%mass / (count * cell_volume)
+            patch_ib(ib_marker)%moment = moment*patch_ib(ib_marker)%mass/(count*cell_volume)
             $:GPU_UPDATE(device='[patch_ib(ib_marker)%moment]')
         end if
 

src/simulation/m_time_steppers.fpp

@@ -632,9 +632,9 @@ contains
                                 patch_ib(i)%vel = patch_ib(i)%vel + rk_coef(s, 3)*dt*(patch_ib(i)%force/patch_ib(i)%mass)/rk_coef(s, 4)
 
                                 ! update the angular velocity with the torque value
-                                patch_ib(i)%angular_vel = (patch_ib(i)%angular_vel * patch_ib(i)%moment) + (rk_coef(s, 3)*dt*patch_ib(i)%torque/rk_coef(s, 4)) ! add the torque to the angular momentum
+                                patch_ib(i)%angular_vel = (patch_ib(i)%angular_vel*patch_ib(i)%moment) + (rk_coef(s, 3)*dt*patch_ib(i)%torque/rk_coef(s, 4)) ! add the torque to the angular momentum
                                 call s_compute_moment_of_inertia(i, patch_ib(i)%angular_vel) ! update the moment of inertia to be based on the direction of the angular momentum
-                                patch_ib(i)%angular_vel = patch_ib(i)%angular_vel / patch_ib(i)%moment ! convert back to angular velocity with the new moment of inertia
+                                patch_ib(i)%angular_vel = patch_ib(i)%angular_vel/patch_ib(i)%moment ! convert back to angular velocity with the new moment of inertia
                             end if
 
                             ! Update the angle of the IB