sanity check for inverse problem

Author: nqdu

src/shared/read_parameter_file.F90

@@ -849,6 +849,16 @@ subroutine read_parameter_file(BROADCAST_AFTER_READ)
       STACEY_INSTEAD_OF_FREE_SURFACE = .false.
     endif
 
+    ! inverse problem is enabled, disable USE_CMT_AND_FORCE/USE_BINARY_SOURCE_FILE
+    if(INVERSE_FWI_FULL_PROBLEM .and. (USE_CMT_AND_FORCE_SOURCE)) then
+      if(myrank == 0) then
+        write(*,'(a)') 'INVERSE_FWI_FULL_PROBLEM is enabled, disabling USE_CMT_AND_FORCE_SOURCE'
+        write(*,*)
+      endif
+      USE_CMT_AND_FORCE_SOURCE = .false.
+      USE_BINARY_SOURCE_FILE = .false.
+    endif
+
     ! checks parameter settings
     call check_simulation_parameters()
 

src/specfem3D/setup_sources_receivers.f90

@@ -465,6 +465,18 @@ subroutine setup_stf_constants()
   ! convert the half duration for triangle STF to the one for Gaussian STF
   hdur_Gaussian(:) = hdur(:) / SOURCE_DECAY_MIMIC_TRIANGLE
 
+  ! sanity check
+  if(INVERSE_FWI_FULL_PROBLEM) then 
+    if(.not. allocated(is_POINTFORCE)) then 
+      allocate(is_POINTFORCE(NSOURCES))
+      if(USE_FORCE_POINT_SOURCE) then 
+        is_POINTFORCE(:) = .true.
+      else 
+        is_POINTFORCE(:) = .false.
+      endif
+    endif
+  endif
+
   !loop every force source to find t0
   t0 = 0. 
   t0_cmt = 0