Merge branch 'master' of https://github.com/JanThorbecke/OpenSource

marchenko/marchenko.c

@@ -284,7 +284,7 @@ int main (int argc, char **argv)
         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;
+            tapersy[j] =(cos(PI*(j+1-(nxs-ntap))/ntap)+1)/2.0;
         if (verbose) vmess("Taper for operator applied ntap=%d", ntap);
         for (l = 0; l < Nfoc; l++) {
             for (i = 0; i < nxs; i++) {
@@ -319,7 +319,7 @@ int main (int argc, char **argv)
         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;
+            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++) {