From 8f04f5169368de25c89275eac7c2615d50b417ef Mon Sep 17 00:00:00 2001
From: JanThorbecke <janth@xs4all.nl>
Date: Thu, 28 Jan 2021 09:33:35 +0100
Subject: [PATCH] more verbose
---
marchenko/marchenko.c | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/marchenko/marchenko.c b/marchenko/marchenko.c
index 1ed7503..1021ecf 100644
--- a/marchenko/marchenko.c
+++ b/marchenko/marchenko.c
@@ -272,14 +272,16 @@ int main (int argc, char **argv)
for (i=0; i<nxs; i++) muteW[i] = muteW[i]-itmin;
timeShift(G_d, nts, nxs, dt, tshift, fmin, fmax);
}
-/*
-*/
for (i=0; i<nxs; i++) tsynW[i] = NINT(i*dxs*p/dt);
//for (i=0; i<nxs; i++) tsynW[i] = 0.0;
if (Nfoc!=1) verr("For plane-wave focusing only one function can be computed at the same time");
//fprintf(stderr,"itmin=%d tshift=%f =%d \n", itmin, tshift, NINT(tshift/dt));
//for (i=0; i<nxs; i++) fprintf(stderr,"i=%d window=%f\n", i, tsynW[i]*dt);
+/* // TESTING SHIFT
+ tshift=0.3;
+ for (i=0; i<nxs; i++) tsynW[i] = NINT(0.3/dt);
+*/
}
else { /* just fill with zero's */
//itmin=0;
@@ -618,11 +620,14 @@ int main (int argc, char **argv)
/* Apply mute with window for Gmin */
if ( plane_wave==1 ) {
/* for plane wave with angle tshift downward */
- if (verbose>1) vmess("Gmin planewave tshift=%f", tshift);
applyMute_tshift(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, 1, tsynW);
if (src_angle > 0.0) {
+ if (verbose>1) vmess("Gmin planewave tshift=%f", tshift);
timeShift(Gmin, nts, npos, dt, tshift, fmin, fmax);
}
+ else {
+ if (verbose>1) vmess("Gmin NO planewave tshift");
+ }
}
else {
applyMute(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, tsynW);
--
GitLab