code cleaning (whitespace)

danielpeter · Jun 5, 2024 · 0e64075 · 0e64075
1 parent ca6d753
commit 0e64075
Show file tree

Hide file tree

Showing 18 changed files with 97 additions and 97 deletions.
diff --git a/src/shared/read_material_table.f90 b/src/shared/read_material_table.f90
@@ -45,7 +45,7 @@ subroutine read_material_table()
                                rho_s_read,rho_f_read, &
                                phi_read,tortuosity_read, &
                                permxx_read,permxz_read,permzz_read,kappa_s_read,kappa_f_read,kappa_fr_read, &
-                               eta_f_read,mu_fr_read,mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,&
+                               eta_f_read,mu_fr_read,mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, &
                                Qe11_read,Qe33_read,Qs11_read,Qs33_read, &
                                compaction_grad
 
@@ -167,7 +167,7 @@ subroutine read_material_table()
         call stop_the_code('Invalid poroelastic material: non-positive value of Qmu')
 
     else if (icodemat(i) == ELECTROMAGNETIC_MATERIAL) then
-      ! electromagnetic material 
+      ! electromagnetic material
       mu0_read(i) = val0read
       e0_read(i) = val1read
       e11_read(i) = val2read
@@ -248,7 +248,7 @@ subroutine initialize_material_properties()
                                eta_f_read,mu_fr_read
 
   ! electromagnetic
-  use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,&
+  use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, &
                                Qe11_read,Qe33_read,Qs11_read,Qs33_read
 
   implicit none
@@ -289,15 +289,15 @@ subroutine initialize_material_properties()
            eta_f_read(nbmodels), &
            mu_fr_read(nbmodels))
 
-  allocate( mu0_read(nbmodels),&
-            e0_read(nbmodels),&
-            e11_read(nbmodels),&
-            e33_read(nbmodels),&
-            sig11_read(nbmodels),&
-            sig33_read(nbmodels),&
-            Qe11_read(nbmodels),&
-            Qe33_read(nbmodels),&
-            Qs11_read(nbmodels),&
+  allocate( mu0_read(nbmodels), &
+            e0_read(nbmodels), &
+            e11_read(nbmodels), &
+            e33_read(nbmodels), &
+            sig11_read(nbmodels), &
+            sig33_read(nbmodels), &
+            Qe11_read(nbmodels), &
+            Qe33_read(nbmodels), &
+            Qs11_read(nbmodels), &
             Qs33_read(nbmodels))
 
   allocate(compaction_grad(nbmodels))
@@ -379,7 +379,7 @@ subroutine print_materials_info()
                                eta_f_read,mu_fr_read
 
   ! electromagnetic
-  use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,&
+  use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, &
                                Qe11_read,Qe33_read,Qs11_read,Qs33_read
 
   implicit none

diff --git a/src/shared/read_parameter_file.F90 b/src/shared/read_parameter_file.F90
@@ -564,14 +564,14 @@ subroutine read_parameter_file_only()
     write(*,'(a)') 'ATTENUATION_PERMITTIVITY        = .false.'
     write(*,*)
   endif
-    
+
   call read_value_logical_p(ATTENUATION_CONDUCTIVITY, 'ATTENUATION_CONDUCTIVITY')
   if (err_occurred() /= 0) then
     some_parameters_missing_from_Par_file = .true.
     write(*,'(a)') 'ATTENUATION_CONDUCTIVITY       = .false.'
     write(*,*)
   endif
-    
+
   call read_value_double_precision_p(f0_electromagnetic, 'f0_electromagnetic')
   if (err_occurred() /= 0) then
     some_parameters_missing_from_Par_file = .true.

diff --git a/src/shared/shared_par.F90 b/src/shared/shared_par.F90
@@ -416,7 +416,7 @@ module shared_parameters
        permzz_read,kappa_s_read,kappa_f_read,kappa_fr_read,eta_f_read,mu_fr_read
 
   ! electromagnetic (EM)
-  double precision, dimension(:),allocatable :: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,&
+  double precision, dimension(:),allocatable :: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, &
                               Qe11_read,Qe33_read,Qs11_read,Qs33_read
 
   ! compaction gradient

diff --git a/src/specfem2D/check_grid.F90 b/src/specfem2D/check_grid.F90
@@ -218,7 +218,7 @@ subroutine check_grid()
            cpIIloc = 0.d0
          endif
          cEloc = 0.d0
- 
+
         else
          cpIloc = 0.d0
          cpIIloc = 0.d0
@@ -241,7 +241,7 @@ subroutine check_grid()
         ! ignore fluid regions with Vs = 0
         if (csloc > TINYVAL) vsmin = min(vsmin,csloc)
         vsmax = max(vsmax,csloc)
-        
+
         if (cEloc > TINYVAL) vEmin = min(vEmin,cEloc)
         vEmax = max(vEmax,cEloc)
 

diff --git a/src/specfem2D/compute_add_sources_electromagnetic.f90 b/src/specfem2D/compute_add_sources_electromagnetic.f90
@@ -38,7 +38,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage)
   use constants, only: CUSTOM_REAL,NGLLX,NGLLZ,NDIM,HALF,myrank
 
   use specfem_par, only: P_SV,ispec_is_electromagnetic,nglob_electromagnetic, &
