search

seismology / mc_kernel

Calculate seismic sensitivity kernels using Axisem
GNU General Public License v3.0
24 stars 8 forks source link

k_a in kernel.f90 #60

Closed kasra-hosseini closed 6 years ago

kasra-hosseini commented 6 years ago

Alex spotted the following mistake:

https://github.com/seismology/mc_kernel/blob/9d2981da26d05f1ac786ed07cba8301e28705c59/src/kernel.f90#L556-L558

Compared to Fichtner's book, p.168 eq 9.20d, it should be:

conv_field_fd = ( bw_field_fd(:,2,:) + bw_field_fd(:,1,:) ) * &
                          ( fw_field_fd(:,2,:) + fw_field_fd(:,1,:) )

Does this make sense?