From 151ee82ba555069a42993d7c0adb1187839003a7 Mon Sep 17 00:00:00 2001
From: Jan Thorbecke <janth@xs4all.nl>
Date: Thu, 3 May 2018 16:21:39 +0200
Subject: [PATCH] made fdemmodc up-to-date with latest source changes in
fdelmodc
---
fdemmodc/demo/em.scr | 2 +-
fdemmodc/em4.c | 6 +++---
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/fdemmodc/demo/em.scr b/fdemmodc/demo/em.scr
index 441956f..058a81a 100755
--- a/fdemmodc/demo/em.scr
+++ b/fdemmodc/demo/em.scr
@@ -45,7 +45,7 @@ makewave fp=250e+6 dt=$dt file_out=wave.su nt=4096 t0=6e-9 verbose=1
src_injectionrate=1 \
dtrcv=1e-10 \
rec_delay=6e-9 \
- verbose=6 \
+ verbose=3 \
tmod=2e-7 \
dt=$dt \
dxrcv=0.02 \
diff --git a/fdemmodc/em4.c b/fdemmodc/em4.c
index 9867327..fb6e647 100644
--- a/fdemmodc/em4.c
+++ b/fdemmodc/em4.c
@@ -92,7 +92,7 @@ int em4(modPar mod, srcPar src, wavPar wav, bndPar bnd, int itime, int ixsrc, in
c1*(Ey[ix*n1+iz] - Ey[(ix-1)*n1+iz]) +
c2*(Ey[(ix+1)*n1+iz] - Ey[(ix-2)*n1+iz]));
//if (hz[ix*n1+iz] > 0.1*FLT_MAX) fprintf(stderr,"%d: hz[%d %d] = %e\n", itime, ix, iz, hz[ix*n1+iz]);
- if (hz[ix*n1+iz] != 0.0) fprintf(stderr,"%d: hz[%d %d] = %e\n", itime, ix, iz, hz[ix*n1+iz]);
+// if (hz[ix*n1+iz] != 0.0) fprintf(stderr,"%d: hz[%d %d] = %e\n", itime, ix, iz, hz[ix*n1+iz]);
}
}
@@ -105,7 +105,7 @@ int em4(modPar mod, srcPar src, wavPar wav, bndPar bnd, int itime, int ixsrc, in
c1*(Ey[ix*n1+iz] - Ey[ix*n1+iz-1]) +
c2*(Ey[ix*n1+iz+1] - Ey[ix*n1+iz-2]));
//if (hx[ix*n1+iz] > 0.1*FLT_MAX) fprintf(stderr,"%d: hx[%d %d] = %e\n", itime, ix, iz, hx[ix*n1+iz]);
- if (fabs(hx[ix*n1+iz]) != 0.0) fprintf(stderr,"%d: hx[%d %d] = %e\n", itime, ix, iz, hx[ix*n1+iz]);
+// if (fabs(hx[ix*n1+iz]) != 0.0) fprintf(stderr,"%d: hx[%d %d] = %e\n", itime, ix, iz, hx[ix*n1+iz]);
}
}
@@ -140,7 +140,7 @@ int em4(modPar mod, srcPar src, wavPar wav, bndPar bnd, int itime, int ixsrc, in
for (iz=mod.ioPz; iz<mod.iePz; iz++) {
Ey[ix*n1+iz] -= eprs[ix*n1+iz]*(dzhx[iz]+dxhz[iz]) + ksigma[ix*n1+iz]*Ey[ix*n1+iz];
//if (Ey[ix*n1+iz] > 0.1*FLT_MAX) fprintf(stderr,"%d: Ey[%d %d] = %e\n", itime, ix, iz, Ey[ix*n1+iz]);
- if (Ey[ix*n1+iz] != 0.0) fprintf(stderr,"%d: Ey[%d %d] = %e\n", itime, ix, iz, Ey[ix*n1+iz]);
+// if (Ey[ix*n1+iz] != 0.0) fprintf(stderr,"%d: Ey[%d %d] = %e\n", itime, ix, iz, Ey[ix*n1+iz]);
}
}
--
GitLab