explicitly specifies size(**,kind=4) in function arguments (to avoid compiler problems where size() returns either integer(kind=4) or integer(kind=8) depending on flags like -mcmedium)

Author: danielpeter

src/generate_databases/model_external_values.F90

@@ -94,7 +94,7 @@ subroutine model_external_broadcast()
   if (myrank == 0) call read_external_model()
 
   ! broadcast the information read on the main to the nodes
-  call bcast_all_dp(MEXT_V%dvs, size(MEXT_V%dvs))
+  call bcast_all_dp(MEXT_V%dvs, size(MEXT_V%dvs,kind=4))
 
   end subroutine model_external_broadcast
 

src/generate_databases/model_salton_trough.f90

@@ -71,7 +71,7 @@ subroutine model_salton_trough_broadcast()
   if (myrank == 0) call read_salton_sea_model()
 
   ! broadcast the information read on the main to the nodes
-  call bcast_all_r(vp_array, size(vp_array))
+  call bcast_all_r(vp_array, size(vp_array,kind=4))
 
   end subroutine model_salton_trough_broadcast
 

src/generate_databases/model_tomography.f90

@@ -121,7 +121,7 @@ subroutine model_tomography_broadcast(myrank)
   ! allocate( vp_tomography(1:nrecord) ,stat=ier)
   ! if (ier /= 0) stop 'error allocating array vp_tomography'
   !endif
-  !call bcast_all_cr(vp_tomography,size(vp_tomography))
+  !call bcast_all_cr(vp_tomography,size(vp_tomography,kind=4))
 
   ! synchronizes processes
   call synchronize_all()

src/generate_databases/setup_mesh_adjacency.f90

@@ -377,7 +377,7 @@ subroutine setup_mesh_adjacency()
   endif
 
   ! check if element has neighbors
-  ! note: in case of a fault in this slice (splitting nodes) and/or scotch paritioning
+  ! note: in case of a fault in this slice (splitting nodes) and/or scotch partitioning
   !       it can happen that an element has no neighbors
   if (NPROC == 1 .and. (.not. ANY_FAULT_IN_THIS_PROC)) then
     ! checks if neighbors were found

src/meshfem3D/determine_cavity.f90

@@ -366,13 +366,13 @@ subroutine cmm_determine_cavity(nglob)
         allocate(tmp_all(1,1),stat=ier)
         if (ier /= 0) call exit_MPI_without_rank('error allocating array 1330')
       endif
-      call sum_all_1Darray_dp(cavity_boundary,tmp_all,size(cavity_boundary))
+      call sum_all_1Darray_dp(cavity_boundary,tmp_all,size(cavity_boundary,kind=4))
       if (myrank == 0) then
         cavity_boundary(:,:) = tmp_all(:,:)
       endif
       deallocate(tmp_all)
       ! broadcasts to all others
-      call bcast_all_dp(cavity_boundary,size(cavity_boundary))
+      call bcast_all_dp(cavity_boundary,size(cavity_boundary,kind=4))
 
       !print *,'cavity boundary after:',myrank,'array:',cavity_boundary(:,:)
 

src/specfem3D/iterate_time.F90

@@ -418,9 +418,9 @@ subroutine it_transfer_from_GPU()
 
       if (ATTENUATION) then
         call transfer_rmemory_from_device(Mesh_pointer,R_xx,R_yy,R_xy,R_xz,R_yz, &
-                                          R_trace,size(R_xx))
+                                          R_trace,size(R_xx,kind=4))
         call transfer_strain_from_device(Mesh_pointer,epsilondev_xx,epsilondev_yy,epsilondev_xy,epsilondev_xz,epsilondev_yz, &
-                                         epsilondev_trace,size(epsilondev_xx))
+                                         epsilondev_trace,size(epsilondev_xx,kind=4))
 
       endif
     endif

src/specfem3D/prepare_attenuation.f90

@@ -114,12 +114,12 @@ subroutine prepare_attenuation()
 
   ! broadcasts
   call bcast_all_i_for_database(ispec, 1)
-  if (size(factor_common) > 0) &
-    call bcast_all_cr_for_database(factor_common(1,1,1,1,1), size(factor_common))
-  if (size(scale_factor) > 0) &
-    call bcast_all_cr_for_database(scale_factor(1,1,1,1), size(scale_factor))
-  call bcast_all_cr_for_database(factor_common_kappa(1,1,1,1,1), size(factor_common_kappa))
-  call bcast_all_cr_for_database(scale_factor_kappa(1,1,1,1), size(scale_factor_kappa))
+  if (size(factor_common,kind=4) > 0) &
+    call bcast_all_cr_for_database(factor_common(1,1,1,1,1), size(factor_common,kind=4))
+  if (size(scale_factor,kind=4) > 0) &
+    call bcast_all_cr_for_database(scale_factor(1,1,1,1), size(scale_factor,kind=4))
+  call bcast_all_cr_for_database(factor_common_kappa(1,1,1,1,1), size(factor_common_kappa,kind=4))
+  call bcast_all_cr_for_database(scale_factor_kappa(1,1,1,1), size(scale_factor_kappa,kind=4))
 
   ! gets stress relaxation times tau_sigma, i.e.
   ! precalculates tau_sigma depending on period band (constant for all Q_mu), and

src/specfem3D/prepare_gpu.f90

@@ -145,7 +145,7 @@ subroutine prepare_GPU()
                                        epsilondev_xx,epsilondev_yy,epsilondev_xy, &
                                        epsilondev_xz,epsilondev_yz, &
                                        ATTENUATION, &
-                                       size(R_xx), &
+                                       size(R_xx,kind=4), &
                                        R_xx,R_yy,R_xy,R_xz,R_yz, &
                                        factor_common, &
                                        R_trace,epsilondev_trace, &
@@ -171,7 +171,7 @@ subroutine prepare_GPU()
                                           b_epsilondev_xx,b_epsilondev_yy,b_epsilondev_xy, &
                                           b_epsilondev_xz,b_epsilondev_yz, &
                                           b_epsilon_trace_over_3, &
-                                          size(R_xx), &
+                                          size(R_xx,kind=4), &
                                           b_R_xx,b_R_yy,b_R_xy,b_R_xz,b_R_yz, &
                                           b_R_trace,b_epsilondev_trace, &
                                           b_alphaval,b_betaval,b_gammaval, &

src/specfem3D/read_forward_arrays.f90

@@ -110,11 +110,11 @@ subroutine read_forward_arrays()
       if (ATTENUATION) then
         call transfer_b_rmemory_to_device(Mesh_pointer, &
                                           b_R_xx,b_R_yy,b_R_xy,b_R_xz,b_R_yz, &
-                                          b_R_trace,size(b_R_xx))
+                                          b_R_trace,size(b_R_xx,kind=4))
         call transfer_b_strain_to_device(Mesh_pointer, &
                                          b_epsilondev_xx,b_epsilondev_yy,b_epsilondev_xy, &
                                          b_epsilondev_xz,b_epsilondev_yz, &
-                                         b_epsilondev_trace,size(b_epsilondev_xx))
+                                         b_epsilondev_trace,size(b_epsilondev_xx,kind=4))
       endif
     endif
 
@@ -217,7 +217,7 @@ subroutine read_forward_arrays_undoatt()
       if (ATTENUATION) then
         call transfer_b_rmemory_to_device(Mesh_pointer, &
                                           b_R_xx,b_R_yy,b_R_xy,b_R_xz,b_R_yz, &
-                                          b_R_trace,size(b_R_xx))
+                                          b_R_trace,size(b_R_xx,kind=4))
       endif
     endif