[cig-commits] r18762 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[dkomati1 at geodynamics.org]

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 *************