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

#define?pi?3.141592653

main（int?argcchar?*argv[]）
{
void?model_vti（int?nxint?nzint?vnxint?vnzint?ntint?npdfloat?dxfloat?dz
?????????????float?vdxfloat?vdzfloat?favgfloat?tmaxfloat?dtfloat?dtout
?????????????float?pfacchar?FN1[]char?FN2[]char?FN3[]int?ns_sxdint?ds_sxdint?fs_sxdint?zs_sxd
?????????????int?isfloat?**p_calfloat?Circle_Radius
?????????????int?mmint?wtypeint?hsxint?myidfloat?*mu_v）;???????
void?mute_directwave（int?flag_muint?nxint?ntfloat?dtfloat?favg
?????????????????????float?dxfloat?dzint?fs_sxdint?ds_sxdint?zs_sxdint?is
?????????????????????float?mu_vfloat?**p_calint?tt）;????????
	int?ijkisnxnzntvnxvnzi_starti_endmmwtypehsx;
	int?ns_sxdds_sxdfs_sxdzs_sxdnpd;
	float?dxdzvdxvdztmaxdtdtoutpfacfavg;
	int?myidnumprocsCircle_Radiusflag_mu;
??????float?mu_v;
??????float?**p_cal;


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

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

????????nx=601;?????????npd=50;??????tmax=1.5;
	??nz=301;?????????favg=40;?????pfac=1000.0;

?	??dx=5.0;???
????????dz=5.0;???
?????
	??nt=1501;????
????????dt=1.0;
?????
????????ns_sxd=1;???????
????????fs_sxd=300;??????
????????ds_sxd=100;
????????zs_sxd=1;?????


????????Circle_Radius=15;//for?get?rid?of?SV
/*************************v*****************************/
??????vdz=dz;vdx=dx;vnx=nx;vnz=nz;dtout=dt;???
/************************Loop?start**************************/

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

	p_cal=alloc2float（ntnx）;


/*******************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（“--------------------------------------------------------\n“）;
????????printf（“--------------------------------------------------------\n“）;
????????printf（“---???\n“）;
????}??????????????????????????????????????????????????????
/********************IS?Loop?start**************/
???for（is=1+myid;is<=ns_sxd;is=is+numprocs）	
????{?????
??????if（myid==0）
????????{
?????????printf（“---???IS========%d??\n“is）

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

?????文件??????25050??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\alloc.c

?????文件???????4455??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\alloc.h

?????文件???????1826??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\bhdr.h

?????文件??????10918??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\cjbsegy.h

?????文件???????4278??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\complex.c

?????文件???????2310??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\complex.h

?????文件???????1355??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\cwp.h

?????文件?????120423??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\fft.c

?????文件???????1390??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\fft.h

?????文件???????2369??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\hdr.h

?????文件????????587??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\hc\mute_direct.h

?????文件????????113??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\Makefile

?????文件????????218??2018-05-27?13:36??二維VTI介質擬聲波正演模擬與逆時偏移\README.md

?????文件??????28135??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\Toa_model_vti.c

?????文件??????66636??2018-05-27?13:34??二維VTI介質擬聲波正演模擬與逆時偏移\Toa_rtm_vti_adcig_cdp.c

?????目錄??????????0??2018-05-27?17:36??二維VTI介質擬聲波正演模擬與逆時偏移\hc

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

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

???????????????270063????????????????????17