diff --git a/fdemmodc/demo/em.scr b/fdemmodc/demo/em.scr index 441956ffa2e97a4faf858e7294d42ed2722d7e41..058a81a6c2a7dac25477912d16e63aaa7a86af15 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 98673278a357ea891e0035e204c0008cc711ee57..fb6e6475bc17aa5fb373ff6fef485ed666113015 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]); } }