code cleaning

Author: danielpeter

src/generate_databases/model_gll.f90

@@ -180,7 +180,7 @@ subroutine model_gll(myrank,nspec,LOCAL_PATH)
   rho_vs(:,:,:,:) = rhostore(:,:,:,:) * vs_read(:,:,:,:)
 
   if (ANISOTROPY) then
-    
+
     ! c11
     if (myrank == 0) write(IMAIN,*) '     reading in: c11.bin'
 
@@ -284,7 +284,7 @@ subroutine model_gll(myrank,nspec,LOCAL_PATH)
 
     read(IIN) c23store
     close(IIN)
-    
+
     ! c24
     if (myrank == 0) write(IMAIN,*) '     reading in: c24.bin'
 
@@ -349,7 +349,7 @@ subroutine model_gll(myrank,nspec,LOCAL_PATH)
 
     read(IIN) c34store
     close(IIN)
-    
+
     ! c35
     if (myrank == 0) write(IMAIN,*) '     reading in: c35.bin'
 

src/gpu/compute_stacey_acoustic_cuda.cu

@@ -282,7 +282,7 @@ void FC_FUNC_(compute_stacey_acoustic_undoatt_cuda,
   GPU_ERROR_CHECKING("compute_stacey_acoustic_undoatt_cuda");
 }
 
-extern EXTERN_LANG 
+extern EXTERN_LANG
 void FC_FUNC_(compute_coupled_injection_contribution_ac_device,
               COMPUTE_COUPLED_INJECTION_CONTRIBUTION_AC_DEVICE) (long *Mesh_pointer,
                                                       realw* b_boundary_injection_potential,
@@ -317,7 +317,7 @@ void FC_FUNC_(compute_coupled_injection_contribution_ac_device,
         SAVE_STACEY,mp->d_num_abs_boundary_faces,mp->d_b_boundary_injection_potential
       );
     }
