fixes kernel runs w/out saving seismos

danielpeter · Dec 20, 2023 · 37c8eb5 · 37c8eb5
1 parent f4ef86b
commit 37c8eb5
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 5 deletions.
diff --git a/src/specfem3D/setup_sources_receivers.f90 b/src/specfem3D/setup_sources_receivers.f90
@@ -1247,9 +1247,6 @@ subroutine setup_receivers()
     do_save_seismograms = .true.
   endif
 
-  ! sets local receivers to zero if no seismogram needs to be saved
-  if (.not. do_save_seismograms) nrec_local = 0
-
   ! seismogram array length (to write out time portions of the full seismograms)
   nlength_seismogram = NTSTEP_BETWEEN_OUTPUT_SEISMOS / NTSTEP_BETWEEN_OUTPUT_SAMPLE
 

diff --git a/src/specfem3D/write_seismograms.f90 b/src/specfem3D/write_seismograms.f90
@@ -40,7 +40,7 @@ subroutine write_seismograms()
     NTSTEP_BETWEEN_OUTPUT_SEISMOS,NTSTEP_BETWEEN_OUTPUT_SAMPLE, &
     do_save_seismograms, &
     WRITE_SEISMOGRAMS_BY_MAIN,OUTPUT_SEISMOS_ASDF, &
-    SAVE_SEISMOGRAMS_IN_ADJOINT_RUN,SAVE_SEISMOGRAMS_STRAIN, &
+    SAVE_SEISMOGRAMS_STRAIN, &
     moment_der,sloc_der,shdur_der,stshift_der, &
     scale_displ
 
@@ -173,7 +173,7 @@ subroutine write_seismograms()
       select case (SIMULATION_TYPE)
       case (1,3)
         ! forward/reconstructed wavefields
-        if (.not. ( SIMULATION_TYPE == 3 .and. (.not. SAVE_SEISMOGRAMS_IN_ADJOINT_RUN) ) ) &
+        if (do_save_seismograms) &
           call write_seismograms_to_file()
         if (SAVE_SEISMOGRAMS_STRAIN) &
           call write_seismograms_strain()