<!-- Taylor Green's problem - Quasi 2D -->
<?xml version="1.0" encoding="utf-8"?>

<NEKTAR>
    <EXPANSIONS>
        <E COMPOSITE="C[0]" NUMMODES="9" FIELDS="u,v,w,p" TYPE="MODIFIED" />
    </EXPANSIONS>
    <FILTERS>
        <FILTER TYPE="Checkpoint">
            <PARAM NAME="OutputFrequency">100</PARAM>
        </FILTER>
        <FILTER TYPE="ModalEnergy">
            <PARAM NAME="OutputFile"> TGV64MEnergy </PARAM>
            <PARAM NAME="OutputFrequency"> 10 </PARAM>
        </FILTER>
        <FILTER TYPE="Energy">
            <PARAM NAME="OutputFile"> TGV64Energy </PARAM>
            <PARAM NAME="OutputFrequency"> 10 </PARAM>
        </FILTER>
    </FILTERS>
    <CONDITIONS>
        <SOLVERINFO>
            <I PROPERTY="SolverType"            VALUE="VelocityCorrectionScheme"   />
            <I PROPERTY="EqType"                VALUE="UnsteadyNavierStokes"       />
            <I PROPERTY="AdvectionForm"         VALUE="Convective"                 />
            <I PROPERTY="Projection"            VALUE="Galerkin"                   />
            <I PROPERTY="Homogeneous"           VALUE="1D"                         />
            <!-- <I PROPERTY="UseFFT"                VALUE="FFTW"                       /> -->
            <I PROPERTY="GlobalSysSoln"         VALUE="DirectMultiLevelStaticCond" />
            <I PROPERTY="SpectralHPDealiasing" 	VALUE="True" />
            <I PROPERTY="SpectralVanishingViscosity"         VALUE="True" />
        </SOLVERINFO>
        <TIMEINTEGRATIONSCHEME> 
            <METHOD> IMEX </METHOD> 
            <ORDER> 2 </ORDER> 
        </TIMEINTEGRATIONSCHEME>
        <PARAMETERS>
            <P> TimeStep       = 0.01    </P>
            <P> FinalTime      = 20.0     </P>
            <P> NumSteps       = FinalTime/TimeStep </P>
            <P> IO_InfoSteps   = 10       </P>
            <P> IO_CFLSteps    = 50       </P>
            <P> Re             = 1600     </P>
            <P> Kinvis         = 1/Re     </P>
            <P> V0             = 1        </P>
            <P> L              = 1        </P>
            <P> LZ             = 2*PI     </P>
            <P> HomModesZ      = 64      </P>
            <P> SVVCutoffRatio      = 0.7      </P>
            <P> SVVDiffCoeff      = 0.1      </P>
        </PARAMETERS>
        <VARIABLES>
            <V ID="0"> u </V>
            <V ID="1"> v </V>
            <V ID="2"> w </V>
            <V ID="3"> p </V>
        </VARIABLES>
        <BOUNDARYREGIONS>
            <B ID="0"> C[1] </B>
            <B ID="1"> C[2] </B>
            <B ID="2"> C[3] </B>
            <B ID="3"> C[4] </B>
        </BOUNDARYREGIONS>
        <BOUNDARYCONDITIONS>
            <REGION REF="0">
                <P VAR="u" VALUE="[2]" />
                <P VAR="v" VALUE="[2]" />
                <P VAR="w" VALUE="[2]" />
                <P VAR="p" VALUE="[2]" />
            </REGION>
            <REGION REF="1">
                <P VAR="u" VALUE="[3]" />
                <P VAR="v" VALUE="[3]" />
                <P VAR="w" VALUE="[3]" />
                <P VAR="p" VALUE="[3]" />
            </REGION>
            <REGION REF="2">
                <P VAR="u" VALUE="[0]" />
                <P VAR="v" VALUE="[0]" />
                <P VAR="w" VALUE="[0]" />
                <P VAR="p" VALUE="[0]" />
            </REGION>
            <REGION REF="3">
                <P VAR="u" VALUE="[1]" />
                <P VAR="v" VALUE="[1]" />
                <P VAR="w" VALUE="[1]" />
                <P VAR="p" VALUE="[1]" />
            </REGION>
        </BOUNDARYCONDITIONS>
        <FUNCTION NAME="InitialConditions">
            <E VAR="u" VALUE="V0 * sin(x/L)*cos(y/L)*cos(z/L)" />
            <E VAR="v" VALUE="-1 * V0 * cos(x/L)*sin(y/L)*cos(z/L)" />
            <E VAR="w" VALUE="0" />
            <E VAR="p" VALUE="0" />
        </FUNCTION>
    </CONDITIONS>
</NEKTAR>
