Merge branch 'master' of https://github.com/ahay/src

Author: sfomel

user/chenyk/Mefd3dmt.c

@@ -423,11 +423,11 @@ int main(int argc, char* argv[])
         if (!sf_histfloat(Fvp,"d3",&dy)) sf_error("No d3= in input");/* veloctiy model: dy */
     	if (!sf_getint("nb",&nb)) nb=30; /* thickness of ABC layer */
     	if (!sf_getint("nt",&nt)) sf_error("nt required");/* number of time steps */
-		if (jsnap>nt) sf_error("make sure jsnap<=nt");
     	if (!sf_getfloat("dt",&dt)) sf_error("dt required");/* time sampling interval */
     	if (!sf_getfloat("fm",&fm)) fm=20.0; /*dominant freq of Ricker wavelet */
    	if (!sf_getint("ft",&ft)) ft=0; /* first recorded time */
     	if (!sf_getint("jsnap",&jsnap)) jsnap=1;	/* interval for snapshots  */
+		if (jsnap>nt) sf_error("make sure jsnap<=nt");
       if(!sf_getfloat("ct",&ct)) ct=0.01;/*for absorbing boundary*/
 
 	if(ifwfd)
@@ -647,24 +647,25 @@ int main(int argc, char* argv[])
 		uvz[sy+nb-1][sx+nb  ][sz+nb+1]-=Mzy*h3/dy/4*wlt[it];*/	
 
 		/*add moment tensor source on stress*/
-		txx[sy+nb  ][sx+nb  ][sz+nb] +=h3*Mxx*wlt[it];
-		tyy[sy+nb  ][sx+nb  ][sz+nb] +=h3*Myy*wlt[it];
-		tzz[sy+nb  ][sx+nb  ][sz+nb] +=h3*Mzz*wlt[it];
+		/*positive and negative +- changed on 01/21/2025*/
+		txx[sy+nb  ][sx+nb  ][sz+nb] -=h3*Mxx*wlt[it];
+		tyy[sy+nb  ][sx+nb  ][sz+nb] -=h3*Myy*wlt[it];
+		tzz[sy+nb  ][sx+nb  ][sz+nb] -=h3*Mzz*wlt[it];
 
-		txy[sy+nb  ][sx+nb  ][sz+nb] +=0.25*h3*Mxy*wlt[it];
-		txy[sy+nb  ][sx+nb-1][sz+nb] +=0.25*h3*Mxy*wlt[it]; 
-		txy[sy+nb-1][sx+nb  ][sz+nb] +=0.25*h3*Mxy*wlt[it];
-		txy[sy+nb-1][sx+nb-1][sz+nb] +=0.25*h3*Mxy*wlt[it];
+		txy[sy+nb  ][sx+nb  ][sz+nb] -=0.25*h3*Mxy*wlt[it];
+		txy[sy+nb  ][sx+nb-1][sz+nb] -=0.25*h3*Mxy*wlt[it]; 
+		txy[sy+nb-1][sx+nb  ][sz+nb] -=0.25*h3*Mxy*wlt[it];
+		txy[sy+nb-1][sx+nb-1][sz+nb] -=0.25*h3*Mxy*wlt[it];
 
-		tyz[sy+nb  ][sx+nb  ][sz+nb ] +=0.25*h3*Myz*wlt[it];
-		tyz[sy+nb  ][sx+nb  ][sz+nb-1]+=0.25*h3*Myz*wlt[it]; 
-		tyz[sy+nb-1][sx+nb  ][sz+nb ] +=0.25*h3*Myz*wlt[it];
-		tyz[sy+nb-1][sx+nb  ][sz+nb-1]+=0.25*h3*Myz*wlt[it];
+		tyz[sy+nb  ][sx+nb  ][sz+nb ] -=0.25*h3*Myz*wlt[it];
+		tyz[sy+nb  ][sx+nb  ][sz+nb-1]-=0.25*h3*Myz*wlt[it]; 
+		tyz[sy+nb-1][sx+nb  ][sz+nb ] -=0.25*h3*Myz*wlt[it];
+		tyz[sy+nb-1][sx+nb  ][sz+nb-1]-=0.25*h3*Myz*wlt[it];
 
-		txz[sy+nb  ][sx+nb  ][sz+nb] +=0.25*h3*Mxz*wlt[it];
-		txz[sy+nb  ][sx+nb-1][sz+nb] +=0.25*h3*Mxz*wlt[it]; 
-		txz[sy+nb  ][sx+nb  ][sz+nb-1] +=0.25*h3*Mxz*wlt[it];
-		txz[sy+nb  ][sx+nb-1][sz+nb-1] +=0.25*h3*Mxz*wlt[it];
+		txz[sy+nb  ][sx+nb  ][sz+nb] -=0.25*h3*Mxz*wlt[it];
+		txz[sy+nb  ][sx+nb-1][sz+nb] -=0.25*h3*Mxz*wlt[it]; 
+		txz[sy+nb  ][sx+nb  ][sz+nb-1] -=0.25*h3*Mxz*wlt[it];
+		txz[sy+nb  ][sx+nb-1][sz+nb-1] -=0.25*h3*Mxz*wlt[it];
 
 		forward_uvy_uvx_uvz(uvx, uvy, uvz, txx, tyy, tzz, txz, txy, tyz, rho);
 		forward_txx_tyy_tzz_txz_txy_tyz(uvx, uvy, uvz, txx, tyy, tzz, txz, txy, tyz, vp, vs);

user/chenyk/SConstruct

@@ -10,7 +10,7 @@ zerotrace hilbert cdivdir cabs cemd1 cemd2 snr2 snr3 snrcyk preerror
 threshold1 tsmf fxdecon simidenoise bin2rsf rsf2bin halfthr hradon
 suplane simidenoise1 rmtrace prekirch psmig pocssemb afd2d
 prekirchsr simpostkirch presr nlm1 nlm2 dblendseis svmf genshotscyk dip_fb dipn_fb dipl dipln selfblend ycvelxf ycvelinvww afd2domp low addtrace matlr pnmomf expd ltfts ltftn fxynpre fxynpre2
-orthol orthon orthonc orthoa smoothn warpscann dipn efd3dmt sint2 sint3 eikonal_rtp eikonal_surf eikonal_surf_dv1d txt2rsf rsf2txt npef5_recon
+orthol orthon orthonc orthoa smoothn warpscann dipn efd3dmt sint2 sint3 eikonal_rtp eikonal_surf eikonal_surf_dv1d txt2rsf rsf2txt npef5_recon aps3d
 '''
 ccprogs = 'diffcxx'