- added reuse option for external source time function file to CMTSOL…

…UTION
SPECFEM · Nov 22, 2023 · b3bf93e · b3bf93e
1 parent 3eb4e22
commit b3bf93e
Showing 1 changed file with 17 additions and 3 deletions.
diff --git a/src/specfem3D/get_cmt.f90 b/src/specfem3D/get_cmt.f90
@@ -55,6 +55,7 @@ subroutine get_cmt(CMTSOLUTION,yr,jda,mo,da,ho,mi,sec, &
 
   ! local variables below
   integer :: julian_day,isource
+  integer :: ishift,nright
   integer :: i,itype,istart,iend,ier,ipos
   double precision :: t_shift(NSOURCES)
   character(len=256) :: string
@@ -372,9 +373,22 @@ subroutine get_cmt(CMTSOLUTION,yr,jda,mo,da,ho,mi,sec, &
       ! gets external STF file name
       read(IIN,"(a)") string
       external_source_time_function_filename = trim(string)
-
-      ! reads in stf values
-      call read_external_source_time_function(isource,user_source_time_function,external_source_time_function_filename)
+      if (trim(external_source_time_function_filename) == 'REUSE' .or. &
+         trim(external_source_time_function_filename) == 'Reuse' .or. &
+         trim(external_source_time_function_filename) == 'reuse') then
+        ! Reuse the source time function of the first source.
+        if (isource == 1) then
+          stop 'Error: "reuse" option cannot be used for the first source!'
+        endif
+        ishift=nint(t_shift(isource)/DT)
+        nright=NSTEP_STF-ishift
+        user_source_time_function(1:ishift,isource)=0.0_CUSTOM_REAL
+        user_source_time_function(ishift+1:NSTEP_STF,isource)=user_source_time_function(1:nright,1)
+      else
+
+        ! reads in stf values
+        call read_external_source_time_function(isource,user_source_time_function,external_source_time_function_filename)
+      endif
     endif
 
   enddo