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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-06-12 08:38:58 -0700 (Wed, 12 Jun 2013)
New Revision: 22232

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90
Log:
remove write seismogram for backward wave field simulation


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 14:51:13 UTC (rev 22231)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/part2_undo_att.F90	2013-06-12 15:38:58 UTC (rev 22232)
@@ -991,33 +991,33 @@
   endif ! nrec_local
 
   ! 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
-      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, &
-            yr_SAC,jda_SAC,ho_SAC,mi_SAC,sec_SAC,t_cmt_SAC,t_shift_SAC, &
-            elat_SAC,elon_SAC,depth_SAC,event_name_SAC,cmt_lat_SAC,cmt_lon_SAC,&
-            cmt_depth_SAC,cmt_hdur_SAC,NPROCTOT_VAL, &
-            OUTPUT_SEISMOS_ASCII_TEXT,OUTPUT_SEISMOS_SAC_ALPHANUM, &
-            OUTPUT_SEISMOS_SAC_BINARY,ROTATE_SEISMOGRAMS_RT,NTSTEP_BETWEEN_OUTPUT_SEISMOS, &
-            seismo_offset,seismo_current,WRITE_SEISMOGRAMS_BY_MASTER, &
-            SAVE_ALL_SEISMOS_IN_ONE_FILE,USE_BINARY_FOR_LARGE_FILE,MODEL)
-      if(myrank==0) then
-        write(IMAIN,*)
-        write(IMAIN,*) ' Total number of time steps written: ', it-it_begin+1
-        write(IMAIN,*)
-      endif
-!   else ! case of SIMULATION_TYPE == 2
-!     if( nrec_local > 0 ) &
-!       call write_adj_seismograms(seismograms,number_receiver_global, &
-!                                 nrec_local,it,nit_written,DT, &
-!                                 NSTEP,NTSTEP_BETWEEN_OUTPUT_SEISMOS,t0,LOCAL_PATH)
-!       nit_written = it
-    endif
-    seismo_offset = seismo_offset + seismo_current
-    seismo_current = 0
-  endif
+!ZN  if(seismo_current == NTSTEP_BETWEEN_OUTPUT_SEISMOS .or. it == it_end) then
+!ZN    if (SIMULATION_TYPE == 1 .or. SIMULATION_TYPE == 3) then
+!ZN      call write_seismograms(myrank,seismograms,number_receiver_global,station_name, &
+!ZN            network_name,stlat,stlon,stele,stbur, &
+!ZN            nrec,nrec_local,ANGULAR_WIDTH_XI_IN_DEGREES,NEX_XI,DT,t0,it_end, &
+!ZN            yr_SAC,jda_SAC,ho_SAC,mi_SAC,sec_SAC,t_cmt_SAC,t_shift_SAC, &
+!ZN            elat_SAC,elon_SAC,depth_SAC,event_name_SAC,cmt_lat_SAC,cmt_lon_SAC,&
+!ZN            cmt_depth_SAC,cmt_hdur_SAC,NPROCTOT_VAL, &
+!ZN            OUTPUT_SEISMOS_ASCII_TEXT,OUTPUT_SEISMOS_SAC_ALPHANUM, &
+!ZN            OUTPUT_SEISMOS_SAC_BINARY,ROTATE_SEISMOGRAMS_RT,NTSTEP_BETWEEN_OUTPUT_SEISMOS, &
+!ZN            seismo_offset,seismo_current,WRITE_SEISMOGRAMS_BY_MASTER, &
+!ZN            SAVE_ALL_SEISMOS_IN_ONE_FILE,USE_BINARY_FOR_LARGE_FILE,MODEL)
+!ZN      if(myrank==0) then
+!ZN        write(IMAIN,*)
+!ZN        write(IMAIN,*) ' Total number of time steps written: ', it-it_begin+1
+!ZN        write(IMAIN,*)
+!ZN      endif
+!ZN!   else ! case of SIMULATION_TYPE == 2
+!ZN!     if( nrec_local > 0 ) &
+!ZN!       call write_adj_seismograms(seismograms,number_receiver_global, &
+!ZN!                                 nrec_local,it,nit_written,DT, &
+!ZN!                                 NSTEP,NTSTEP_BETWEEN_OUTPUT_SEISMOS,t0,LOCAL_PATH)
+!ZN!       nit_written = it
+!ZN    endif
+!ZN    seismo_offset = seismo_offset + seismo_current
+!ZN    seismo_current = 0
+!ZN  endif
 
 !
 !-------------------------------------------------------------------------------------------------

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 14:51:13 UTC (rev 22231)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/specfem3D.F90	2013-06-12 15:38:58 UTC (rev 22232)
@@ -2203,8 +2203,8 @@
   if(SIMULATION_TYPE == 1)then
     it = 0
     do iteration_on_subset = 1, NSTEP / NT_500
-
-      call save_forward_arrays_undoatt(myrank,SIMULATION_TYPE,SAVE_FORWARD, &
+      if(SAVE_FORWARD)then
+        call save_forward_arrays_undoatt(myrank,SIMULATION_TYPE,SAVE_FORWARD, &
                     NUMBER_OF_RUNS,NUMBER_OF_THIS_RUN, &
                     displ_crust_mantle,veloc_crust_mantle,accel_crust_mantle, &
                     displ_inner_core,veloc_inner_core,accel_inner_core, &
@@ -2212,6 +2212,7 @@
                     R_memory_crust_mantle,R_memory_inner_core, &
                     epsilondev_crust_mantle,epsilondev_inner_core, &
                     A_array_rotation,B_array_rotation,LOCAL_PATH,iteration_on_subset)
+      endif
 
       do it_of_this_subset = 1, NT_500