time {
  steady          = false
  number of steps = 30000
  dt              = 9.12e-5
}

reference conditions {
  Re          = 150
  temperature = 277.7777778
  mach        = 0.2
  alpha       = 0
  beta        = 0
  length      = 1.
  area        = 1.044e-3
}

fluid {
  species = air_perfect
}

unst {

  #restart {
  #  file = restart/Unst_restartLatest.hdf5
  #}

  convergence {
    relative tolerance = 1.e-3;
    absolute tolerance = 1.e-12;
  }

  grid {
    file = voro_cylinder_2d.cgns
  }

  implicit {
    jacobian method = numerical
    method          = bdf2
    number of steps = 16
    adaptive {
      min CFL = 5.
      max CFL = 100.
    }
  }

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

  navier stokes {
    viscous                  = true
    relative tolerance       = 1e-3
    q change threshold       = 0.5
    density change threshold = 0.5
  }

  output {

    loads {
      moment center = [0., 0., 0.]
      interval      = 1
      normalize     = true
      gauge         = true
      group_1 {
        name    = cylinder
        include = [1]
      }
    }

    volume {
      interval  = 500
      variables = [q, vorticity, mach, density]
    }

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

    restart {
      interval = 1000
    }
  }

  #Cylinder
  BC_1 {type = wall; kind = noslip}
  #Z min
  BC_2 {type = symmetry}
  #Y min
  BC_3 {type = farfield}
  #X min
  BC_4 {type = farfield}
  #Y max
  BC_5  {type = farfield}
  #X max
  BC_6  {type = farfield}
  #Z max
  BC_7  {type = symmetry}
}