From 9cdaedd2c0a07d33d1d1468c848ea8cb7c2d9637 Mon Sep 17 00:00:00 2001
From: Jan at TU-Delft <J.W.Thorbecke@tudelft.nl>
Date: Thu, 13 Sep 2018 14:20:55 +0200
Subject: [PATCH] corecting time-shifts to match kxmod results
---
fdelmodc/defineSource.c | 2 +-
fdelmodc/getParameters.c | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/fdelmodc/defineSource.c b/fdelmodc/defineSource.c
index c1fca21..bc9059c 100644
--- a/fdelmodc/defineSource.c
+++ b/fdelmodc/defineSource.c
@@ -145,7 +145,7 @@ int defineSource(wavPar wav, srcPar src, modPar mod, recPar rec, float **src_nwa
}
if (src.type < 6) { // shift wavelet with +1/2 DeltaT due to staggered in time
- tshift=-(0.5*rec.skipdt+0.5)*wav.dt;
+ tshift=-(0.5*rec.skipdt+1.5)*wav.dt;
for (iw=1;iw<iwmax;iw++) {
om = deltom*iw*tshift;
tmp.r = ctrace[iw].r*cos(-om) - ctrace[iw].i*sin(-om);
diff --git a/fdelmodc/getParameters.c b/fdelmodc/getParameters.c
index 0473ec0..dca4f52 100644
--- a/fdelmodc/getParameters.c
+++ b/fdelmodc/getParameters.c
@@ -174,7 +174,7 @@ int getParameters(modPar *mod, recPar *rec, snaPar *sna, wavPar *wav, srcPar *sr
if (!getparfloat("rec_delay",&rdelay)) rdelay=0.0;
rec->delay=NINT(rdelay/mod->dt);
// mod->tmod += rdelay;
- mod->nt = NINT(mod->tmod/mod->dt)+1;
+ mod->nt = NINT(mod->tmod/mod->dt);
dt = mod->dt;
if (!getparint("src_type",&src->type)) src->type=1;
--
GitLab