-                         NSOURCES,source_time_function,&
+                         NSOURCES,source_time_function, &
                          islice_selected_source,ispec_selected_source,sourcearrays, &
                          ibool
   implicit none
@@ -58,7 +58,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage)
 
       ! element containing source
       ispec = ispec_selected_source(i_source)
-      
+
       ! source element is electromagnetic
       if (ispec_is_electromagnetic(ispec)) then
         ! source time function
@@ -83,7 +83,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage)
            ! 2D: y-component uses array(1,..)
            accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - sourcearrays(1,i,j,i_source) * stf_used
             enddo
-          enddo         
+          enddo
        endif
 
       endif ! source element is electromagnetic

diff --git a/src/specfem2D/compute_attenuation_electromagnetic.f90 b/src/specfem2D/compute_attenuation_electromagnetic.f90
@@ -66,7 +66,7 @@ subroutine compute_attenuation_electromagnetic()
           permlxx = spermittivitystore(1,i,j,ispec)
           permlzz = spermittivitystore(2,i,j,ispec)
 
-          permx_loc(i,j) = veloc_electromagnetic(1,iglob) * permlxx 
+          permx_loc(i,j) = veloc_electromagnetic(1,iglob) * permlxx
           permz_loc(i,j) = veloc_electromagnetic(2,iglob) * permlzz
 
           ! time stepping

diff --git a/src/specfem2D/compute_forces_electromagnetic.F90 b/src/specfem2D/compute_forces_electromagnetic.F90
@@ -31,14 +31,14 @@
 !========================================================================
 
 
-  subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic,&
+  subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic, &
                                             rx_permattenuation,rz_permattenuation,iphase)
 
   ! compute forces for the electromagnetic elements
   use constants, only: CUSTOM_REAL,NGLLX,NGLLZ,NGLJ,NDIM, &
     ONE,HALF,PI,TINYVAL,FOUR_THIRDS
 
-  use specfem_par, only: nglob,&
+  use specfem_par, only: nglob, &
                          ibool,kmato,ispec_is_electromagnetic, &
                          xix,xiz,gammax,gammaz,jacobian, &
                          hprime_xx,hprimewgll_xx,hprime_zz,hprimewgll_zz,wxgll,wzgll, &
@@ -177,7 +177,7 @@ subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromag
         !--- if external medium, get electromagnetic parameters of current grid point
         inv_mag_permeability = inv_magpermeabilitystore(i,j,ispec)
 
-        if(P_SV) then
+        if (P_SV) then
         sigma_xx = inv_mag_permeability*(dux_dxl(i,j)+duz_dzl(i,j)) !
         sigma_zz = inv_mag_permeability*(dux_dxl(i,j)+duz_dzl(i,j)) !
         sigma_xz = inv_mag_permeability*(duz_dxl(i,j) - dux_dzl(i,j))  !  d_t \tilde{H}_xz = d_t H_y = mu^-1(dxEz-dzEx)
@@ -255,7 +255,7 @@ subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromag
             accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - (wzgll(j) * tempx1l + wxgll(i) * tempx2l) &
                               - wxgll(i) * wzgll(j) * jacobian(i,j,ispec) * rx_permattenuation(i,j,ispec)
             accel_electromagnetic(2,iglob) = 0.d0
-           endif 
+           endif
          else
            if (P_SV) then
             accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - (wzgll(j) * tempx1l + wxgll(i) * tempx2l)

diff --git a/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90 b/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90
@@ -55,7 +55,7 @@ subroutine compute_forces_electromagnetic_main()
 
     ! main solver for the electromagnetic elements
     ! electromagnetic term
-    call compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic,&
+    call compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic, &
                                         rx_permattenuation,rz_permattenuation,iphase)
     !                                 PML_BOUNDARY_CONDITIONS,e1,e11,e13,iphase)
 
@@ -73,7 +73,7 @@ subroutine compute_forces_electromagnetic_main()
 
 !      ! PML boundary
 !      if (PML_BOUNDARY_CONDITIONS) then
