[cig-commits] r12698 - seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2008-08-21 16:48:07 -0700 (Thu, 21 Aug 2008)
New Revision: 12698

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/read_compute_parameters.F90
   seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
Log:
impose a few parameters and save seismograms every 10,000 time steps in the case of Gordon Bell runs


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/read_compute_parameters.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/read_compute_parameters.F90	2008-08-21 23:25:16 UTC (rev 12697)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/read_compute_parameters.F90	2008-08-21 23:48:07 UTC (rev 12698)
@@ -1111,7 +1111,19 @@
   call read_value_logical(PRINT_SOURCE_TIME_FUNCTION, 'solver.PRINT_SOURCE_TIME_FUNCTION')
 
 !! DK DK temporary patch for the large Gordon Bell runs
-  if(PATCH_FOR_GORDON_BELL) RECEIVERS_CAN_BE_BURIED = .false.
+  if(PATCH_FOR_GORDON_BELL) then
+    NTSTEP_BETWEEN_OUTPUT_INFO      = 100
+    NTSTEP_BETWEEN_OUTPUT_SEISMOS   = 5000000
+    NTSTEP_BETWEEN_READ_ADJSRC      = 1000
+    OUTPUT_SEISMOS_ASCII_TEXT       = .true.
+    OUTPUT_SEISMOS_SAC_ALPHANUM     = .false.
+    OUTPUT_SEISMOS_SAC_BINARY       = .false.
+    ROTATE_SEISMOGRAMS_RT           = .false.
+    WRITE_SEISMOGRAMS_BY_MASTER     = .true.
+    SAVE_ALL_SEISMOS_IN_ONE_FILE    = .true.
+    USE_BINARY_FOR_LARGE_FILE       = .false.
+    RECEIVERS_CAN_BE_BURIED         = .false.
+  endif
 
   if(err_occurred() /= 0) stop 'an error occurred while reading the parameter file'
 

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90	2008-08-21 23:25:16 UTC (rev 12697)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/specfem3D.F90	2008-08-21 23:48:07 UTC (rev 12698)
@@ -2511,6 +2511,20 @@
 
   endif ! nrec_local
 
+!! DK DK added this for Gordon Bell runs, to save seismograms every 10000 time steps
+!! DK DK just in case the full simulation does not finish
+  if(PATCH_FOR_GORDON_BELL .and. mod(it,10000) == 0) then
+      call write_seismograms(myrank,uxdstore,uydstore,uzdstore,number_receiver_global,station_name, &
+            network_name,stlat,stlon,stele,nrec,nrec_local,DT,t0,it_end, &
+            yr_SAC,jda_SAC,ho_SAC,mi_SAC,sec_SAC,t_cmt_SAC, &
+            elat_SAC,elon_SAC,depth_SAC,mb_SAC,ename_SAC,cmt_lat_SAC,cmt_lon_SAC,&
+            cmt_depth_SAC,cmt_hdur_SAC,NSOURCES_SAC,NPROCTOT, &
+            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,one_seismogram,scale_displ,nu)
+  endif
+
 ! write the current or final seismograms
   if(COMPUTE_STORE_SEISMOGRAMS .and. (seismo_current == NTSTEP_BETWEEN_OUTPUT_SEISMOS .or. it == it_end)) then