I figured that the segmentation fault in convmPy.f90 occurred in the loop computing the seismograms in the time domain. It seems that it's just that you forgot to correct for the python vs fortran index convention when using rindex (which starts at 0) to index fortran arrays.
I added one line to explicitly show the end of the OMP DO ORDERED.
…and Fortran different indexing
I figured that the segmentation fault in convmPy.f90 occurred in the loop computing the seismograms in the time domain. It seems that it's just that you forgot to correct for the python vs fortran index convention when using
rindex
(which starts at 0) to index fortran arrays.I added one line to explicitly show the end of the OMP DO ORDERED.