From e30857569c7db5fb02c01f7865f9bd93adf2c5be Mon Sep 17 00:00:00 2001
From: "joeri.brackenhoff" <joeri.brackenhoff@login0.ogun.local>
Date: Wed, 27 Mar 2019 12:37:54 -0300
Subject: [PATCH] 3D
---
fdelmodc3D/getParameters3D.c | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/fdelmodc3D/getParameters3D.c b/fdelmodc3D/getParameters3D.c
index 47f5db5..520ef3b 100644
--- a/fdelmodc3D/getParameters3D.c
+++ b/fdelmodc3D/getParameters3D.c
@@ -218,8 +218,8 @@ long getParameters3D(modPar *mod, recPar *rec, snaPar *sna, wavPar *wav, srcPar
stabfactor = 1.0/sqrt(2.0);
}
else {
- dispfactor = 5;
- stabfactor = 0.606; /* courant number */
+ dispfactor = 6;
+ stabfactor = 0.496; /* courant number */
}
@@ -377,7 +377,7 @@ long getParameters3D(modPar *mod, recPar *rec, snaPar *sna, wavPar *wav, srcPar
}
for (j=0; j<bnd->ntap; j++) {
for (i=0; i<bnd->ntap; i++) {
- for (l=0; l<bnd->ntap; i++) {
+ for (l=0; l<bnd->ntap; l++) {
wfct = (scl*sqrt(i*i+j*j+l*l));
bnd->tapxyz[l*bnd->ntap*bnd->ntap+j*bnd->ntap+i] = exp(-(wfct*wfct));
}
@@ -1047,7 +1047,7 @@ long getParameters3D(modPar *mod, recPar *rec, snaPar *sna, wavPar *wav, srcPar
vmess("*******************************************");
vmess("tsnap1 = %f tsnap2 = %f ", tsnap1, tsnap2);
vmess("dtsnap = %f Nsnap = %li ", dtsnap, sna->nsnap);
- vmess("nzsnap = %li nysnap = %li nxsnap = %li ", sna->nz, sna->nz, sna->nx);
+ vmess("nzsnap = %li nysnap = %li nxsnap = %li ", sna->nz, sna->ny, sna->nx);
vmess("dzsnap = %f dysnap = %f dxsnap = %f ", dzsnap, dysnap, dxsnap);
vmess("zmin = %f zmax = %f ", sub_z0+dz*sna->z1, sub_z0+dz*sna->z2);
vmess("ymin = %f ymax = %f ", sub_y0+dy*sna->y1, sub_y0+dy*sna->y2);
--
GitLab