// **************************************************
// 2D Mach = 1.5 flow over a 4.995 degree diamond
// **************************************************

fluid {
  species = air_perfect
}

reference conditions {
  pressure    = 101325.0
  temperature = 288.2
  mach        = 1.5
  alpha       = 0
  beta        = 0
  length      = 1.0
  area        = 1.e-2
}

unst {

  grid {
    file = voro_diamond_airfoil.cgns
  }

  convergence {
    relative tolerance = 1.e-12
  }

  implicit {
    do line search    = false
    number of steps   = 1000
    jacobian method   = numerical
    jacobian interval = 1
    adaptive {
      min cfl = 1.
      max cfl = 1000.
    }
  }

  convective terms {
    scheme        = MUSCL
    order         = 2
    limiter       = none
    flux function = ausmpw+
  }

  navierstokes {
    viscous = false
  }

  output {
    volume {
      format    = cgns
      interval  = 50
      variables = [q, density, mach, residuals]
    }

    surface {
      interval  = 50
      variables = [q, density, mach]
    }

    restart {
      interval = 500
    }

    loads {
      moment center = [0, 0, 0]; interval = 1; normalize = true; gauge = true;
      group_1 {name = diamond;  include = [2]}
    }

    line probe_1 {
      start            = [18, 20, 0]
      end              = [27, 20, 0]
      number of points = 1000
      interval         = 50
    }


  }

  # BC_1, BC_2, BC_3, BC_4 : Four edges of the diamond.
  # BC_5  : zmin (Symmtric since case is 2D)
  # BC_6  : ymin
  # BC_7  : xmin
  # BC_8  : ymax
  # BC_9  : xmax
  # BC_10 : zmax

  BC_1  {type  = wall; kind = slip}
  BC_2  {type  = wall; kind = slip}
  BC_3  {type  = wall; kind = slip}
  BC_4  {type  = wall; kind = slip}
  BC_5  {type  = symmetry}
  BC_6  {type  = farfield}
  BC_7  {type  = farfield}
  BC_8  {type  = farfield}
  BC_9  {type  = farfield}
  BC_10  {type = symmetry}


}