diff --git a/fdelmodc/fdelmodc.c b/fdelmodc/fdelmodc.c index ae1e91baa54c5ba4eece610286a1c13629ba7024..17fcb9f9bff0864098b02dc9f63bb33429429a80 100644 --- a/fdelmodc/fdelmodc.c +++ b/fdelmodc/fdelmodc.c @@ -599,8 +599,8 @@ shared (shot, bnd, mod, src, wav, rec, ixsrc, izsrc, it, src_nwav, verbose) writeToFile = ! ( (((it-rec.delay)/rec.skipdt)+1)%rec.nt ); itwritten = fileno*(rec.nt)*rec.skipdt; - isam = (it-rec.delay-itwritten)/rec.skipdt+1; /* Note that time step it=0 (t=0 for t**-fields t=-1/2 dt for v*-field) is not recorded */ + isam = (it-rec.delay-itwritten)/rec.skipdt+1; /* store time at receiver positions */ getRecTimes(mod, rec, bnd, it, isam, vx, vz, tzz, txx, txz, l2m, rox, roz, diff --git a/fdelmodc/getParameters.c b/fdelmodc/getParameters.c index de3b10bf1f2adf60f02d68eacefcdb0e034bf047..34244b3f624e1522dbfe375ab86330ef263553ff 100644 --- a/fdelmodc/getParameters.c +++ b/fdelmodc/getParameters.c @@ -1118,7 +1118,7 @@ int getParameters(modPar *mod, recPar *rec, snaPar *sna, wavPar *wav, srcPar *sr if (!getparint("scale",&rec->scale)) rec->scale=0; if (!getparfloat("dxspread",&dxspread)) dxspread=0; if (!getparfloat("dzspread",&dzspread)) dzspread=0; - rec->nt=MIN(rec->nt, NINT((mod->tmod-rdelay)/dtrcv)); + rec->nt=MIN(rec->nt, NINT((mod->tmod-rdelay)/dtrcv)+1); /* allocation of receiver arrays is done in recvPar */ /*