* OMEGA DATA *
         1        20       0.2 : IOM(INITIAL),NOM(FINAL),DOM(Delta Omega)
* MEDIUM DATA *
         3           : NL (NUMBER OF LAYERS)      t/m3,m/s,  ,m/s,  , m
    2.00  2000.0    50.0  1000.0    50.0  1000.0 : DNS, VP,QP, VS,QS,THICK 1
    2.20  3000.0   100.0  2000.0   100.0  1000.0 : DNS, VP,QP, VS,QS,THICK 2
    2.50  5000.0   200.0  3000.0   200.0     0.0 : DNS, VP,QP, VS,QS,THICK 3
* WAVENUMBER INTEGRATION DATA FOR SIMPSON'S RULE *
         1 : PATTERN 1 (C1,N1,C2,N2) OR 2 (WB,N1,WC,N2)
     200.0       200     800.0       100 : C1,N1,C2,N2 (C1<C2) for Pattern 1
    0.0006       200     0.010       200 : WB,N1,WC,N2 (WB<WC) for Pattern 2
* SOURCE POINT DATA *
         2     900.0 : SP LAYER, SP DEPTH FROM UPPER BOUNDARY OF THE LAYER (m)
       0.0       1.0       1.0 : QX, QY, QZ
         1 : RICKER WAVELET (=1) OR ORIGINAL WAVE (=2)
       3.0       3.0 : DOMINANT PERIOD AND DELAY TIME FOR RICKER WAVE (sec)
       0.0 : N OF AA**N WHERE AA=DN1*VS1*VS1 (AMPLIFICATION FACTOR FOR RICKER)
* OBSERVATION POINT DATA *
         2     800.0 : OP LAYER, OP DEPTH FROM UPPER BOUNDARY OF THE LAYER (m)
         2       0.0 : NUMBER OF OP, and Angle from the x-axis (Unit: degree)
    5000.0   10000.0       0.0 : R1,R2, AND R3 (unit: m)
* DATA FOR THE ASYMPTOTIC SOLUTIONS *
         1 : BE USED (=1), OR NOT USED (=0)
* CHANGE OF SIGNS OF IMAGINARY PARTS OF FINAL RESULTS (FOR FFT) *
         1 : CHANGE SIGN (=1), NOT CHANGE SIGN (=0) 


