code cleaning (fortran linter fortitude)

Author: danielpeter

src/meshfem3D/model_1066a.f90

@@ -35,6 +35,8 @@
 
   module model_1066a_par
 
+  implicit none
+
   ! number of layers in DATA/1066a/1066a.dat
   integer, parameter :: NR_1066A = 160
 

src/meshfem3D/model_1dberkeley.f90

@@ -579,7 +579,7 @@ subroutine determine_1dberkeley_moho_depth(moho_depth)
 
   ! local parameters
   double precision :: moho_radius
-  double precision :: earthradius = 6371.d0
+  double precision, parameter :: earthradius = 6371.d0
 
   ! get moho radius (in km)
   call determine_1dberkeley_moho_radius(moho_radius)

src/meshfem3D/model_1dref.f90

@@ -53,6 +53,8 @@
 
   module model_1dref_par
 
+  implicit none
+
   ! number of layers in DATA/s362ani/REF
   integer, parameter :: NR_REF = 750
 

src/meshfem3D/model_ak135.f90

@@ -54,6 +54,8 @@
 
   module model_ak135_par
 
+  implicit none
+
   ! number of layers in DATA/ak135/ak135.dat
   integer, parameter :: NR_AK135F_NO_MUD = 136
 

src/meshfem3D/model_aniso_mantle.f90

@@ -42,6 +42,8 @@
 
   module model_aniso_mantle_par
 
+  implicit none
+
   ! model_aniso_mantle_variables
   double precision,dimension(:,:,:,:),allocatable :: AMM_V_beta
   double precision,dimension(:),allocatable :: AMM_V_pro

src/meshfem3D/model_atten3D_QRFSI12.f90

@@ -49,6 +49,8 @@
 
   module model_atten3D_QRFSI12_par
 
+  implicit none
+
   ! QRFSI12 constants
   integer,parameter :: NKQ = 8,MAXL_Q = 12
   integer,parameter :: NSQ=(MAXL_Q+1)**2,NDEPTHS_REFQ = 913
@@ -110,7 +112,7 @@ subroutine model_atten3D_QRFSI12_broadcast()
   if (ier /= 0) stop 'Error allocating helper array'
   xlmvec(:) = 0.0
 
-  end subroutine
+  end subroutine model_atten3D_QRFSI12_broadcast
 
 !
 !-------------------------------------------------------------------------------------------------

src/meshfem3D/model_attenuation_gll.f90

@@ -37,6 +37,8 @@ module model_gll_qmu_par
 
   use constants, only: CUSTOM_REAL
 
+  implicit none
+
   ! GLL model_variables
   type model_gll_qmu_variables
     !TODO: check if `sequence` is needed

src/meshfem3D/model_case65TAY.f90

@@ -46,6 +46,8 @@
 
   module model_case65tay_par
 
+  implicit none
+
   ! number of layers in DATA/case65TAY/case65TAY.dat
   integer, parameter :: NR_case65TAY = 681
 

src/meshfem3D/model_crust_1_0.f90

@@ -55,6 +55,8 @@
 
   module model_crust_1_0_par
 
+  implicit none
+
   ! crustal_model_constants
   ! crustal model parameters for crust1.0
   integer, parameter :: CRUST_NP  = 9

src/meshfem3D/model_crust_2_0.f90

@@ -48,6 +48,8 @@
 
   module model_crust_2_0_par
 
+  implicit none
+
   ! crustal_model_constants
   ! crustal model parameters for crust2.0
   integer, parameter :: CRUST_NP = 8