Spelling and formatting run

Author: danieljvickers

src/common/m_derived_types.fpp

@@ -329,9 +329,9 @@ module m_derived_types
         !! Patch conditions for moving imersed boundaries
         integer :: moving_ibm ! 0 for no moving, 1 for moving, 2 for moving on forced path
         real(wp), dimension(1:3) :: vel
-        real(wp), dimension(1:3) :: step_vel ! velcoity array used to store intermediate steps in the time_stepper module
+        real(wp), dimension(1:3) :: step_vel ! velocity array used to store intermediate steps in the time_stepper module
         real(wp), dimension(1:3) :: angular_vel
-        real(wp), dimension(1:3) :: step_angular_vel ! velcoity array used to store intermediate steps in the time_stepper module
+        real(wp), dimension(1:3) :: step_angular_vel ! velocity array used to store intermediate steps in the time_stepper module
 
     end type ib_patch_parameters
 

src/common/m_ib_patches.fpp

@@ -529,9 +529,9 @@ contains
         ! Computing the beginning and the end x- and y-coordinates of the
         ! rectangle based on its centroid and lengths
         x_boundary%beg = -0.5_wp*length_x
-        x_boundary%end =  0.5_wp*length_x
+        x_boundary%end = 0.5_wp*length_x
         y_boundary%beg = -0.5_wp*length_y
-        y_boundary%end =  0.5_wp*length_y
+        y_boundary%end = 0.5_wp*length_y
 
         ! Since the rectangular patch does not allow for its boundaries to
         ! be smoothed out, the pseudo volume fraction is set to 1 to ensure
@@ -545,7 +545,7 @@ contains
         ! variables of the current patch are assigned to this cell.
         do j = 0, n
             do i = 0, m
