From a395c0fa43522542c3db992670aa863db21480d6 Mon Sep 17 00:00:00 2001
From: Jan at TU-Delft <J.W.Thorbecke@tudelft.nl>
Date: Thu, 22 Nov 2018 07:39:58 +0100
Subject: [PATCH] snaphot timing
---
fdelmodc/fdelmodc.c | 2 +-
fdelmodc/getParameters.c | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/fdelmodc/fdelmodc.c b/fdelmodc/fdelmodc.c
index ae1e91b..17fcb9f 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 de3b10b..34244b3 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 */
/*
--
GitLab