[cig-commits] r22088 - seismo/3D/SPECFEM3D_GLOBE/branches/undo_att

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-05-16 14:40:06 -0700 (Thu, 16 May 2013)
New Revision: 22088

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part1.f90
Log:
added the case SIMULATION_TYPE == 2 to Part 1 (only SIMULATION_TYPE == 3 will be handled by Part 2, since it is the only case in which we need to use the b_* backward arrays)


Modified: seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part1.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part1.f90	2013-05-16 21:34:14 UTC (rev 22087)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part1.f90	2013-05-16 21:40:06 UTC (rev 22088)
@@ -1,4 +1,8 @@
 
+!! DK DK
+!! DK DK this first part handles the cases SIMULATION_TYPE == 1 and SIMULATION_TYPE == 2
+!! DK DK
+
     ! Newmark time scheme update
 
     ! mantle
@@ -909,8 +913,21 @@
                                 islice_selected_source,ispec_selected_source,it, &
                                 hdur,xi_source,eta_source,gamma_source,nu_source)
 
-!   ! add adjoint sources
-!   if (SIMULATION_TYPE == 2 .or. SIMULATION_TYPE == 3) then
+    ! add adjoint sources only if adjoint simulation is performed for source inversion only
+    if (SIMULATION_TYPE == 2) then
+      if( nadj_rec_local > 0 ) &
+        call compute_add_sources_adjoint(myrank,nrec, &
+                                nadj_rec_local,NSTEP,NTSTEP_BETWEEN_READ_ADJSRC, &
+                                accel_crust_mantle,adj_sourcearrays, &
+                                nu,xi_receiver,eta_receiver,gamma_receiver, &
+                                xigll,yigll,zigll,ibool_crust_mantle, &
+                                islice_selected_rec,ispec_selected_rec, &
+                                NSTEP_SUB_ADJ,iadjsrc_len,iadjsrc,iadj_vec, &
+                                it,it_begin,station_name,network_name,DT)
+    endif
+
+!   ! add adjoint sources and add sources for backward/reconstructed wavefield
+!   if (SIMULATION_TYPE == 3) then
 !     if( nadj_rec_local > 0 ) &
 !       call compute_add_sources_adjoint(myrank,nrec, &
 !                               nadj_rec_local,NSTEP,NTSTEP_BETWEEN_READ_ADJSRC, &
@@ -920,15 +937,12 @@
 !                               islice_selected_rec,ispec_selected_rec, &
 !                               NSTEP_SUB_ADJ,iadjsrc_len,iadjsrc,iadj_vec, &
 !                               it,it_begin,station_name,network_name,DT)
-!   endif
-
-!   ! add sources for backward/reconstructed wavefield
-!   if (SIMULATION_TYPE == 3) &
 !     call compute_add_sources_backward(myrank,NSOURCES,NSTEP, &
 !                               b_accel_crust_mantle,sourcearrays, &
 !                               DT,t0,tshift_cmt,hdur_gaussian,ibool_crust_mantle, &
 !                               islice_selected_source,ispec_selected_source,it, &
 !                               hdur,xi_source,eta_source,gamma_source,nu_source)
+!   endif
 
     ! NOISE_TOMOGRAPHY
     if ( NOISE_TOMOGRAPHY == 1 ) then