<?xml version="1.0" encoding="utf-8"?>

<NEKTAR xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
    xsi:noNamespaceSchemaLocation="http://www.nektar.info/nektar.xsd">

    <EXPANSIONS>
        <E COMPOSITE="C[10-11]" NUMMODES="3" FIELDS="rho,rhou,rhov,E" TYPE="MODIFIED" />
    </EXPANSIONS>

    <CONDITIONS>
        <SOLVERINFO>
	      <I PROPERTY="EQType"                VALUE="NavierStokesCFE"     />
	      <I PROPERTY="Projection"            VALUE="DisContinuous"       />
	      <I PROPERTY="AdvectionType"         VALUE="WeakDG"              />
	      <I PROPERTY="DiffusionType"         VALUE="LDGNS"               />
	      <I PROPERTY="TimeIntegrationMethod" VALUE="ForwardEuler"        /> <!-- ClassicalRungeKutta4 -->
	      <I PROPERTY="UpwindType"            VALUE="HLLC"                />
	      <I PROPERTY="ProblemType"           VALUE="General"             />
	      <I PROPERTY="EquationOfState"       VALUE="IdealGas"            />
	      <!--<I PROPERTY="ViscosityType"         VALUE="Variable"            /> -->
        </SOLVERINFO>

	    <PARAMETERS>
	      <!--Computational parameters-->
	      <P> TimeStep            = 0.33  		                </P>  
	      <P> NumSteps            = 7500       			</P>
	      <P> IO_CheckSteps       = 2500         			</P> 
	      <P> IO_InfoSteps        = 200            		        </P>
	      <P> IO_CFLSteps         = 200            		        </P>
	      <!-- Fluid properties -->
	      <P> Gamma               = 1.4          			</P>
	      <P> Pr                  = 0.72         			</P>
	      <P> Mach                = 0.5         			</P> <!-- 0.6 -->
	      <P> Re                  = 2000.0      			</P>
	      <P> GasConstant         = 1.0/(Gamma*Mach*Mach)		</P>
	      <P> mu                  = 1.0/Re     			</P>
	      <!-- Boundary condition parameters-->
	      <P> rhoInf              = 1.0            			</P>
	      <P> TInf                = 1.0            			</P>
	      <P> pInf                = rhoInf*GasConstant*TInf	        </P>
	      <P> cInf                = sqrt(Gamma*GasConstant*TInf)    </P>
	      <P> uInf	              = Mach*cInf	       		</P>
	      <P> vInf	              = 0.0	       		        </P>
	    </PARAMETERS>
	    
	    <VARIABLES>
	      <V ID="0"> rho  </V>
	      <V ID="1"> rhou </V>
	      <V ID="2"> rhov </V>
	      <V ID="3"> E    </V>     
	    </VARIABLES>


        <BOUNDARYREGIONS>
            <B ID="0"> C[9] </B>  <!-- wall !-->
            <B ID="1"> C[7] </B>    <!-- outlet !-->
            <B ID="2"> C[8] </B>    <!-- farfield -->
            <B ID="3"> C[6] </B>    <!-- inlet !-->
        </BOUNDARYREGIONS>

        <BOUNDARYCONDITIONS>
      		<!-- Wall -->    
	      <REGION REF="0">
		<D VAR="rho"  USERDEFINEDTYPE="WallAdiabatic" VALUE="0" />
		<D VAR="rhou" USERDEFINEDTYPE="WallAdiabatic" VALUE="0" />
		<D VAR="rhov" USERDEFINEDTYPE="WallAdiabatic" VALUE="0" />
		<D VAR="E"    USERDEFINEDTYPE="WallAdiabatic" VALUE="0" />
	      </REGION>
	      <!-- Inflow -->
	      <REGION REF="3">  
		<D VAR="rho"  VALUE="rhoInf"      />
		<D VAR="rhou" VALUE="rhoInf*uInf" />
		<D VAR="rhov" VALUE="rhoInf*vInf" />
		<D VAR="E"    VALUE="pInf/(Gamma-1) + 0.5*rhoInf*(uInf*uInf + vInf*vInf)" />
	      </REGION>
	      <!-- Outflow -->
	      <REGION REF="1">
		<D VAR="rho"  VALUE="rhoInf"      />
		<D VAR="rhou" VALUE="rhoInf*uInf" />
		<D VAR="rhov" VALUE="rhoInf*vInf" />
		<D VAR="E"    VALUE="pInf/(Gamma-1) + 0.5*rhoInf*(uInf*uInf + vInf*vInf)" />
	      </REGION>
	      <!-- Farfield -->
	      <REGION REF="2">  
		<D VAR="rho"  VALUE="rhoInf"      />
		<D VAR="rhou" VALUE="rhoInf*uInf" />
		<D VAR="rhov" VALUE="rhoInf*vInf" />
		<D VAR="E"    VALUE="pInf/(Gamma-1) + 0.5*rhoInf*(uInf*uInf + vInf*vInf)" />               
	      </REGION>
        </BOUNDARYCONDITIONS>

        <FUNCTION NAME="InitialConditions">
	      <E VAR="rho"  VALUE="rhoInf" />	         
	      <E VAR="rhou" VALUE="rhoInf*uInf" />	        
	      <E VAR="rhov" VALUE="rhoInf*vInf" />	        
	      <E VAR="E"    VALUE="pInf/(Gamma-1) + 0.5*rhoInf*(uInf*uInf + vInf*vInf)" />
        </FUNCTION>


    </CONDITIONS>

</NEKTAR>
