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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-06-13 10:11:36 -0700 (Thu, 13 Jun 2013)
New Revision: 22251

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att_changed_from_forward.F90
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90
Log:
try to fix the time error 


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-13 16:54:36 UTC (rev 22250)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90	2013-06-13 17:11:36 UTC (rev 22251)
@@ -1012,6 +1012,17 @@
   ! write the current or final seismograms
   if(seismo_current == NTSTEP_BETWEEN_OUTPUT_SEISMOS .or. it == it_end) then
     if (SIMULATION_TYPE == 1 .or. SIMULATION_TYPE == 3) then
+      do irec_local = 1,nrec_local
+        do i = 1,seismo_current/NT_500
+           do j = 1,NT_500/2
+              do k = 1,3
+                seismograms_temp(k) = seismograms(k,irec_local,(i-1)*NT_500 + j)
+                seismograms(k,irec_local,(i-1)*NT_500 + j)  = seismograms(k,irec_local,(i-1)*NT_500 + (NT_500-j+1))
+                seismograms(k,irec_local,(i-1)*NT_500 + (NT_500-j+1)) = seismograms_temp(k)  
+              enddo          
+           enddo
+        enddo
+      enddo
       call write_seismograms(myrank,seismograms,number_receiver_global,station_name, &
             network_name,stlat,stlon,stele,stbur, &
             nrec,nrec_local,ANGULAR_WIDTH_XI_IN_DEGREES,NEX_XI,DT,t0,it_end, &

Modified: seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att_changed_from_forward.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att_changed_from_forward.F90	2013-06-13 16:54:36 UTC (rev 22250)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att_changed_from_forward.F90	2013-06-13 17:11:36 UTC (rev 22251)
@@ -795,6 +795,17 @@
   ! write the current or final seismograms
 if(UNDO_ATT_WITH_STORE)then
   if(seismo_current == NTSTEP_BETWEEN_OUTPUT_SEISMOS .or. it == it_end) then
+    do irec_local = 1,nrec_local
+      do i = 1,seismo_current/NT_500
+         do j = 1,NT_500/2
+            do k = 1,3
+              seismograms_temp(k) = seismograms(k,irec_local,(i-1)*NT_500 + j)
+              seismograms(k,irec_local,(i-1)*NT_500 + j)  = seismograms(k,irec_local,(i-1)*NT_500 + (NT_500-j+1))
+              seismograms(k,irec_local,(i-1)*NT_500 + (NT_500-j+1)) = seismograms_temp(k)  
+            enddo          
+         enddo
+      enddo
+    enddo
     if (SIMULATION_TYPE == 3) then
       call write_seismograms(myrank,seismograms,number_receiver_global,station_name, &
             network_name,stlat,stlon,stele,stbur, &

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-13 16:54:36 UTC (rev 22250)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90	2013-06-13 17:11:36 UTC (rev 22251)
@@ -922,9 +922,10 @@
   real(kind=CUSTOM_REAL), dimension(:,:), allocatable :: displ_outer_core_store_store_as_bwf
   real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: displ_inner_core_store_as_bwf
 
-  integer :: iteration_on_subset,it_of_this_subset,j
+  integer :: iteration_on_subset,it_of_this_subset,j,irec_local,k
   integer :: it_temp,seismo_current_temp
   logical :: UNDO_ATT_WITH_STORE
+  real(kind=CUSTOM_REAL), dimension(3) :: seismograms_temp
 
   include "declaration_part_for_backward_wavefield_simulation.f90"