* Data for Dt, Duration, and Period *
    0.15      64 : Delta Time (sec), and Number of Time (must be Power  of 2)
     0.3    0.00 : Minimum Period (sec), and Imaginary Omega for Phinney
* MEDIUM DATA *
         2           : NL (NUMBER OF LAYERS)     t/m3,m/s,   ,m/s,   ,   m
 2.2 3000.0  100.0 0.0 1000.0  50.0 0.0  500.0 :DN,VP,Qp,Fp,Vs,Qs,Fs,Thick
 2.5 5000.0  200.0 0.0 3000.0 100.0 0.0    0.0 :DN,VP,Qp,Fp,Vs,Qs,Fs,Thick
* FAULT DATA :
    2000.0     500.0 : Total Length (STR) and Width (DIP) (m), 20
    4    1    6  : Number of Sub-Faults (NSTR,NDIP),Gauss Points per Wavelength
       0.0    1 : Time Delay(sec), Source Type(0:Ramp; 1:Smoothed; 2:exp)
       0.0       0.0     500.0 : Location of Hypocenter (X, Y and Z: m)
       0.0       0.0     500.0 : Location of Fault Origin (X, Y and Z: m)
       0.0      90.0 : Strike, Dip (deg)
    2000.0 : Rupture Velocity (m/s)
       1     1.0 : Number of Time Windows and Time Interval (sec) 20
     1.0     1.0 : Half-Rise time 1 & 2 for 1st Time Window (sec)
* dislocations (m) on sub-faults up to 6 time windows
 1.0
 1.0
 1.0
 1.0
* rake angles (degree) on sub-faults for 6 time windows
 0.0
 0.0
 0.0
 0.0
* Static Wavenumber Integ. Data for Greenfield's quadrature *
       2.0 : The first corner w*k on real axis (ex. 2.0)
        16 : Initial Num. of Intg. Points for Adaptive Newton-Cotes Quadrature
      10.0 : The second corner w*k on imag. axis (ex. 10.0)
        16 : Initial Num. of Intg. Points for Adaptive Newton-Cotes Quadrature
* Dynamic Wavenumber Integ. Data for Simpson's and Filon's quadratures *
       200 : Number of Integration Points from 0 to om/Ryleigh(min)
        50 : Number of Integration Points from om/Ryl(min) to om/c(final)
      10.0 : Factor for c(final): c(final)=Ryl(min)/Factor (Use 10, usually)
* CHANGE OF SIGNS OF IMAGINARY PARTS OF FINAL RESULTS (FOR FFT) *
         1 : CHANGE SIGN (3D1), NOT CHANGE SIGN (3D0)20
* OBSERVATION POINT DATA: System Origin *
         1 : NUMBER OF OP
    1500.0      50.0       0.0 : (x,y,z)
