source code cleaning

danielpeter · danielpeter · commit d693a9c94121 · 2025-09-16T11:27:16.000+02:00
diff --git a/src/generate_databases/read_parameters.f90 b/src/generate_databases/read_parameters.f90
@@ -196,16 +196,16 @@ subroutine read_parameters()
     ! else
     !    write(IMAIN,*) 'using a CMTSOLUTION source'
     ! endif
-    if(.not. USE_CMT_AND_FORCE_SOURCE) then 
+    if (.not. USE_CMT_AND_FORCE_SOURCE) then
       if (USE_FORCE_POINT_SOURCE) then
         write(IMAIN,*) 'using a FORCESOLUTION source instead of a CMTSOLUTION source'
       else
         write(IMAIN,*) 'using a CMTSOLUTION source'
       endif
-    else 
+    else
       write(IMAIN,*) 'using both a FORCESOLUTION and a CMTSOLUTION source'
     endif
-    
+
     if (USE_RICKER_TIME_FUNCTION) then
        write(IMAIN,*) '  with a Ricker source time function'
     else
@@ -215,13 +215,13 @@ subroutine read_parameters()
       !  else
       !     write(IMAIN,*) '  with a Gaussian source time function'
       !  endif
-      if(USE_CMT_AND_FORCE_SOURCE) then
+      if (USE_CMT_AND_FORCE_SOURCE) then
        write(IMAIN,*) ' CMTSOLUTION with a quasi-Heaviside source time function'
        write(IMAIN,*) ' FORCESOLUTION with a Gaussian source time function'
-      else 
+      else
         if (USE_FORCE_POINT_SOURCE) then
           write(IMAIN,*) ' with a quasi-Heaviside source time function'
-        else 
+        else
           write(IMAIN,*) ' with a Gaussian source time function'
         endif
       endif
diff --git a/src/gpu/smooth_pde_cuda.cu b/src/gpu/smooth_pde_cuda.cu
@@ -17,7 +17,7 @@ void FC_FUNC_(prepare_smooth_pde_gpu,
   Mesh* mp = (Mesh*)(*Mesh_pointer);
   Smooth_pde_data* sp = (Smooth_pde_data*) malloc(sizeof(Smooth_pde_data));
   *Container_smooth_pde = (long)sp;
-  
+
   gpuCreateCopy_todevice_realw((void**)&sp->d_dat_smooth_glob, dat_smooth_glob, mp->NGLOB_AB);
   gpuMalloc_field((void**)&sp->d_ddat_smooth_glob, mp->NGLOB_AB);
   // initializes values to zero
@@ -70,8 +70,8 @@ void FC_FUNC_(compute_update_element_smooth_pde_gpu,
                                                                mp->xix_regular,mp->jacobian_regular,
                                                                mp->d_hprime_xx,
                                                                mp->d_hprimewgll_xx,
-                                                               mp->d_wgllwgll_xy, 
-                                                               mp->d_wgllwgll_xz, 
+                                                               mp->d_wgllwgll_xy,
+                                                               mp->d_wgllwgll_xz,
                                                                mp->d_wgllwgll_yz,
                                                                sp->cv, sp->ch
                                                                );
@@ -94,8 +94,8 @@ void FC_FUNC_(compute_update_element_smooth_pde_gpu,
                                                                mp->xix_regular,mp->jacobian_regular,
                                                                mp->d_hprime_xx,
                                                                mp->d_hprimewgll_xx,
-                                                               mp->d_wgllwgll_xy, 
-                                                               mp->d_wgllwgll_xz, 
+                                                               mp->d_wgllwgll_xy,
+                                                               mp->d_wgllwgll_xz,
                                                                mp->d_wgllwgll_yz,
                                                                sp->cv, sp->ch
                                                                );
@@ -114,7 +114,7 @@ void FC_FUNC_(transfer_boun_dat_smooth_pde_from_device,
   TRACE("transfer_boun_dat_smooth_pde_from_device");
   Mesh* mp = (Mesh*)(*Mesh_pointer);
   Smooth_pde_data* sp = (Smooth_pde_data*)(*Container_smooth_pde);
-  
+
   if (sp->size_mpi_buffer_smooth > 0) {
     int blocksize = BLOCKSIZE_TRANSFER;
     int size_padded = ((int)ceil(((double)(mp->max_nibool_interfaces_ext_mesh))/((double)blocksize)))*blocksize;
@@ -152,7 +152,7 @@ void FC_FUNC_(transfer_boun_dat_smooth_pde_from_device,
     // copies buffer to CPU
     gpuMemcpy_tohost_field(send_buffer_smooth,sp->d_send_buffer,sp->size_mpi_buffer_smooth);
   }
-  
+
   GPU_ERROR_CHECKING("transfer_boun_dat_smooth_pde_from_device");
 }
 
@@ -162,7 +162,7 @@ void FC_FUNC_(transfer_asmbl_dat_smooth_pde_from_device,
                                                          long * Container_smooth_pde,
                                                          field * recv_buffer_smooth){
   TRACE("transfer_asmbl_dat_smooth_pde_from_device");
-  
+
   Mesh* mp = (Mesh*)(*Mesh_pointer);
   Smooth_pde_data* sp = (Smooth_pde_data*)(*Container_smooth_pde);
 
@@ -244,7 +244,7 @@ void FC_FUNC_(kernel_3_smooth_pde_cuda,
                                                                               sp->d_rvol,
                                                                               size);
   }
-#endif  
+#endif
 
   GPU_ERROR_CHECKING("kernel_3_smooth_pde_cuda");
 }
@@ -295,7 +295,7 @@ void FC_FUNC_(zero_pml_smooth_pde_cuda,
               ZERO_PML_SMOOTH_PDE_CUDA)(long * Mesh_pointer,
                                         long * Container_smooth_pde) {
   TRACE("zero_pml_smooth_pde_cuda");
-  
+
   Mesh* mp = (Mesh*)(*Mesh_pointer);
   Smooth_pde_data* sp = (Smooth_pde_data*)(*Container_smooth_pde);
 
@@ -325,7 +325,7 @@ void FC_FUNC_(zero_pml_smooth_pde_cuda,
                                                                          mp->d_ibool,
                                                                          sp->d_CPML_to_spec);
   }
-#endif 
+#endif
 
   GPU_ERROR_CHECKING("zero_pml_smooth_pde_cuda");
 }
