二维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 “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