From 3f5121511e0e66562f54936ac60c88f34a0c30e5 Mon Sep 17 00:00:00 2001
From: Jan Thorbecke <janth@xs4all.nl>
Date: Thu, 12 Apr 2018 10:55:42 +0200
Subject: [PATCH] update to latest getRecTimes
---
fdemmodc/fdemmodc.c | 19 +++++++++++--------
1 file changed, 11 insertions(+), 8 deletions(-)
diff --git a/fdemmodc/fdemmodc.c b/fdemmodc/fdemmodc.c
index 8003702..fbd2ffb 100644
--- a/fdemmodc/fdemmodc.c
+++ b/fdemmodc/fdemmodc.c
@@ -22,9 +22,11 @@ int writeSrcRecPos(modPar *mod, recPar *rec, srcPar *src, shotPar *shot);
int em4(modPar mod, srcPar src, wavPar wav, bndPar bnd, int itime, int ixsrc, int izsrc, float **src_nwav, float *hz, float *hx, float *Ey, float *eprs, float *ksigma, float *mu, int verbose);
-int getRecTimes(modPar mod, recPar rec, bndPar bnd, int itime, int isam, float *vx, float *vz, float *tzz, float *txx,
- float *txz, float *rec_vx, float *rec_vz, float *rec_txx, float *rec_tzz, float *rec_txz,
- float *rec_p, float *rec_pp, float *rec_ss, int verbose);
+
+int getRecTimes(modPar mod, recPar rec, bndPar bnd, int itime, int isam, float *vx, float *vz, float *tzz, float *txx,
+ float *txz, float *l2m, float *rox, float *roz,
+ float *rec_vx, float *rec_vz, float *rec_txx, float *rec_tzz, float *rec_txz,
+ float *rec_p, float *rec_pp, float *rec_ss, float *rec_udp, float *rec_udvz, int verbose);
int writeEmRec(recPar rec, modPar mod, int ixsrc, int izsrc, int nsam, int ishot, int fileno,
float *rec_vx, float *rec_vz, float *rec_txx, float *rec_tzz, float *rec_txz,
@@ -283,7 +285,7 @@ int main(int argc, char **argv)
if (rec.type.vz) rec_vz = (float *)calloc(size,sizeof(float));
if (rec.type.vx) rec_vx = (float *)calloc(size,sizeof(float));
if (rec.type.p) rec_p = (float *)calloc(size,sizeof(float));
-
+
if(sna.beam) {
size = sna.nz*sna.nx;
if (sna.type.vz) beam_vz = (float *)calloc(size,sizeof(float));
@@ -377,10 +379,11 @@ shared (shot, bnd, mod, src, wav, rec, ixsrc, izsrc, it, src_nwav, verbose)
isam = (it-rec.delay-itwritten)/rec.skipdt;
/* store time at receiver positions */
- getRecTimes(mod, rec, bnd, it, isam, vx, vz, tzz, txx, txz,
- rec_vx, rec_vz, rec_txx, rec_tzz, rec_txz,
- rec_p, rec_pp, rec_ss, verbose);
-
+ getRecTimes(mod, rec, bnd, it, isam, vx, vz, tzz, txx, txz,
+ mu, eprs, ksigma,
+ rec_vx, rec_vz, rec_txx, rec_tzz, rec_txz,
+ rec_p, rec_pp, rec_ss, NULL, NULL, verbose);
+
/* at the end of modeling a shot, write receiver array to output file(s) */
if (writeToFile && (it+rec.skipdt <= it1-1) ) {
fileno = ( ((it-rec.delay)/rec.skipdt)+1)/rec.nt;
--
GitLab