[cig-commits] r22233 - in seismo/3D/SPECFEM3D_GLOBE/branches/undo_att: . src/specfem3D

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-06-12 09:10:48 -0700 (Wed, 12 Jun 2013)
New Revision: 22233

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90
Log:
set PARTIAL_PHYS_DISPERSION_ONLY = .false. when using UNDO_ATT


Modified: seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90	2013-06-12 15:38:58 UTC (rev 22232)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90	2013-06-12 16:10:48 UTC (rev 22233)
@@ -591,6 +591,7 @@
 !       ! that's to say, the ensemble forward source is kind of a surface force density, not a body force density
 !       ! therefore, we must add it here, before applying the inverse of mass matrix
 !   else
+#ifndef UNDO_ATT
     if ( NOISE_TOMOGRAPHY == 3 ) then
         ! third step of noise tomography, i.e., read the surface movie saved at every timestep
         ! use the movie to reconstruct the ensemble forward wavefield
@@ -602,6 +603,7 @@
                               NSPEC2D_TOP(IREGION_CRUST_MANTLE),noise_surface_movie, &
                               it,jacobian2D_top_crust_mantle,wgllwgll_xy)
     endif
+#endif
 
     ! ****************************************************
     ! **********  add matching with fluid part  **********
@@ -940,6 +942,7 @@
     ! note: this is done here after the Newmark time scheme, otherwise the indexing for sources
     !          and adjoint sources will become more complicated
     !          that is, index it for adjoint sources will match index NSTEP - 1 for backward/reconstructed wavefields
+#ifndef UNDO_ATT
     if(SIMULATION_TYPE == 3 .and. it == 1) then
       call read_forward_arrays(myrank, &
                     b_displ_crust_mantle,b_veloc_crust_mantle,b_accel_crust_mantle, &
@@ -949,6 +952,7 @@
                     b_epsilondev_crust_mantle,b_epsilondev_inner_core, &
                     b_A_array_rotation,b_B_array_rotation,LOCAL_PATH)
     endif
+#endif
 
 ! write the seismograms with time shift
 

Modified: seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90	2013-06-12 15:38:58 UTC (rev 22232)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90	2013-06-12 16:10:48 UTC (rev 22233)
@@ -2254,6 +2254,12 @@
 
   if(SIMULATION_TYPE == 3)then
 
+#ifdef UNDO_ATT
+      if(ATTENUATION_VAL)then
+          PARTIAL_PHYS_DISPERSION_ONLY = .false.
+      endif
+#endif
+
     allocate(displ_crust_mantle_store_as_bwf(NDIM,NGLOB_CRUST_MANTLE,NT_500),stat=ier)
     if( ier /= 0 ) call exit_MPI(myrank,'error allocating displ_crust_mantle_store_as_bwf')
     allocate(displ_outer_core_store_store_as_bwf(NGLOB_OUTER_CORE,NT_500),stat=ier)