-#endif 
+#endif
 #ifdef USE_HIP
     if (run_hip){
       hipLaunchKernelGGL(compute_stacey_acoustic_injection_kernel,dim3(grid), dim3(threads), 0, mp->compute_stream,

src/gpu/kernels/compute_stacey_acoustic_kernel.cu

@@ -321,11 +321,11 @@ __global__ void compute_stacey_acoustic_injection_kernel(const realw* veloc_inj,
       nx = abs_boundary_normal[INDEX3(NDIM,NGLL2,0,igll,iface)];
       ny = abs_boundary_normal[INDEX3(NDIM,NGLL2,1,igll,iface)];
       nz = abs_boundary_normal[INDEX3(NDIM,NGLL2,2,igll,iface)];
-      
-      // n \cdot u 
+
+      // n \cdot u
       un = ux * nx + uy * ny + uz * nz;
       un += chi_dot /(cpl * rhol);
-    
+
       // gets associated, weighted jacobian
       jacobianw = abs_boundary_jacobian2Dw[INDEX2(NGLL2,igll,iface)];
 

src/gpu/kernels/kernel_proto.cu.h

@@ -911,7 +911,7 @@ __global__ void compute_stacey_elastic_undoatt_kernel(realw* veloc,
                                                       int num_abs_boundary_faces) ;
 
 // nqdu added
-__global__ void 
+__global__ void
 compute_stacey_elastic_injection_kernel(const realw* veloc_inj,
                                         const realw* tract_inj,
                                         realw* accel,
@@ -944,7 +944,7 @@ __global__ void compute_stacey_acoustic_injection_kernel(const realw* veloc_inj,
                                        int SAVE_STACEY,
                                        int num_abs_boundary_faces,
                                        field* b_boundary_injection_potential
-                                       );                                       
+                                       );
 
 //
 // src/gpu/kernels/enforce_free_surface_cuda_kernel.cu

src/specfem3D/compute_stacey_acoustic.f90

@@ -566,7 +566,7 @@ subroutine compute_coupled_injection_contribution_ac(NGLOB_AB,potential_dot_dot_
           ! point index using table lookup
           ipt = ipt_table(igll,iface)
 
-          ! interpolates velocity/stress  
+          ! interpolates velocity/stress
           ! nqdu NOT CORRECT!!!
           ! vx_FK = cs1 * Veloc_FK(1,ipt,iim1) + cs2 * Veloc_FK(1,ipt,ii) + cs3 * Veloc_FK(1,ipt,iip1) + cs4 * Veloc_FK(1,ipt,iip2)
           ! vy_FK = cs1 * Veloc_FK(2,ipt,iim1) + cs2 * Veloc_FK(2,ipt,ii) + cs3 * Veloc_FK(2,ipt,iip1) + cs4 * Veloc_FK(2,ipt,iip2)
@@ -575,7 +575,7 @@ subroutine compute_coupled_injection_contribution_ac(NGLOB_AB,potential_dot_dot_
           ! tx_FK = cs1 * Tract_FK(1,ipt,iim1) + cs2 * Tract_FK(1,ipt,ii) + cs3 * Tract_FK(1,ipt,iip1) + cs4 * Tract_FK(1,ipt,iip2)
           ! ty_FK = cs1 * Tract_FK(2,ipt,iim1) + cs2 * Tract_FK(2,ipt,ii) + cs3 * Tract_FK(2,ipt,iip1) + cs4 * Tract_FK(2,ipt,iip2)
           ! tz_FK = cs1 * Tract_FK(3,ipt,iim1) + cs2 * Tract_FK(3,ipt,ii) + cs3 * Tract_FK(3,ipt,iip1) + cs4 * Tract_FK(3,ipt,iip2)
-          
+
           ! now displ/chi_dot are saved in veloc_fk/tract_fk
           ux_FK = cs1 * Veloc_FK(1,ipt,iim1) + cs2 * Veloc_FK(1,ipt,ii) + &
                   cs3 * Veloc_FK(1,ipt,iip1) + cs4 * Veloc_FK(1,ipt,iip2)
@@ -585,7 +585,7 @@ subroutine compute_coupled_injection_contribution_ac(NGLOB_AB,potential_dot_dot_
                   cs3 * Veloc_FK(3,ipt,iip1) + cs4 * Veloc_FK(3,ipt,iip2)
           chi_dot_FK = cs1 * Tract_FK(1,ipt,iim1) + cs2 * Tract_FK(1,ipt,ii) + &
                        cs3 * Tract_FK(1,ipt,iip1) + cs4 * Tract_FK(1,ipt,iip2)
-          
+
           ! velocity
           ! vx = - vx_FK
           ! vy = - vy_FK
@@ -642,7 +642,7 @@ subroutine compute_coupled_injection_contribution_ac(NGLOB_AB,potential_dot_dot_
         !       here, it is added to first undo the factor (rhol*cpl) used to add the velocity contribution to the traction together
         !       with another (rhol*cpl) used in the expressions of the Sommerfeld condition, where
         !          absorbl = potential_dot_acoustic(iglob) * jacobianw / cpl / rhol
-        
+
         un = ux*nx + uy*ny + uz*nz
         un = un + chi_dot / (rhol * cpl)
 
@@ -673,25 +673,25 @@ end subroutine compute_coupled_injection_contribution_ac
   subroutine compute_coupled_injection_contribution_ac_GPU(iphase,Mesh_pointer)
 
     use constants
-  
+
     use specfem_par, only: SAVE_STACEY,SIMULATION_TYPE
-  
-    use specfem_par, only:  num_abs_boundary_faces
-  
+
+    use specfem_par, only: num_abs_boundary_faces
+
     ! boundary coupling
     use shared_parameters, only: COUPLE_WITH_INJECTION_TECHNIQUE
     ! boundary injection wavefield parts for saving together with b_absorb_field
     use specfem_par_coupling, only: b_boundary_injection_potential
-  
+
     implicit none
-  
+
     ! communication overlap
     integer,intent(in) :: iphase
-  
+
     ! GPU_MODE variables
     integer(kind=8),intent(in) :: Mesh_pointer
 
-  
+
   ! safety checks
   if (.not. COUPLE_WITH_INJECTION_TECHNIQUE) return
 
@@ -703,13 +703,13 @@ subroutine compute_coupled_injection_contribution_ac_GPU(iphase,Mesh_pointer)
 
   ! only for forward wavefield
   if (SIMULATION_TYPE /= 1) return
-  
-  
+
+
     ! compute contribution in device
-    call compute_coupled_injection_contribution_ac_device(Mesh_pointer,&
-                                                          b_boundary_injection_potential,&
+    call compute_coupled_injection_contribution_ac_device(Mesh_pointer, &
+                                                          b_boundary_injection_potential, &
                                                           SAVE_STACEY)
-  
-  
+
+
     end subroutine compute_coupled_injection_contribution_ac_GPU
-  
+

src/specfem3D/compute_stacey_viscoelastic.F90

@@ -616,7 +616,7 @@ subroutine compute_coupled_injection_contribution_el(NGLOB_AB,accel,iphase,it)
 
   ! gets velocity & stress for boundary points
   select case(INJECTION_TECHNIQUE_TYPE)
-  !NQDU don't need 
+  !NQDU don't need
   ! case (INJECTION_TECHNIQUE_IS_DSM)
   !   ! DSM coupling
   !   if (old_DSM_coupling_from_Vadim) then
@@ -799,7 +799,7 @@ subroutine compute_coupled_injection_contribution_el_GPU(iphase,Mesh_pointer)
 
   use specfem_par, only: SAVE_STACEY,SIMULATION_TYPE
 
-  use specfem_par, only:  num_abs_boundary_faces
+  use specfem_par, only: num_abs_boundary_faces
 
   ! boundary coupling
   use shared_parameters, only: COUPLE_WITH_INJECTION_TECHNIQUE
@@ -829,7 +829,7 @@ subroutine compute_coupled_injection_contribution_el_GPU(iphase,Mesh_pointer)
   if (SIMULATION_TYPE /= 1) return
 
   ! compute contribution in device
-  call compute_coupled_injection_contribution_el_device(Mesh_pointer,b_boundary_injection_field,&
+  call compute_coupled_injection_contribution_el_device(Mesh_pointer,b_boundary_injection_field, &
                                                         SAVE_STACEY)