[cig-commits] [commit] devel, master: Fixed viscous attenuation for poroelastic media to suppress the error on inversion of the permeability matrix in specfem2D.f90. Modified some parameters in Par_file_acoustic_poroelastic so users can turn ATTENUATION_PORO_FLUID_PART on or off to use or not viscous attenuation. (d5f0734)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branches: devel,master
Link       : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c

>---------------------------------------------------------------

commit d5f0734f06d39c819da81e8d8487754f536ec70a
Author: Christina Morency <morency1 at llnl.gov>
Date:   Wed Feb 27 19:23:06 2013 +0000

    Fixed viscous attenuation for poroelastic media to suppress the error on inversion of the permeability matrix in
     specfem2D.f90.
    Modified some parameters in Par_file_acoustic_poroelastic so users can turn ATTENUATION_PORO_FLUID_PART on or off to use or not viscous attenuation.


>---------------------------------------------------------------

d5f0734f06d39c819da81e8d8487754f536ec70a
 acoustic_poroelastic/Par_file_acoustic_poroelastic | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index b345abe..94259d9 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -19,14 +19,14 @@ READ_EXTERNAL_SEP_FILE          = .false.        # Read external SEP file from D
 ATTENUATION_VISCOELASTIC_SOLID  = .false.        # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
 ATTENUATION_PORO_FLUID_PART     = .false.        # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
 Q0                              =  1             # quality factor for viscous attenuation
-freq0                           =  1000           # frequency for viscous attenuation
+freq0                           =  100           # frequency for viscous attenuation
 p_sv                            = .true.         # set the type of calculation (P-SV or SH/membrane waves)
 
 # time step parameters
 nt                              = 2500           # total number of time steps
 deltat                          = 4.d-4         # duration of a time step (see section "How to choose the time step" of the manual for how to do this)
 USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
-time_stepping_scheme            = 2   # 1 = Newmark (2nd order),     2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta),     3 = classical 4th-order 4-stage Runge-Kutta
+time_stepping_scheme            = 1   # 1 = Newmark (2nd order),     2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta),     3 = classical 4th-order 4-stage Runge-Kutta
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -109,7 +109,7 @@ nbmodels                        = 2              # nb of different models
 # III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
 # For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
 # and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
-1 3 2650.d0 880.d0 0.3d0 1.2 1d-9 0.d0 1d-9 1.22d10 1.985d9 9.6d9 10d-4 0.0001 9999999999
+1 3 2650.d0 880.d0 0.3d0 1.2 1d-9 0.d0 1d-9 1.22d10 1.985d9 9.6d9 10d-4 0.001 9999999999
 2 1 1020.d0 1500.d0 0.d0 0 0 10.d0 10.d0 0 0 0 0 0 0
 
 # external mesh or not