add sensitivity kernel for love wave

xumi1993 · Nov 26, 2024 · b8e9c34 · b8e9c34
1 parent e412985
commit b8e9c34
Show file tree

Hide file tree

Showing 5 changed files with 967 additions and 6 deletions.
diff --git a/src/model.f90 b/src/model.f90
@@ -234,7 +234,6 @@ subroutine inv1d(this, vsinv, niter, misfits)
       vsinv = vsinv * (1-update_total)
     enddo
     niter = iter
-
   end subroutine inv1d
 
   subroutine add_pert(this,nx,ny,nz,pert_vel,hmarg,anom_size)

diff --git a/src/para.f90 b/src/para.f90
@@ -107,7 +107,7 @@ subroutine read_par_file(this, fname)
         ! if (associated(io_err)) call exit_mpi(myrank, trim(io_err%message))
         this%data%iwave = data_sec%get_integer('iwave', error=io_err)
         if (associated(io_err)) call exit_mpi(myrank, trim(io_err%message))
-        if (this%data%iwave == 2) then
+        if (this%data%iwave == RAYLEIGH) then
           this%data%wave_name = 'rayleigh'
         else
           this%data%wave_name = 'love'

diff --git a/src/shared/shared_par.f90 b/src/shared/shared_par.f90
@@ -4,6 +4,8 @@ module constants
 
   ! if NUMBER_OF_SIMULTANEOUS_RUNS > 1
   character(len=MAX_STRING_LEN) :: OUTPUT_FILES = 'OUTPUT_FILES'
+  integer, parameter :: GROUP_VEL = 1, PHASE_VEL = 0
+  integer, parameter :: RAYLEIGH = 2, LOVE = 1
 
   ! we put a default value here
   real(kind=dp), parameter :: pi = 3.14159265359