/* -------------------------------------------------------------
    This is a sample DG problem definition file                           
    Made by Jean-Francois Remacle in 2002.
    (C) 2002 Scientific Computation Research Center

    This is the classical backward facing step problem,
    a mach 3 flow in which a step is impulsively inserted.

    Comments can be defined the same way as they are in a C++
    file. Note that comments into comments will make the parser
    to crash !!!
   ------------------------------------------------------------- */

/*
    Reserved word : Gamma is the gas constant, default = 1.4
*/

Gamma = 1.4;

/* 
    This file defines a external flow problem
    The mesh has a certain number of parametres
    
    Mach  : Mach number
    Angle : Angle that the flow has with the X axis     
*/

Mach =  3;
Angle = 0;

/*
    We defined some variables to help us
*/

rho   = 1.0;
p     = 1.0;
v     = Sqrt [Gamma] * Mach;
v_x   = Cos  [Angle] * v   ;
v_y   = Sin  [Angle] * v   ;

Function
{
  /// Initial Conditions : those names required !!!
  /// It means DG user should always use thos names
  /// to define the initial pressures, velocity and
  /// densities !!!
  PressureInit []  = Const [p];
  VelocityInit []  = Vector [ v_x , v_y , 0.0 ];
  DensityInit  []  = Const [rho];

  /// Same for, Boundary Conditions : those names required !!!
  /// Boundary Conditions, as well as initial conditions, 
  /// may be defined as piecewise functions, the pieces
  /// being model regions or faces (we call everything region
  /// here, but the word Face or Edge is allowed).
  PressureBnd [Region[{1,8,9}]]  = Const [p];
  VelocityBnd [Region[{1,8,9}]]  = Vector [ v_x , v_y , 0.0 ];
  DensityBnd  [Region[{1,8,9}]]  = Const [rho];
  /// for post processing , the user can define any function
  /// he needs. Some variables preceeded by $ are current
  /// quantities : 
  /// $X, $Y and $Z for coordinates, $Time for current time,
  /// $Pressure, $Density , $VX, $VY, $VZ as velocity vector
  NrmVelocity [] = Sqrt [ $VX * $VX + $VY * $VY + $VZ * $VZ ];
  Pressure    [] = $Pressure;
  Density     [] = $Density;
}


/// Simulation parameters

/// The mesh
Mesh				= "step.msh";
/// Where you want to put the results
ChunkSpace                      = ".";
/// We use a limiter (supersonic flow with shocks)
Limiter                         = 1;
/// No shock detector, limit everywhere 
/// SENSITIVE PARAMETER USE WITH CARE !!!
Detector                        = 0.0;
/// Time parameters
/// Security coefficient for the time step
CFL_Limit         	       	= 0.15;
/// Initial time
Time                        	= 0.0;
/// Final Time
Tend                        	= 4.0;
/// Use local time stepping ( 1 : no local time stepping)
LocalTimeStepping               = 2;
/// Goal time step for the local time stepping
GoalTimeStep                	= 2.e-3;
/// conforming refinement parameters
H_Min_CellSize              	= 1.e-2;
H_Max_CellSize              	= 1.e-1;
H_Max_AspectRatio           	= 1.0;
H_SamplePeriod              	= 2.e-2;
/// we use conforming H refinement ...
H_Refinement                    = 1;   
H_Conforming                    = 0;
H_Anisotropic                   = 0;
/// 2 levels of refinement
/// THIS IS A SENSITIVE PARAMETER
/// IF = 5 e.g., simulation may take a while ;-)
H_MaxLevelOfRefinement          = 2;

/// Sensors can be added to the simulation
/// for atomatically generating post-processing
/// files. The quantities that are exported are
/// defined into the Function section.
/// Several kind of sensors are available
/// Gmsh for gmsh output format (colormaps)
/// Tecplot for tecplot outputs (colormaps)
/// Ensight for ensight outputs (colormaps, to come ...)
/// Gnuplot for gnuplot outputs (curves)

Sensor 
{
  Type       = "Gmsh";
  SampleRate = 1.e-1;
  Quantity   = "Pressure";
  Name       = "BFS_Pressure";
}
Sensor 
{
  Type       = "Gmsh";
  SampleRate = 1.e-1;
  Quantity   = "Density";
  Name       = "BFS_Density";
}