2D Pseudo-spectral algorithm of a Multiparameter Wave Equation

This code calculates the acoustic pressure in a two layer squared domain of 200x200m .

There are two media:   

The first one is defined with c0=1510 m/s rho=1 and absor=0.

The second medium is defined according to the variable type:

• type=0: c=c0     rho=1 and abs=0
• type=1: c=2*c0 rho=1 and abs=0
• type=2: c=2*c0 rho=2 and abs=0.0000000000538
   

Detailed with the following scheme:

      _________________________________________________
   1  |                                               |
      |                                               |
      |         c0=1510 rho=1 and abs=0               |
      |                                               |
  100 |-----------------------------------------------|
      | type=0: c=c0 rho=1 and abs=0                  |
      | type=1: c=2*c0 rho=1 and abs=0                |
      | type=2: c=2*c0 rho=2 and abs=0.0000000000538  |
  200 |_______________________________________________|

A Ricker wavelet of f=1MHz is emitted at the position (Nx/2,8).

In all the cases, we also fix the spatial sampling Dx=0.0001m, the temporal step Dt=0.00000001s and the total simulation time 0.00002s.

We use the MULTIPARAMETER WAVE EQUATION formulation to solve this problem.