adds check to compute source derivatives only for CMT sources (in pure adjoint runs)

Author: danielpeter

src/specfem3D/compute_seismograms.f90

@@ -534,7 +534,7 @@ subroutine compute_seismograms_moment_adjoint()
     hprime_xx,hprime_yy,hprime_zz, &
     hxir_store,hetar_store,hgammar_store, &
     hpxir_store,hpetar_store,hpgammar_store, &
-    ELASTIC_SIMULATION
+    ELASTIC_SIMULATION,USE_FORCE_POINT_SOURCE
 
   use specfem_par, only: GPU_MODE, Mesh_pointer
 
@@ -570,6 +570,10 @@ subroutine compute_seismograms_moment_adjoint()
   if (SIMULATION_TYPE /= 2) return
   if (.not. ELASTIC_SIMULATION) return
 
+  ! only for CMT sources
+  if (USE_FORCE_POINT_SOURCE) return
+  if (nrec_local == 0) return
+
   ! strain and moment derivatives are computed here on CPU routines,
   ! thus transfers displacement to the CPU
   if (GPU_MODE) call transfer_displ_from_device(NDIM*NGLOB_AB, displ, Mesh_pointer)

src/specfem3D/save_adjoint_kernels.f90

@@ -1282,13 +1282,15 @@ subroutine save_kernels_source_derivatives()
   integer :: irec_local,ier
   character(len=MAX_STRING_LEN) :: outputname
 
+  ! only for CMT sources
+  if (USE_FORCE_POINT_SOURCE) return
+
   ! checks
-  if (ADIOS_FOR_KERNELS ) stop 'Source derivative kernels not implemented yet for ADIOS'
+  if (ADIOS_FOR_KERNELS) stop 'Source derivative kernels not implemented yet for ADIOS'
 
   ! writes out derivative kernels
   do irec_local = 1, nrec_local
-    write(outputname,'(a,i6.6)') OUTPUT_FILES(1:len_trim(OUTPUT_FILES)) // &
-        '/src_frechet.',number_receiver_global(irec_local)
+    write(outputname,'(a,i6.6)') trim(OUTPUT_FILES) // '/src_frechet.',number_receiver_global(irec_local)
 
     open(unit=IOUT,file=trim(outputname),status='unknown',iostat=ier)
     if (ier /= 0) then