@@ -429,5 +429,5 @@ void FC_FUNC_(transfer_dat_smooth_pde_from_device,
   gpuMemcpy_tohost_realw(dat_smooth_glob, sp->d_dat_smooth_glob, mp->NGLOB_AB);
 
   GPU_ERROR_CHECKING("transfer_dat_smooth_pde_from_device");
-  
+
 }
diff --git a/src/shared/count_number_of_sources.f90 b/src/shared/count_number_of_sources.f90
@@ -197,7 +197,7 @@ subroutine count_number_of_sources_by_type(NSOURCES,sources_filename,source_is_f
 
   ! get the number of lines describing the sources
   !if (USE_FORCE_POINT_SOURCE) then
-  if(source_is_force) then 
+  if (source_is_force) then
     ! FORCESOLUTION
     nlines_per_source = NLINES_PER_FORCESOLUTION_SOURCE
   else
@@ -233,7 +233,7 @@ subroutine count_number_of_sources_by_type(NSOURCES,sources_filename,source_is_f
   ! checks lines are a multiple
   if (mod(icounter,nlines_per_source) /= 0) then
     !if (USE_FORCE_POINT_SOURCE) then
-    if(source_is_force) then 
+    if (source_is_force) then
       print *,'Error: total number of lines in FORCESOLUTION file should be a multiple of ',nlines_per_source
       stop 'Error total number of lines in FORCESOLUTION file should be a multiple of NLINES_PER_FORCESOLUTION_SOURCE'
     else
@@ -261,15 +261,15 @@ end subroutine count_number_of_sources_by_type
 
   ! NQDU for binary source file
   subroutine count_number_of_sources_binary(ns_cmt,ns_force,sources_filename)
-    use constants, only : MAX_STRING_LEN
+    use constants, only: MAX_STRING_LEN
     implicit none
     integer,intent(inout) :: ns_cmt,ns_force
-    character(len=MAX_STRING_LEN),intent(in) :: sources_filename 
+    character(len=MAX_STRING_LEN),intent(in) :: sources_filename
 
     ! local
     integer,parameter :: IO_SOURCE = 12341
 
-    open(IO_SOURCE,file=trim(sources_filename),form='unformatted',&
+    open(IO_SOURCE,file=trim(sources_filename),form='unformatted', &
               action='read',access='stream')
     read(IO_SOURCE) ns_cmt
     read(IO_SOURCE) ns_force
diff --git a/src/shared/read_parameter_file.F90 b/src/shared/read_parameter_file.F90
@@ -426,14 +426,14 @@ subroutine read_parameter_file(BROADCAST_AFTER_READ)
     ! read two flags in CMT + FORCE simulation
     call read_value_logical(USE_CMT_AND_FORCE_SOURCE, 'USE_CMT_AND_FORCE_SOURCE', ier)
     ier = 0
-    if (USE_CMT_AND_FORCE_SOURCE) then 
+    if (USE_CMT_AND_FORCE_SOURCE) then
       ! write(*,'(a)') 'cmt + force simulation is enabled'
       USE_FORCE_POINT_SOURCE  = .true.
     endif
 
     call read_value_logical(USE_BINARY_SOURCE_FILE, 'USE_BINARY_SOURCE_FILE', ier)
     ier = 0
-    if(.not. USE_CMT_AND_FORCE_SOURCE) USE_BINARY_SOURCE_FILE = .false.  ! binary file is disabled 
+    if (.not. USE_CMT_AND_FORCE_SOURCE) USE_BINARY_SOURCE_FILE = .false.  ! binary file is disabled
 
     call read_value_logical(USE_RICKER_TIME_FUNCTION, 'USE_RICKER_TIME_FUNCTION', ier)
     if (ier /= 0) then
@@ -850,8 +850,8 @@ subroutine read_parameter_file(BROADCAST_AFTER_READ)
     endif
 
     ! inverse problem is enabled, disable USE_CMT_AND_FORCE/USE_BINARY_SOURCE_FILE
-    if(INVERSE_FWI_FULL_PROBLEM .and. (USE_CMT_AND_FORCE_SOURCE)) then
-      if(myrank == 0) then
+    if (INVERSE_FWI_FULL_PROBLEM .and. (USE_CMT_AND_FORCE_SOURCE)) then
+      if (myrank == 0) then
         write(*,'(a)') 'INVERSE_FWI_FULL_PROBLEM is enabled, disabling USE_CMT_AND_FORCE_SOURCE'
         write(*,*)
       endif
@@ -1258,8 +1258,8 @@ subroutine read_compute_parameters()
   else
     ! gets number of sources
     !NQDU call count_number_of_sources(NSOURCES,sources_filename)
-    if(USE_CMT_AND_FORCE_SOURCE) then 
-      if(.not. USE_BINARY_SOURCE_FILE) then
+    if (USE_CMT_AND_FORCE_SOURCE) then
+      if (.not. USE_BINARY_SOURCE_FILE) then
         ! Force
         sources_filename = path_to_add(1:len_trim(path_to_add)) // &
                            IN_DATA_FILES(1:len_trim(IN_DATA_FILES))//'FORCESOLUTION'
@@ -1269,17 +1269,17 @@ subroutine read_compute_parameters()
         sources_filename = path_to_add(1:len_trim(path_to_add)) // &
                            IN_DATA_FILES(1:len_trim(IN_DATA_FILES))//'CMTSOLUTION'
         call count_number_of_sources_by_type(NSOURCES_CMT,sources_filename,.false.)
-      else 
+      else
         sources_filename = path_to_add(1:len_trim(path_to_add)) // &
                           IN_DATA_FILES(1:len_trim(IN_DATA_FILES))//'SOLUTION.bin'
         call count_number_of_sources_binary(NSOURCES_CMT,NSOURCES_FORCE,sources_filename)
       endif
       NSOURCES = NSOURCES_CMT + NSOURCES_FORCE
-    else 
+    else
       NSOURCES_CMT = 0
       NSOURCES_FORCE = 0
       call count_number_of_sources_by_type(NSOURCES,sources_filename,USE_FORCE_POINT_SOURCE)
-      if(USE_FORCE_POINT_SOURCE) then 
+      if (USE_FORCE_POINT_SOURCE) then
         NSOURCES_FORCE = NSOURCES
       else
         NSOURCES_CMT = NSOURCES
diff --git a/src/shared/shared_par.F90 b/src/shared/shared_par.F90
@@ -169,7 +169,7 @@ module shared_input_parameters
   logical :: HAS_FINITE_FAULT_SOURCE
 
   ! cmt + point force simulation
-  logical :: USE_CMT_AND_FORCE_SOURCE = .false. 
+  logical :: USE_CMT_AND_FORCE_SOURCE = .false.
   logical :: USE_BINARY_SOURCE_FILE = .false.
 
   ! external source time function
@@ -234,7 +234,7 @@ module shared_compute_parameters
   integer :: NSOURCES
 
   !number of sources in CMTSOLUTION/FORCESOLUTION
-  integer :: NSOURCES_CMT, NSOURCES_FORCE 
+  integer :: NSOURCES_CMT, NSOURCES_FORCE
 
   ! anchor points
   integer :: NGNOD2D
diff --git a/src/specfem3D/comp_source_time_function.f90 b/src/specfem3D/comp_source_time_function.f90
@@ -400,7 +400,7 @@ end function comp_source_time_function_mono
   double precision function comp_source_time_function_ext(it_index,isource)
 
   use specfem_par, only: myrank, NSTEP, user_source_time_function
-  use specfem_par, only : USE_BINARY_SOURCE_FILE,isource_glob2loc
+  use specfem_par, only: USE_BINARY_SOURCE_FILE,isource_glob2loc
 
   implicit none
 
@@ -418,7 +418,7 @@ double precision function comp_source_time_function_ext(it_index,isource)
   endif
 
   ! gets stored STF
-  if(.not. USE_BINARY_SOURCE_FILE) then 
+  if (.not. USE_BINARY_SOURCE_FILE) then
     comp_source_time_function_ext = user_source_time_function(it_index,isource)
   else
     comp_source_time_function_ext = user_source_time_function(it_index,isource_glob2loc(isource))
diff --git a/src/specfem3D/compute_add_sources_acoustic.f90 b/src/specfem3D/compute_add_sources_acoustic.f90
@@ -437,7 +437,7 @@ subroutine compute_add_sources_acoustic_GPU()
   use specfem_par, only: FAULT_SIMULATION
 
   !
-  use specfem_par,only : islice_selected_source
+  use specfem_par, only: islice_selected_source
 
   implicit none
 
@@ -483,7 +483,7 @@ subroutine compute_add_sources_acoustic_GPU()
 
       ! determines source time function value
       stf = 0.
-      if(myrank == islice_selected_source(isource)) &
+      if (myrank == islice_selected_source(isource)) &
         stf = get_stf_acoustic(time_source_dble,isource,it)
 
       ! stores precomputed source time function factor
@@ -602,7 +602,7 @@ subroutine compute_add_sources_acoustic_backward_GPU()
   ! faults
   use specfem_par, only: FAULT_SIMULATION
 
-  use specfem_par,only: islice_selected_source
+  use specfem_par, only: islice_selected_source
 
   implicit none
 
@@ -703,8 +703,8 @@ subroutine compute_add_sources_acoustic_backward_GPU()
     time_source_dble = time_t - tshift_src(isource)
 
     ! determines source time function value
-    stf = 0.0d0 
-    if(myrank == islice_selected_source(isource)) &
+    stf = 0.0d0
+    if (myrank == islice_selected_source(isource)) &
       stf = get_stf_acoustic(time_source_dble,isource,NSTEP-it_tmp+1)
 
     ! stores precomputed source time function factor
@@ -762,7 +762,7 @@ double precision function get_stf_acoustic(time_source_dble,isource,it_tmp_ext)
 
   ! determines source time function value
   !NQDU if (USE_FORCE_POINT_SOURCE) then
-  if(is_POINTFORCE(isource)) then 
+  if (is_POINTFORCE(isource)) then
     ! single point force
     select case(force_stf(isource))
     case (0)
diff --git a/src/specfem3D/compute_add_sources_poroelastic.f90 b/src/specfem3D/compute_add_sources_poroelastic.f90
@@ -135,7 +135,7 @@ subroutine compute_add_sources_poroelastic()
                 ! we distinguish between a single force which can be applied both in fluid and solid
                 ! and a moment-tensor source which only makes sense for a solid
                 !if (USE_FORCE_POINT_SOURCE) then
-                if(is_POINTFORCE(isource)) then 
+                if (is_POINTFORCE(isource)) then
                   ! single point force
                   ! the source is applied to both solid and fluid phase: bulk source.
                   fac_s = 1._CUSTOM_REAL - phil/tortl
@@ -363,7 +363,7 @@ subroutine compute_add_sources_poroelastic()
                 ! we distinguish between a single force which can be applied both in fluid and solid
                 ! and a moment-tensor source which only makes sense for a solid
                 !if (USE_FORCE_POINT_SOURCE) then
-                if(is_POINTFORCE(isource)) then 
+                if (is_POINTFORCE(isource)) then
                   ! single point force
                   ! the source is applied to both solid and fluid phase: bulk source.
                   fac_s = 1._CUSTOM_REAL - phil/tortl
@@ -434,7 +434,7 @@ double precision function get_stf_poroelastic(time_source_dble,isource,it_tmp_ex
 
   ! determines source time function value
   !if (USE_FORCE_POINT_SOURCE) then
-  if(is_POINTFORCE(isource)) then 
+  if (is_POINTFORCE(isource)) then
     ! single point force
     select case(force_stf(isource))
     case (0)
diff --git a/src/specfem3D/compute_add_sources_viscoelastic.F90 b/src/specfem3D/compute_add_sources_viscoelastic.F90
@@ -495,8 +495,8 @@ subroutine compute_add_sources_viscoelastic_GPU()
     islice_selected_rec, &
     nadj_rec_local,source_adjoint,nadj_rec_local,number_adjsources_global, &
     Mesh_pointer
-  
-  use specfem_par,only : islice_selected_source
+
+  use specfem_par, only: islice_selected_source
 
   use specfem_par_noise, only: irec_main_noise,noise_surface_movie
 
@@ -572,7 +572,7 @@ subroutine compute_add_sources_viscoelastic_GPU()
 
         ! determines source time function value
         stf = 0.0d0
-        if(myrank == islice_selected_source(isource)) &
+        if (myrank == islice_selected_source(isource)) &
           stf = get_stf_viscoelastic(time_source_dble,isource,it)
 
         ! LTS
@@ -733,7 +733,7 @@ subroutine compute_add_sources_viscoelastic_backward_GPU()
   ! faults
   use specfem_par, only: FAULT_SIMULATION
 
-  use specfem_par,only : islice_selected_source
+  use specfem_par, only: islice_selected_source
 
   implicit none
 
@@ -807,7 +807,7 @@ subroutine compute_add_sources_viscoelastic_backward_GPU()
 
       ! determines source time function value
       stf = 0.0d0
-      if(myrank == islice_selected_source(isource)) &
+      if (myrank == islice_selected_source(isource)) &
         stf = get_stf_viscoelastic(time_source_dble,isource,NSTEP-it_tmp+1)
 
       ! stores precomputed source time function factor
@@ -880,7 +880,7 @@ double precision function get_stf_viscoelastic(time_source_dble,isource,it_tmp_e
 
   ! determines source time function value
   !if (USE_FORCE_POINT_SOURCE) then
-  if(is_POINTFORCE(isource)) then 
+  if (is_POINTFORCE(isource)) then
     ! single point force
     select case(force_stf(isource))
     case (0)
diff --git a/src/specfem3D/finalize_simulation.f90 b/src/specfem3D/finalize_simulation.f90
@@ -184,8 +184,8 @@ subroutine finalize_simulation_cleanup()
     deallocate(utm_x_source,utm_y_source)
     deallocate(nu_source)
     deallocate(user_source_time_function)
-    if(allocated(is_POINTFORCE)) deallocate(is_POINTFORCE)
-    if(allocated(isource_glob2loc)) deallocate(isource_glob2loc)
+    if (allocated(is_POINTFORCE)) deallocate(is_POINTFORCE)
+    if (allocated(isource_glob2loc)) deallocate(isource_glob2loc)
   endif
 
   ! receivers
diff --git a/src/specfem3D/get_cmt.f90 b/src/specfem3D/get_cmt.f90
diff --git a/src/specfem3D/locate_source.F90 b/src/specfem3D/locate_source.F90
diff --git a/src/specfem3D/setup_sources_receivers.f90 b/src/specfem3D/setup_sources_receivers.f90
diff --git a/src/tomography/postprocess_sensitivity_kernels/smooth_sem_pde.F90 b/src/tomography/postprocess_sensitivity_kernels/smooth_sem_pde.F90
