miguelplate.afl
$GLOBAL,
	TITLE    = 'Shock-BL Interaction',
	SIMTYPE  = 2,
	FOLDER   = 'shockbl',
	MACH     = 2.0,
	REYNOLDS = 300000.0,
	ORDER    = 2,
	BCORDER  = 1,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE TURBULENCE MODEL
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

COMMENT:~~~~~~~~~~~~SPALART-ALLMARAS~~~~~~~~~~~~~~~
$TURB,
    TURBTYPE = 0,
    $END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE SPATIAL SCHEME
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



$SPATIAL,
    SCHEME   = 5, 5, -999,
    	COMPACTI = 'C4-AC4-C6-AC4-C4',
    	COMPACTJ = 'C4-AC4-C6-AC4-C4',
    	COMPACTK = 'C4-AC4-C6-AC4-C4',
    	COMPACTG = 'C4-AC4-C6-AC4-C4',	
    METRIC   = 0,
    VISCOUS  = 1,  1,  0,
    VISCHEME = 1,
    VISCMON  = 1,
    CUTOFFI =   5.0000001E-02 ,
    ISOTROPI =            0 ,
    ICUT =            0 ,           0 ,           0 ,
    CUTOFFJ =   5.0000001E-02 ,
    ISOTROPJ =            0 ,
    JCUT =            0 ,           0 ,           0 ,
    SUTHERLN = 0.38,
    $END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE TIME INTEGRATION SCHEME
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
$TIMESTEP,
	MAXITER  = 80000,
	SCHEME   = 'BW2',
	KTVDRK   = 0,
	GLOBALDT = 1.E-3,
	SUBDT    = 2.E-4,
	SUBON    = 1,
	MAXSUB   = 2,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE DAMPING/FILTERING
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

 $DAMPING,
        TYPE =            2 ,
        ES2 =   0.2000000     ,
        ES4 =   0.02 ,
        FES2 =    10.000000     ,
        FES4 =    2.000000     ,
        MAXRED =            1 ,
        OMGAV =   0.5 ,
        SRCONST =   0.0000000E+00 ,
        FILTER =            1 ,           0 ,           0 ,
        NFILTER =            1 ,           0 ,           0 ,
        ALPHAN =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA1 =   0.0E+00     ,  0.0E+00     ,  0.0E+00 ,
        ALPHA2 =   0.49000     ,  0.49000     ,  0.49000     ,
        ALPHA3 =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA4 =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA5 =   0.45000     ,  0.45000     ,  0.45000     ,
        ORDERN =            8 ,           8 ,           8 ,
        ORDER1 =            0 ,           0 ,           0 ,
        ORDER2 =            2 ,           2 ,           2 ,
        ORDER3 =            6 ,           6 ,           6 ,
        ORDER4 =            6 ,           6 ,           6 ,
        ORDER5 =            8 ,           8 ,           8 ,
      $END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    OUTPUT AND DEBUG SETUP
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

$OUTPUT,
	FORMAT = 'TECP',
	PRINTFRQ = 100,
    	TECANIM  = 0,
    	TRANGE   = 0,1000000,10,
	$END

$DEBUG,
	POINT = 4, 4, 2,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE GRID
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
$BLOCK,
	NAME     = 'block01',
	MESHFILE = 'myres.cgns',
	FILETYPE = 'CGNS',
	ORDER    = 3,
	ISIZE    = 286,
	JSIZE    = 131,
	KSIZE    = 1,
	PERIODIC = 0, 0, 3,
	$END
    MESHFILE = 'flatplate_1.tmp',
    FILETYPE = 'AEROFLO',

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:      DEFINE BOUNDARY CONDITIONS
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

COMMENT:~~~~~~~~~~~~I-DIRECTION~~~~~~~~~~~~~~~~~~~~

comment:	Outflow condition at i= IE
$FACEBC,
	BLOCK = 'block01',
	TYPE = 31,
	SIDE = 2,
	$END

COMMENT:~~~~~~~~~~~~J-DIRECTION~~~~~~~~~~~~~~~~~~~~

comment:	Free stream condition at J= JE
$FACEBC,
	BLOCK = 'block01',
	TYPE = 32,
	SIDE = 6,
	$END

comment:    Equation Rho = 1.4*Ma*Ma*P/1.8 at J = 1, 
$FACEBC,
    BLOCK    = 'block01',
    TYPE     = 80,
    SIDE     = 5,
    VARIABLE = 1,
    VAR1     = 5,
    OPERATO1 = '*',
    VAL1     = 3.11111111,
    $END

comment:    Wall at J = 1
$FACEBC,
    BLOCK    = 'block01',
    TYPE     = 1,
    SIDE     = 5,
    $END