From 97b34a1d248fb71e743982c8d5e86734ee9a0e89 Mon Sep 17 00:00:00 2001
From: JanThorbecke <janth@xs4all.nl>
Date: Thu, 9 Jul 2020 14:39:39 +0200
Subject: [PATCH] taper correction
---
marchenko/marchenko.c | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/marchenko/marchenko.c b/marchenko/marchenko.c
index 18f36d2..8bc284a 100644
--- a/marchenko/marchenko.c
+++ b/marchenko/marchenko.c
@@ -283,9 +283,11 @@ int main (int argc, char **argv)
tapersy[j] = (cos(PI*(j-ntap)/ntap)+1)/2.0;
for (j = ntap; j < nxs-ntap; j++)
tapersy[j] = 1.0;
- for (j = nxs-ntap; j < nxs; j++)
- tapersy[j] =(cos(PI*(j-(nxs-ntap))/ntap)+1)/2.0;
+ for (j = nxs-ntap-1; j < nxs; j++)
+ tapersy[j] =(cos(PI*(j+1-(nxs-ntap))/ntap)+1)/2.0;
if (verbose) vmess("Taper for operator applied ntap=%d", ntap);
+ for (j = 0; j < nxs; j++)
+ fprintf(stderr,"taper [%d] = %f \n", j, tapersy[j]);
for (l = 0; l < Nfoc; l++) {
for (i = 0; i < nxs; i++) {
for (j = 0; j < nts; j++) {
@@ -318,8 +320,8 @@ int main (int argc, char **argv)
tapersh[j] = (cos(PI*(j-ntap)/ntap)+1)/2.0;
for (j = ntap; j < nx-ntap; j++)
tapersh[j] = 1.0;
- for (j = nx-ntap; j < nx; j++)
- tapersh[j] =(cos(PI*(j-(nx-ntap))/ntap)+1)/2.0;
+ for (j = nx-ntap-1; j < nx; j++)
+ tapersh[j] =(cos(PI*(j+1-(nx-ntap))/ntap)+1)/2.0;
if (verbose) vmess("Taper for shots applied ntap=%d", ntap);
for (l = 0; l < nshots; l++) {
for (j = 1; j < nw; j++) {
--
GitLab