Skip to content
Snippets Groups Projects
Commit 3f512151 authored by Jan Thorbecke's avatar Jan Thorbecke
Browse files

update to latest getRecTimes

parent 50ee0bfe
No related branches found
No related tags found
No related merge requests found
...@@ -22,9 +22,11 @@ int writeSrcRecPos(modPar *mod, recPar *rec, srcPar *src, shotPar *shot); ...@@ -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 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, int getRecTimes(modPar mod, recPar rec, bndPar bnd, int itime, int isam, float *vx, float *vz, float *tzz, float *txx,
float *rec_p, float *rec_pp, float *rec_ss, int verbose); 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, 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, 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) ...@@ -283,7 +285,7 @@ int main(int argc, char **argv)
if (rec.type.vz) rec_vz = (float *)calloc(size,sizeof(float)); 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.vx) rec_vx = (float *)calloc(size,sizeof(float));
if (rec.type.p) rec_p = (float *)calloc(size,sizeof(float)); if (rec.type.p) rec_p = (float *)calloc(size,sizeof(float));
if(sna.beam) { if(sna.beam) {
size = sna.nz*sna.nx; size = sna.nz*sna.nx;
if (sna.type.vz) beam_vz = (float *)calloc(size,sizeof(float)); 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) ...@@ -377,10 +379,11 @@ shared (shot, bnd, mod, src, wav, rec, ixsrc, izsrc, it, src_nwav, verbose)
isam = (it-rec.delay-itwritten)/rec.skipdt; isam = (it-rec.delay-itwritten)/rec.skipdt;
/* store time at receiver positions */ /* store time at receiver positions */
getRecTimes(mod, rec, bnd, it, isam, vx, vz, tzz, txx, txz, getRecTimes(mod, rec, bnd, it, isam, vx, vz, tzz, txx, txz,
rec_vx, rec_vz, rec_txx, rec_tzz, rec_txz, mu, eprs, ksigma,
rec_p, rec_pp, rec_ss, verbose); 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) */ /* at the end of modeling a shot, write receiver array to output file(s) */
if (writeToFile && (it+rec.skipdt <= it1-1) ) { if (writeToFile && (it+rec.skipdt <= it1-1) ) {
fileno = ( ((it-rec.delay)/rec.skipdt)+1)/rec.nt; fileno = ( ((it-rec.delay)/rec.skipdt)+1)/rec.nt;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment