Try

Author: bottero

EXAMPLES/axisymmetric_case_AXISYM_option/Par_file_axisym

@@ -10,7 +10,7 @@ NOISE_TOMOGRAPHY                = 0  # 0 = earthquake simulation, 1/2/3 = noise
 SAVE_FORWARD                    = .false. # save the last frame, needed for adjoint simulation
 
 # parameters concerning partitioning
-nproc                           = 20              # number of processes
+nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
 
 ngnod                           = 4              # number of control nodes per element (4 or 9)

src/specfem2D/compute_forces_viscoelastic.F90

@@ -696,7 +696,7 @@ subroutine compute_forces_viscoelastic(accel_elastic,veloc_elastic,displ_elastic
             ! When implementing viscoelasticity according to the Carcione 1993 paper, attenuation is
             ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
             ! 2004 paper and his 2007 book. See also file doc/problem_attenuation_reference_Specfem2D_fixed_by_Xie_Zhinan.pdf
-            
+
             if(AXISYM) then
               if (is_on_the_axis(ispec)) then
                 if (abs(coord(1,ibool(i,j,ispec))) < TINYVAL) then ! First GLJ point
@@ -734,7 +734,7 @@ subroutine compute_forces_viscoelastic(accel_elastic,veloc_elastic,displ_elastic
               sigma_xz = mul_unrelaxed_elastic*(duz_dxl + dux_dzl)
               sigma_zz = lambdaplus2mu_unrelaxed_elastic*duz_dzl + lambdal_unrelaxed_elastic*dux_dxl
             endif
-            
+
             ! add the memory variables using the relaxed parameters (Carcione 2007 page 125)
             ! beware: there is a bug in Carcione's equation (2c) of his 1993 paper for sigma_zz, we fixed it in the code below.
             ! When implementing viscoelasticity according to the Carcione 1993 paper, attenuation is

src/specfem2D/compute_gradient_attenuation.f90

@@ -52,8 +52,6 @@ subroutine compute_gradient_attenuation(displ_elastic,dux_dxl,duz_dxl,dux_dzl,du
 
   include "constants.h"
 
-  !logical :: AXISYM
-  
   integer :: nspec,nglob
 
   integer, dimension(NGLLX,NGLLZ,nspec) :: ibool