time
{
    steady = false
    #number of steps = 10000
    number of steps = 1500
    cfl = 2.5
}

solver defaults
{
    explicit {method = ssprk53}
}

reference conditions
{
    pressure    = 101327.
    temperature = 300.0
    mach = 1.5
    beta = -45
    length = 1
}

fluid
{
    species = air_perfect
}

cartesian
{


    navier stokes
    {
        viscous = false
    }

    convective terms
    {
       scheme = hwcns                                     # variant of the Hybrid Weighted Compact Non-linear Scheme, https://arc.aiaa.org/doi/pdf/10.2514/6.2011-3857
       flux function = hllec                              # Robust Harten-Lax-Leer-Einsfeldt-Contact approximate Riemann solver (https://doi.org/10.1007/BF01414629)
       fluxorder = 4                                      # formally 3rd or 4th order convective terms via nodal flux blending with face fluxes
       interporder = 3                                    # third order left and right state interpolation to the face
       limiter = weno3op                                  # 3rd order Weighted Essentially Non-Oscillatory interpolation with optimal accuracy, https://epubs.siam.org/doi/pdf/10.1137/19M1260396
       upwind blending sensors = [weno]                   # use WENO-type sensor on density to detect stencils with discontinuities and use upwind there and central otherwise
       central components = [q]                           # blend all primitive variables between upwind and central interpolation based on upwind blending sensors
    }

    output
    {
        timeaveraging {  start ctu = 0.695 }

        monitor{
            interval  = 100
            variables = [nt,time, dt,maxcfl,blocks_1*,numcells,mups,memory,elapsed]
        }

        volume
        {
             interval  = 1000
             variables = [q,mach,density,mask,wenoshocksensor,errormetric]
        }

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

        restart
        {
             interval  = 10000
        }

        line probe_1 {start = [-2.0, 4.0];end = [2.0, 4.0];number of points = 100;interval  = 100}

    }

    bcs
    {
        ibmethod = gcm
        gcmtype  = imagedist
        bc_1
        {
            type = wall
            kind = [slip]
        }

    }


    grid
    {
        finest dx = 0.01
        targetrelativefarfielddistance = [200^3]
        regrid interval = 100
        max box size = 8
        tag all geometry = true

        solution based amr
        {
          method = auto
          max cells = 600000
          error metric time method = ins
          error metric stats interval = 100
          error metric threshold ratio = 0.5
          regridgrowerrormetrictags = 1

          errormetric_1
          {
            type = laplacianp
            weight = 1.0
          }       

        }
    
        #fine box over entire airfoil and nearby
        tagboxes_1
        {
          lo = [-0.55, -0.55]
          hi = [ 0.55, 0.55]
          dx = 0.0
        }


    }

    geometry
    {
        surfaces_1
        {
            file name =diamond.fmt; rotatedeg = [45,0]
        }
    }
}