[cig-commits] r18762 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Fri Jul 15 07:10:39 PDT 2011
Author: dkomati1
Date: 2011-07-15 07:10:39 -0700 (Fri, 15 Jul 2011)
New Revision: 18762
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
added zero initial condition for the backward field explicitly; it was implicit before in the memory allocation statement
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2011-07-14 22:07:49 UTC (rev 18761)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2011-07-15 14:10:39 UTC (rev 18762)
@@ -2691,7 +2691,7 @@
if(anyabs .and. SIMULATION_TYPE == 2) then
- ! reads in absorbing bounday data
+ ! reads in absorbing boundary data
if(any_elastic) then
call prepare_absorb_elastic(NSTEP,p_sv, &
nspec_left,nspec_right,nspec_bottom,nspec_top, &
@@ -3735,6 +3735,12 @@
allocate(coorg_send_ps_vector_field(d1_coorg_send_ps_vector_field,d2_coorg_send_ps_vector_field))
allocate(coorg_recv_ps_vector_field(d1_coorg_recv_ps_vector_field,d2_coorg_recv_ps_vector_field))
+ if(SIMULATION_TYPE == 2) then ! Adjoint calculation: set zero initial condition for the backward field
+ b_displ_elastic = ZERO
+ b_veloc_elastic = ZERO
+ b_accel_elastic = ZERO
+ endif
+
! *********************************************************
!
! ************* MAIN LOOP OVER THE TIME STEPS *************
More information about the CIG-COMMITS
mailing list