-!        call pml_boundary_electromagnetic(accel_electromagnetic,veloc_electromagnetic,displ_electromagnetic,&
+!        call pml_boundary_electromagnetic(accel_electromagnetic,veloc_electromagnetic,displ_electromagnetic, &
 !          displ_electromagnetic_old)
 !      endif
 

diff --git a/src/specfem2D/compute_stacey_electromagnetic.f90 b/src/specfem2D/compute_stacey_electromagnetic.f90
@@ -50,7 +50,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
                          b_absorb_electromagnetic_left,b_absorb_electromagnetic_right, &
                          b_absorb_electromagnetic_bottom,b_absorb_electromagnetic_top, &
                          ib_left,ib_right,ib_bottom,ib_top, &
-                         anyabs,STACEY_ABSORBING_CONDITIONS,Qe11_electromagnetic,Qe33_electromagnetic,&
+                         anyabs,STACEY_ABSORBING_CONDITIONS,Qe11_electromagnetic,Qe33_electromagnetic, &
                          Qs11_electromagnetic,Qs33_electromagnetic,P_SV
 
   implicit none
@@ -101,14 +101,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
         permlxx = spermittivitystore(1,i,j,ispec)  !e11
         permlzz = spermittivitystore(2,i,j,ispec)  !e33
         condlxx = sconductivitystore(1,i,j,ispec)  !sig11
-        condlzz = sconductivitystore(2,i,j,ispec)  !sig33 
+        condlzz = sconductivitystore(2,i,j,ispec)  !sig33
         two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1
-      
+
         call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, &
-         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),&
+         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), &
          Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), &
          permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability)
-  
+
         cpxl = sqrt(cpxsquare)
         cpzl = sqrt(cpzsquare)
 
@@ -117,7 +117,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
 
         if (P_SV) then
           ! P_SV case
-          vx = veloc_electromagnetic(1,iglob) 
+          vx = veloc_electromagnetic(1,iglob)
           vz = veloc_electromagnetic(2,iglob)
           vn = nx*vx + nz*vz
 
@@ -128,7 +128,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
           ! with velocity's tangential component: veloc_tangential = v - vn * normal
           ! total traction
           !    T = T_normal + T_tangential
-          tx = rho_cpx*vn*nx 
+          tx = rho_cpx*vn*nx
           tz = rho_cpz*vn*nz
 
         else
@@ -178,14 +178,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
         permlxx = spermittivitystore(1,i,j,ispec)  !e11
         permlzz = spermittivitystore(2,i,j,ispec)  !e33
         condlxx = sconductivitystore(1,i,j,ispec)  !sig11
-        condlzz = sconductivitystore(2,i,j,ispec)  !sig33 
+        condlzz = sconductivitystore(2,i,j,ispec)  !sig33
         two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1
-      
+
         call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, &
-         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),&
+         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), &
          Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), &
          permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability)
-  
+
         cpxl = sqrt(cpxsquare)
         cpzl = sqrt(cpzsquare)
 
@@ -194,7 +194,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
 
         if (P_SV) then
           ! P_SV case
-          vx = veloc_electromagnetic(1,iglob) 
+          vx = veloc_electromagnetic(1,iglob)
           vz = veloc_electromagnetic(2,iglob)
           vn = nx*vx + nz*vz
           tx = rho_cpx*vn*nx
@@ -215,7 +215,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
           ! SH case
           accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - ty*weight
         endif
- 
+
         if (SAVE_FORWARD .and. SIMULATION_TYPE == 1) then
          if (P_SV) then
           ! P_SV case
@@ -245,14 +245,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
         permlxx = spermittivitystore(1,i,j,ispec)  !e11
         permlzz = spermittivitystore(2,i,j,ispec)  !e33
         condlxx = sconductivitystore(1,i,j,ispec)  !sig11
-        condlzz = sconductivitystore(2,i,j,ispec)  !sig33 
+        condlzz = sconductivitystore(2,i,j,ispec)  !sig33
         two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1
-      
+
         call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, &
-         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),&
+         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), &
          Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), &
          permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability)
-  
+
         cpxl = sqrt(cpxsquare)
         cpzl = sqrt(cpzsquare)
 
@@ -261,7 +261,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
 
         if (P_SV) then
           ! P_SV case
-          vx = veloc_electromagnetic(1,iglob) 
+          vx = veloc_electromagnetic(1,iglob)
           vz = veloc_electromagnetic(2,iglob)
           vn = nx*vx+nz*vz
           tx = rho_cpx*vn*nx
@@ -326,14 +326,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
         permlxx = spermittivitystore(1,i,j,ispec)  !e11
         permlzz = spermittivitystore(2,i,j,ispec)  !e33
         condlxx = sconductivitystore(1,i,j,ispec)  !sig11
-        condlzz = sconductivitystore(2,i,j,ispec)  !sig33 
+        condlzz = sconductivitystore(2,i,j,ispec)  !sig33
         two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1
-      
+
         call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, &
-         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),&
+         ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), &
          Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), &
          permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability)
-  
+
         cpxl = sqrt(cpxsquare)
         cpzl = sqrt(cpzsquare)
 
@@ -342,7 +342,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag
 
         if (P_SV) then
           ! P_SV case
-          vx = veloc_electromagnetic(1,iglob) 
+          vx = veloc_electromagnetic(1,iglob)
           vz = veloc_electromagnetic(2,iglob)
           vn = nx*vx+nz*vz
           tx = rho_cpx*vn*nx