-                ! get the x and y coodinates in the local IB frame
+                ! get the x and y coordinates in the local IB frame
                 xy_local = [x_cc(i) - x_centroid, y_cc(j) - y_centroid, 0._wp]
                 xy_local = matmul(inverse_rotation, xy_local)
                 if (x_boundary%beg <= xy_local(1) .and. &
@@ -1004,36 +1004,36 @@ contains
 
         ! construct the x, y, and z rotation matrices
         if (num_dims == 3) then
-          ! also compute the x and y axes in 3D
-          angle = patch_ib(patch_id)%angles(1)
-          rotation(1, 1, :) = [1._wp, 0._wp     , 0._wp      ]
-          rotation(1, 2, :) = [0._wp, cos(angle), -sin(angle)]
-          rotation(1, 3, :) = [0._wp, sin(angle), cos(angle) ]
-
-          angle = patch_ib(patch_id)%angles(2)
-          rotation(2, 1, :) = [cos(angle) , 0._wp, sin(angle)]
-          rotation(2, 2, :) = [0._wp      , 1._wp, 0._wp     ]
-          rotation(2, 3, :) = [-sin(angle), 0._wp, cos(angle)]
-
-          ! apply the y rotation to the x rotation
-          patch_ib(patch_id)%rotation_matrix(:, :) = matmul(rotation(1, :, :), rotation(2, :, :))
-          patch_ib(patch_id)%rotation_matrix_inverse(:, :) = matmul(transpose(rotation(2, :, :)), transpose(rotation(1, :, :)))
+            ! also compute the x and y axes in 3D
+            angle = patch_ib(patch_id)%angles(1)
+            rotation(1, 1, :) = [1._wp, 0._wp, 0._wp]
+            rotation(1, 2, :) = [0._wp, cos(angle), -sin(angle)]
+            rotation(1, 3, :) = [0._wp, sin(angle), cos(angle)]
+
+            angle = patch_ib(patch_id)%angles(2)
+            rotation(2, 1, :) = [cos(angle), 0._wp, sin(angle)]
+            rotation(2, 2, :) = [0._wp, 1._wp, 0._wp]
+            rotation(2, 3, :) = [-sin(angle), 0._wp, cos(angle)]
+
+            ! apply the y rotation to the x rotation
+            patch_ib(patch_id)%rotation_matrix(:, :) = matmul(rotation(1, :, :), rotation(2, :, :))
+            patch_ib(patch_id)%rotation_matrix_inverse(:, :) = matmul(transpose(rotation(2, :, :)), transpose(rotation(1, :, :)))
         end if
 
         ! z component first, since it applies in 2D and 3D
         angle = patch_ib(patch_id)%angles(3)
         rotation(3, 1, :) = [cos(angle), -sin(angle), 0._wp]
-        rotation(3, 2, :) = [sin(angle), cos(angle) , 0._wp]
-        rotation(3, 3, :) = [0._wp     , 0._wp      , 1._wp]
+        rotation(3, 2, :) = [sin(angle), cos(angle), 0._wp]
+        rotation(3, 3, :) = [0._wp, 0._wp, 1._wp]
 
         if (num_dims == 3) then
-          ! apply the z rotation to the xy rotation in 3D
-          patch_ib(patch_id)%rotation_matrix(:, :) = matmul(patch_ib(patch_id)%rotation_matrix(:, :), rotation(3, :, :))
-          patch_ib(patch_id)%rotation_matrix_inverse(:, :) = matmul(transpose(rotation(3, :, :)), patch_ib(patch_id)%rotation_matrix_inverse(:, :))
+            ! apply the z rotation to the xy rotation in 3D
+            patch_ib(patch_id)%rotation_matrix(:, :) = matmul(patch_ib(patch_id)%rotation_matrix(:, :), rotation(3, :, :))
+            patch_ib(patch_id)%rotation_matrix_inverse(:, :) = matmul(transpose(rotation(3, :, :)), patch_ib(patch_id)%rotation_matrix_inverse(:, :))
         else
-          ! write out only the z rotation in 2D
-          patch_ib(patch_id)%rotation_matrix(:, :) = rotation(3, :, :)
-          patch_ib(patch_id)%rotation_matrix_inverse(:, :) = transpose(rotation(3, :, :))
+            ! write out only the z rotation in 2D
+            patch_ib(patch_id)%rotation_matrix(:, :) = rotation(3, :, :)
+            patch_ib(patch_id)%rotation_matrix_inverse(:, :) = transpose(rotation(3, :, :))
         end if
 
     end subroutine s_update_ib_rotation_matrix

src/pre_process/m_global_parameters.fpp

@@ -545,7 +545,7 @@ contains
             patch_ib(i)%angles(:) = 0._wp
             patch_ib(i)%angular_vel(:) = 0._wp
 
-            ! sets values of a rotation matrix which can be used when calculating rotaitons
+            ! sets values of a rotation matrix which can be used when calculating rotations
             patch_ib(i)%rotation_matrix = 0._wp
             patch_ib(i)%rotation_matrix(1, 1) = 1._wp
             patch_ib(i)%rotation_matrix(2, 2) = 1._wp

src/pre_process/m_initial_condition.fpp

@@ -193,7 +193,7 @@ contains
 
         if (ib) then
             do i = 1, num_ibs
-              ! call s_update_ib_rotation_matrix(i)
+                ! call s_update_ib_rotation_matrix(i)
             end do
             call s_apply_ib_patches(ib_markers%sf, levelset, levelset_norm)
         end if

src/simulation/m_ibm.fpp

@@ -94,7 +94,7 @@ contains
         do i = 1, num_ibs
             if (patch_ib(i)%moving_ibm /= 0) then
                 moving_immersed_boundary_flag = .true.
-                
+
             end if
             call s_update_ib_rotation_matrix(i)
         end do
@@ -270,8 +270,8 @@ contains
                 else
                     ! get the vector that points from the centroid to the ghost
                     radial_vector = physical_loc - [patch_ib(patch_id)%x_centroid, &
-                        patch_ib(patch_id)%y_centroid, patch_ib(patch_id)%z_centroid]
-                    ! convert the angular velcoity from the inertial reference frame to the fluids frame, then convert to linear velocity
+                                                    patch_ib(patch_id)%y_centroid, patch_ib(patch_id)%z_centroid]
+                    ! convert the angular velocity from the inertial reference frame to the fluids frame, then convert to linear velocity
                     rotation_velocity = cross_product(matmul(patch_ib(patch_id)%rotation_matrix, patch_ib(patch_id)%angular_vel), radial_vector)
                     do q = 1, 3
                         ! if mibm is 1 or 2, then the boundary may be moving
@@ -920,7 +920,7 @@ contains
 
         ! recalulcate the rotation matrix based upon the new angles
         do i = 1, num_ibs
-          if (patch_ib(i)%moving_ibm .ne. 0) call s_update_ib_rotation_matrix(i)
+            if (patch_ib(i)%moving_ibm /= 0) call s_update_ib_rotation_matrix(i)
         end do
 
         ! recompute the new ib_patch locations and broadcast them.
@@ -952,13 +952,13 @@ contains
     end subroutine s_finalize_ibm_module
 
     function cross_product(a, b) result(c)
-      implicit none
-      real(8), intent(in) :: a(3), b(3)
-      real(8) :: c(3)
+        implicit none
+        real(8), intent(in) :: a(3), b(3)
+        real(8) :: c(3)
 
-      c(1) = a(2)*b(3) - a(3)*b(2)
-      c(2) = a(3)*b(1) - a(1)*b(3)
-      c(3) = a(1)*b(2) - a(2)*b(1)
+        c(1) = a(2)*b(3) - a(3)*b(2)
+        c(2) = a(3)*b(1) - a(1)*b(3)
+        c(3) = a(1)*b(2) - a(2)*b(1)
     end function cross_product
 
 end module m_ibm