made fdemmodc up-to-date with latest source changes in fdelmodc

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 \

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]);
 		}
 	}