[cig-commits] r21304 - in seismo/3D/SPECFEM3D_GLOBE/trunk/src: shared specfem3D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Tue Jan 29 09:56:35 PST 2013
Author: dkomati1
Date: 2013-01-29 09:56:35 -0800 (Tue, 29 Jan 2013)
New Revision: 21304
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_element.F90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle_Dev.F90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core_Dev.F90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/initialize_simulation.f90
seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/specfem3D.F90
Log:
changed USE_ATTENUATION_MIMIC to USE_PHYSICAL_DISPERSION_ONLY
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/save_header_file.f90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -539,15 +539,15 @@
! mimicking effect of attenuation on apparent velocities, not amplitudes. that is,
! phase shifts should be correctly accounted for, but amplitudes will differ in adjoint simulations
if( ATTENUATION ) then
- write(IOUT,*) 'logical, parameter :: USE_ATTENUATION_MIMIC = .true.'
+ write(IOUT,*) 'logical, parameter :: USE_PHYSICAL_DISPERSION_ONLY = .true.'
else
- write(IOUT,*) 'logical, parameter :: USE_ATTENUATION_MIMIC = .false.'
+ write(IOUT,*) 'logical, parameter :: USE_PHYSICAL_DISPERSION_ONLY = .false.'
endif
else
! calculates full attenuation (phase & amplitude effects) if used
- write(IOUT,*) 'logical, parameter :: USE_ATTENUATION_MIMIC = .false.'
+ write(IOUT,*) 'logical, parameter :: USE_PHYSICAL_DISPERSION_ONLY = .false.'
endif
! attenuation and/or adjoint simulations
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_element.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_element.F90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_element.F90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -219,7 +219,7 @@
sigma_yz = mul*duzdyl_plus_duydzl
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! note: fortran passes pointers to array location, thus R_memory(1,1,...) should be fine
call compute_element_att_stress( R_memory(1,1,i,j,k,ispec), &
@@ -756,7 +756,7 @@
c45*duzdxl_plus_duxdzl + c44*duzdyl_plus_duydzl + c34*duzdzl
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! note: fortran passes pointers to array location, thus R_memory(1,1,...) should be fine
call compute_element_att_stress( R_memory(1,1,i,j,k,ispec), &
@@ -1123,7 +1123,7 @@
c45*duzdxl_plus_duxdzl + c44*duzdyl_plus_duydzl + c34*duzdzl
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! note: fortran passes pointers to array location, thus R_memory(1,1,...) should be fine
call compute_element_att_stress( R_memory(1,1,i,j,k,ispec), &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle.f90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle.f90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -691,7 +691,7 @@
endif ! end of test whether isotropic or anisotropic element
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
do i_SLS = 1,N_SLS
R_xx_val = R_memory(1,i_SLS,i,j,k,ispec)
R_yy_val = R_memory(2,i_SLS,i,j,k,ispec)
@@ -912,7 +912,7 @@
! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. )) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. )) then
! use Runge-Kutta scheme to march in time
do i_SLS = 1,N_SLS
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle_Dev.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle_Dev.F90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_crust_mantle_Dev.F90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -624,7 +624,7 @@
! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! updates R_memory
call compute_element_att_memory_cr(ispec,R_memory, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core.f90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core.f90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -420,7 +420,7 @@
endif
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
do i_SLS = 1,N_SLS
R_xx_val = R_memory(1,i_SLS,i,j,k,ispec)
R_yy_val = R_memory(2,i_SLS,i,j,k,ispec)
@@ -644,7 +644,7 @@
! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. )) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. )) then
do i_SLS = 1,N_SLS
factor_common_use = factor_common(i_SLS,:,:,:,ispec)
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core_Dev.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core_Dev.F90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_inner_core_Dev.F90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -539,7 +539,7 @@
endif
! subtract memory variables if attenuation
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! note: fortran passes pointers to array location, thus R_memory(1,1,...) should be fine
call compute_element_att_stress( R_memory(1,1,i,j,k,ispec), &
@@ -805,7 +805,7 @@
! equation (9.59) page 350): Q_\alpha = Q_\mu * 3 * (V_p/V_s)^2 / 4
! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
- if(ATTENUATION_VAL .and. ( USE_ATTENUATION_MIMIC .eqv. .false. ) ) then
+ if(ATTENUATION_VAL .and. ( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) ) then
! updates R_memory
call compute_element_att_memory_ic(ispec,R_memory, &
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/initialize_simulation.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/initialize_simulation.f90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/initialize_simulation.f90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -391,8 +391,8 @@
if((SIMULATION_TYPE == 1 .and. SAVE_FORWARD) .or. SIMULATION_TYPE == 3) then
if ( ATTENUATION_VAL) then
! checks mimic flag:
- ! attenuation for adjoint simulations must have USE_ATTENUATION_MIMIC set by xcreate_header_file
- if( USE_ATTENUATION_MIMIC .eqv. .false. ) &
+ ! attenuation for adjoint simulations must have USE_PHYSICAL_DISPERSION_ONLY set by xcreate_header_file
+ if( USE_PHYSICAL_DISPERSION_ONLY .eqv. .false. ) &
call exit_MPI(myrank,'error in compiled attenuation parameters, please recompile solver 17b')
! user output
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/specfem3D.F90 2013-01-29 01:51:27 UTC (rev 21303)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/specfem3D.F90 2013-01-29 17:56:35 UTC (rev 21304)
@@ -1639,7 +1639,7 @@
if(ATTENUATION_3D_VAL) write(IMAIN,*) 'using 3D attenuation'
- if(USE_ATTENUATION_MIMIC ) write(IMAIN,*) 'mimicking effects on velocity only'
+ if(USE_PHYSICAL_DISPERSION_ONLY) write(IMAIN,*) 'mimicking physical dispersion effects on velocity only'
else
write(IMAIN,*) 'no attenuation'
endif
More information about the CIG-COMMITS
mailing list