//#########################################################
//##?????????2D?Acoustic?VTI?Medium?Forward???
//##??Ps?:?P?+?sv?wave?and?get?rid?of?sv????????
//##?????????????????????????????????????Rong?Tao?
//########################################################
#include
#include
#include
#include
#include
#include?“mpi.h“
#include?“/home/Toa/hc/cjbsegy.h“
#include?“/home/Toa/hc/fft.c“
#include?“/home/Toa/hc/alloc.c“
#include?“/home/Toa/hc/complex.c“

#define?pi?3.141592653
float?stencil[5]={-205.0/72.08.0/5.0-1.0/5.08.0/315.0-1.0/560.0};
main（int?argcchar?*argv[]）
{
void?model_vti（int?nxint?nzint?ntint?npdfloat?dxfloat?dzfloat?favgfloat?tmaxfloat?dt
?????????????float?pfacchar?FN1[]char?FN2[]char?FN3[]int?nsint?dsint?fsint?zs
?????????????int?isfloat?**p_calfloat?**P_upfloat?**P_dofloat?**P_lefloat?**P_rifloat?**Q_upfloat?**Q_dofloat?**Q_lefloat?**Q_ri
?????????????float?Circle_Radiusint?mmint?wtypeint?hsxint?myidfloat?*mu_vfloat?**Pfloat?**Q）;???????
void?mute_directwave（int?flag_muint?nxint?ntfloat?dtfloat?favg
?????????????????????float?dxfloat?dzint?fsint?dsint?zsint?is
?????????????????????float?mu_vfloat?**p_calint?tt）;????????
	int?ijkisnxnznti_starti_endmmwtypehsx;
	int?nsdsfszsnpd;
	float?dxdztmaxdtpfacfavg;
	int?myidnumprocsCircle_Radiusflag_mu;
????????float?mu_v;
????????float?**p_cal**P_up**P_do**P_le**P_ri**Q_up**Q_do**Q_le**Q_ri;
????????float?**P**Q;


/********?rankswaveletreceiversmute?direct?**********/
????????mm=4;wtype=1;hsx=1;flag_mu=0;
/*******************?dat?document?**********************/
????????char?FN1[250]={“vel100100.dat“};
????????char?FN2[250]={“epsilu100100.dat“};
????????char?FN3[250]={“deta100100.dat“};
	????char?FN4[250]={“shot.dat“};

/*****************?parameters?**************************/

????????nx=100;?????????npd=50;??????tmax=1.5;
	?nz=100;?????????favg=20;?????pfac=1000.0;

?	?dx=5.0;???
????????dz=5.0;???
??????
	?nt=101;????
????????dt=1.0;
?????
????????ns=1;???????
????????fs=nx/ns/2;??????
????????ds=nx/ns;
????????zs=1;?????


????????Circle_Radius=15;//for?get?rid?of?SV

/************************Loop?start**************************/

??????FILE?*fp4;
??????fp4=fopen（FN4“wb“）;

	P=alloc2float（nznx）;	????????Q=alloc2float（nznx）;

	p_cal=alloc2float（ntnx）;
	P_up=alloc2float（ntnx）;	Q_up=alloc2float（ntnx）;
	P_do=alloc2float（ntnx）;	Q_do=alloc2float（ntnx）;
	P_le=alloc2float（ntnz）;	Q_le=alloc2float（ntnz）;
	P_ri=alloc2float（ntnz）;	Q_ri=alloc2float（ntnz）;


/*******************MPI************************/
??????MPI_Init（&argc&argv）;
??????MPI_Comm_rank（MPI_COMM_WORLD&myid）;
??????MPI_Comm_size（MPI_COMM_WORLD&numprocs）;
/*******************MPI***********************/	
??????MPI_Barrier（MPI_COMM_WORLD）;
???if（myid==0）
????{
????????printf（“-----------------------------------

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????17829??2017-10-26?11:49??二維VTI介質擬聲波正演模擬與逆時偏移\Toa_model_VTI_2order.c

?????文件??????28832??2017-02-18?15:18??二維VTI介質擬聲波正演模擬與逆時偏移\Toa_RTM_VTI_2order.c

?????目錄??????????0??2018-05-29?11:13??二維VTI介質擬聲波正演模擬與逆時偏移

-----------?---------??----------?-----??----

????????????????46661????????????????????3