二维二阶VTI介质拟声波正演模拟与逆时偏移

二维VTI介质拟声波正演模拟与逆时偏移，二维VTI介质拟声波正演模拟与逆时偏移，二维VTI介质拟声波正演模拟与逆时偏移，二维VTI介质拟声波正演模拟与逆时偏移

//#########################################################
//##         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