diff --git a/fdelmodc/defineSource.c b/fdelmodc/defineSource.c index c1fca2140b5195e749e36a4aa21ddb7cbdbd2ccd..bc9059c73b518980d35a4b3356796ddaa879b4c4 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 0473ec065a7690fdf0c4ee6e4c227178840026db..dca4f5248913ee2d3458554878fe86becb39e6c8 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;