Ran formatting

Author: danieljvickers

src/common/m_ib_patches.fpp

@@ -139,7 +139,7 @@ contains
 
         integer, intent(in) :: patch_id
         integer, dimension(0:m, 0:n, 0:p), intent(inout) :: ib_markers_sf
-        
+
         real(wp), dimension(1:2) :: center
         real(wp) :: radius
 
@@ -540,9 +540,9 @@ contains
                 xy_local = matmul(inverse_rotation, xy_local)
 
                 if (-0.5_wp*length(1) <= xy_local(1) .and. &
-                    0.5_wp*length(1)  >= xy_local(1) .and. &
+                    0.5_wp*length(1) >= xy_local(1) .and. &
                     -0.5_wp*length(2) <= xy_local(2) .and. &
-                    0.5_wp*length(2)  >= xy_local(2)) then
+                    0.5_wp*length(2) >= xy_local(2)) then
 
                     ! Updating the patch identities bookkeeping variable
                     ib_markers_sf(i, j, 0) = patch_id
@@ -569,7 +569,7 @@ contains
         ! Generic loop iterators
         integer :: i, j, k
         real(wp) :: radius
-        real(wp), dimension(1:3):: center
+        real(wp), dimension(1:3) :: center
 
         !! Variables to initialize the pressure field that corresponds to the
             !! bubble-collapse test case found in Tiwari et al. (2013)

src/post_process/m_data_input.f90

@@ -158,7 +158,7 @@ impure subroutine s_read_ib_data_files(file_loc_base, t_step)
         integer, dimension(MPI_STATUS_SIZE) :: status
         integer(KIND=MPI_OFFSET_KIND) :: disp
         integer :: m_MOK, n_MOK, p_MOK, MOK, WP_MOK, save_index
-        
+
 #endif
         if (.not. ib) return
 
@@ -179,7 +179,7 @@ impure subroutine s_read_ib_data_files(file_loc_base, t_step)
                 p_MOK = int(p_glb + 1, MPI_OFFSET_KIND)
                 MOK = int(1._wp, MPI_OFFSET_KIND)
                 WP_MOK = int(8._wp, MPI_OFFSET_KIND)
-                save_index = t_step / t_step_save ! get the number of saves done to this point
+                save_index = t_step/t_step_save ! get the number of saves done to this point
 
                 data_size = (m + 1)*(n + 1)*(p + 1)
                 var_MOK = int(sys_size + 1, MPI_OFFSET_KIND)

src/simulation/m_data_output.fpp

@@ -1022,16 +1022,16 @@ contains
 
             !Writeib data
             if (ib) then
-              write (file_loc, '(A)') 'ib.dat'
-              file_loc = trim(case_dir)//'/restart_data'//trim(mpiiofs)//trim(file_loc)
-              call MPI_FILE_OPEN(MPI_COMM_WORLD, file_loc, ior(MPI_MODE_WRONLY, MPI_MODE_CREATE), &
-                               mpi_info_int, ifile, ierr)
-                               
+                write (file_loc, '(A)') 'ib.dat'
+                file_loc = trim(case_dir)//'/restart_data'//trim(mpiiofs)//trim(file_loc)
+                call MPI_FILE_OPEN(MPI_COMM_WORLD, file_loc, ior(MPI_MODE_WRONLY, MPI_MODE_CREATE), &
+                                   mpi_info_int, ifile, ierr)
+
                 var_MOK = int(sys_size + 1, MPI_OFFSET_KIND)
-                disp = m_MOK*max(MOK, n_MOK)*max(MOK, p_MOK)*WP_MOK*(var_MOK - 1 + int(t_step / t_step_save))
+                disp = m_MOK*max(MOK, n_MOK)*max(MOK, p_MOK)*WP_MOK*(var_MOK - 1 + int(t_step/t_step_save))
 
                 call MPI_FILE_SET_VIEW(ifile, disp, MPI_INTEGER, MPI_IO_IB_DATA%view, &
-                                        'native', mpi_info_int, ierr)
+                                       'native', mpi_info_int, ierr)
                 call MPI_FILE_WRITE_ALL(ifile, MPI_IO_IB_DATA%var%sf, data_size, &
                                         MPI_INTEGER, status, ierr)
                 call MPI_FILE_CLOSE(ifile, ierr)

src/simulation/m_ibm.fpp

@@ -271,9 +271,9 @@ contains
                     radial_vector = physical_loc - [patch_ib(patch_id)%x_centroid, &
                                                     patch_ib(patch_id)%y_centroid, patch_ib(patch_id)%z_centroid]
                     rotation_velocity = cross_product(matmul(patch_ib(patch_id)%rotation_matrix, patch_ib(patch_id)%angular_vel), radial_vector)
-                    
+
                     ! add only the component of the IB's motion that is normal to the surface
-                    vel_g = vel_g + sum((patch_ib(patch_id)%vel+rotation_velocity)*norm)*norm
+                    vel_g = vel_g + sum((patch_ib(patch_id)%vel + rotation_velocity)*norm)*norm
                 end if
             else
                 if (patch_ib(patch_id)%moving_ibm == 0) then