-----------------------------------------------------------------------------
        <---  Read data in the FIRST  dimension --->

1            Index'  1 (First Index), 2 (Second Index)
1           Method'  0 (Fourier), 1 (Chebyshev), 2 (Legendre)
2        Algorithm'  0 (MXM), 1 (Even-Odd), 2 (Transform)
0         Symmetry'  0 (No), 1 (Yes)

              <---  Grid Size and BlockSize  --->

65              N '  Total number of points to be worked on
0        BlockSize'  Size of the Data Block to be used for Algorithm=1,2

              <---  Grid Setup  --->

-1            Left'  Left  Physical Domain in x
 1           Right'  Right Physical Domain in x
 0           Map_F'  0 (Translation), 0 [0, 2 Pi), 1 (Linear Mapping)
1              Map'  0 (No Mapping), 1 (Angular/Tal-Ezer Mapping)
0           Manual'  0 (Machine computed alpha), 1 (hand coded alpha)
0.999        alpha'  mapping parameter
0            Angle'  Angle of points concentration for the Fourier Method.

              <---  Differentiation Filtering Setup  --->

1       Filtering?'  Pick an Filter (1,2,3,4) for the derivatives, else 0
0      Mode CutOff'  Cut off frequency
0            alpha'  0 if determined by the machine zero for exponential filter
16           Order'  Order of the exponential filter
1   Flux_Smoothing'  Smooth the Flux function BEFORE Differentiated?
1   Flux_Smoothing'  Smooth the Flux function AFTER  Differentiated?

              <---  Solution Filtering Setup  --->

1       Filtering?'  Pick an Filter (1,2,3,4) for the Solution, else 0
0      Mode CutOff'  Cut off frequency
0            alpha'  0 if determined by the machine zero for exponential filter
16           Order'  Order of the exponential filter
-----------------------------------------------------------------------------
        <---  Read data in the SECOND dimension --->

2            Index'  1 (First Index), 2 (Second Index)
0           Method'  0 (Fourier), 1 (Chebyshev), 2 (Legendre)
2        Algorithm'  0 (MXM), 1 (Even-Odd), 2 (Transform)
0         Symmetry'  0 (No), 1 (Yes)

              <---  Grid Size and BlockSize  --->

64              M '  Total number of points to be worked on
0        BlockSize'  Size of the Data Block to be used for Algorithm=1,2

              <---  Grid Setup  --->

-1          Bottom'  Bottom Physical Domain in y
 1             Top'  Top    Physical Domain in y
 1           Map_F'  0 (Translation), 0 [0, 2 Pi), 1 (Linear Mapping)
0              Map'  0 (No Mapping), 1 (Angular/Tal-Ezer Mapping)
0           Manual'  0 (Machine computed alpha), 1 (hand coded alpha)
0.999        alpha'  mapping parameter
0            Angle'  Angle of points concentration for the Fourier Method.

              <---  Differentiation Filtering Setup  --->

1       Filtering?'  Pick an Filter (1,2,3,4) for the derivatives, else 0
0      Mode CutOff'  Cut off frequency
0            alpha'  0 if determined by the machine zero for exponential filter
16           Order'  Order of the exponential filter
1   Flux_Smoothing'  Smooth the Flux function BEFORE Differentiated?
1   Flux_Smoothing'  Smooth the Flux function AFTER  Differentiated?

              <---  Solution Filtering Setup  --->

1       Filtering?'  Pick an Filter (1,2,3,4) for the Solution, else 0
0      Mode CutOff'  Cut off frequency
0            alpha'  0 if determined by the machine zero for exponential filter
16           Order'  Order of the exponential filter
-----------------------------------------------------------------------------
        <---  Read Timing data --->

1.50    Final_Time'  Final Time
1.00           CFL'  CFL condition
-----------------------------------------------------------------------------
        <---  Flux Function --->

0      Flux_Choice'  0-Linear Advection  u_t+a*u_x+b*u_y=0
                     1-Rotation          u_t-y*u_x+x*u_y=0
1 2.5  a, b       '  Wave Speed for the Linear Advection PDE (a,b)
0 0    C_x, C_y   '  Center of the Gaussian (C_x, C_y)
0.25  Radius_Ratio'  Radius ratio of the Gaussian relative to the domain size.
8      Order      '  Order of the Gaussian e**(-alpha*(r/Rad)**Order)