fluid {
  species = air_perfect
}

reference conditions {
  mach        = 0.2
  Re          = 5.9e6
  temperature = 288.15      # K
  alpha       = 19.7
  beta        = 0.00
  length      = 7.00532     # m
  area        = 191.844778  # Half Sref m2
  dx          = 4e-3
}

time {
  steady      = false
  finish ctu  = 90
  initial cfl = 0.05
  cfl         = 2.
}

unst {

  explicit {
    method = ssprk33
  }

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

  grid {
    file = voro_hlpw5_case2p4.cgns
  }

  diagnostics {
    target size filter   = 4e-3
    max percent filtered = 0.05
  }

  convective terms {
    scheme                  = KEEPBLEND
    order                   = 2
    limiter                 = minmod
    flux function           = ausmpw+
    upwind blending sensors = [ducrosmod]
  }

  navierstokes {
    viscous = true
  }

  turbulence {
    sgs model = vreman
    sgs const = 0.08
  }

  wall function {
    type                 = blendedlog
    do transition        = true
    transition threshold = 0.05
  }

  trip time relaxation = 40. # a factor of dt
  # trip_1 : Wing Upper Near Root
  # trip_2 : Wing Upper Near Nacelle Pylon
  # trip_3 : Wing Upper Near Tip
  # trip_4 : Nacelle outer circle
  # trip_5 : Nacelle inner circle
  # trip_6 : Nacelle pylon
  # trip_7 : Htail LE upper
  # trip_8 : Htail LE lower
  # trip_9 : Vtail LE
  # trip_10: Fuse nose zplus
  # trip_11: Fuse nose zminus
  # trip_12: Fuse nose yplus
  trip_line_1  {dxscale = 1.;begin = [2.546070E+01, 3.052050E+00, 4.729215E+00];end = [2.618612E+01, 3.971736E+00, 4.802266E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_2  {dxscale = 1.;begin = [3.000402E+01, 8.994870E+00, 5.215469E+00];end = [3.083125E+01, 1.007270E+01, 5.302925E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_3  {dxscale = 1.;begin = [4.507610E+01, 2.900512E+01, 6.720341E+00];end = [4.536383E+01, 2.937944E+01, 6.747311E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_4  {dxscale = 1.;center = [24.861590843432634, 9.803591119112955, 2.3914368745527694];axis = [-0.998535326, -0.030048636, -0.044992013];radius = 1.8;projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_5  {dxscale = 1.;center = [25.215387395427886, 9.781733682555998, 2.4194252651721597];axis = [-0.993729501, 0.022842868, -0.109452649];radius = 1.4;projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_6  {dxscale = 2.;begin = [2.634552E+01, 9.469491E+00, 4.450000E+00];end = [2.633544E+01, 9.758339E+00, 4.450000E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_7  {dxscale = 2.;begin = [5.739016E+01, 1.680946E+00, 6.745135E+00];end = [6.467291E+01, 1.066800E+01, 7.780111E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_8  {dxscale = 2.;begin = [5.734611E+01, 1.690681E+00, 6.252033E+00];end = [6.465877E+01, 1.066800E+01, 7.645515E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_9  {dxscale = 2.;begin = [5.402829E+01, 4.588661E-01, 8.514268E+00];end = [6.378401E+01, 1.260798E-01, 1.842068E+01];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_10 {dxscale = 2.;begin = [2.519835E+00, 0.000000E+00, 5.436731E+00];end = [2.519835E+00, 4.400000E-01, 5.436731E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_11 {dxscale = 2.;begin = [2.519835E+00, 0.000000E+00, 4.627336E+00];end = [2.519835E+00, 4.400000E-01, 4.627336E+00];projectionmethod = closest;projectionmaxdistance = 3;}
  trip_line_12 {dxscale = 2.;begin = [2.519835E+00, 4.400000E-01, 4.627336E+00];end = [2.519835E+00, 4.400000E-01, 5.436731E+00];projectionmethod = closest;projectionmaxdistance = 3;}

  output {

    time averaging {
      start ctu = 40
    }

    monitor {
      interval  = 10
      variables = [nt, time, ctu, dt, maxcfl, maxmach, cl_airplane, cl_wing_fuse, wallt, mups, elapsed]
    }

    #volume {
    #  write initial = false
    #  format        = cgns
    #  interval     = 10000
    #  variables     = [q,qstats,density,mach,upwind]
    #}

    surface {
      interval  = 10000
      variables = [q, mach, tauw, pcoef, cf, qavg, yplusavg, tauwavg, turbulence];
    }

    restart {
      interval = 10000
    }

    loads {
      moment center = [33.67786, 0.0, 4.51993]
      interval      = 10
      normalize     = true
      gauge         = true
      group_1 {name = airplane;  include = [1-33]}
      group_2 {name = wing;      include = [1, 3-9, 12, 13, 15, 16, 19, 20]}
      group_3 {name = tail;      include = [25-33]}
    }

    isosurfaces_1  { contour = cutplane; point = [0, 5.08, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #A
    isosurfaces_2  { contour = cutplane; point = [0, 7.76, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #B
    isosurfaces_3  { contour = cutplane; point = [0, 10.18, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #C
    isosurfaces_4  { contour = cutplane; point = [0, 12.60, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #D
    isosurfaces_5  { contour = cutplane; point = [0, 15.15, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #E
    isosurfaces_6  { contour = cutplane; point = [0, 17.96, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #F
    isosurfaces_7  { contour = cutplane; point = [0, 20.76, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #G
    isosurfaces_8  { contour = cutplane; point = [0, 23.18, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #H
    isosurfaces_9  { contour = cutplane; point = [0, 26.37, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #I
    isosurfaces_10 { contour = cutplane; point = [0, 28.92, 0]; normal = [0, 1, 0]; format = cgns; interval = 5000; colorfield = q } #J
    #isosurfaces_11 { contour = qcriterion; isovalue = 1e6; format = cgns; interval = 5000; colorfield = umag} # iso-surface of q-criterion colored by velocity magnitude

  }

  # BC_1 : flaps_upper
  # BC_2 : nacelle_inner
  # BC_3 : wing_le_upper
  # BC_4 : wing_te_upper
  # BC_5 : wing_te_lower
  # BC_6 : wing_le_lower
  # BC_7 : wing_flap_cove
  # BC_8 : flap_track_fairings_connection
  # BC_9 : flap_track_fairings_back
  # BC_10: nacelle_te
  # BC_11 : nacelle_le
  # BC_12 : wing_lower
  # BC_13 : wing_upper
  # BC_14 : nacelle_outer
  # BC_15 : slats
  # BC_16 : flaps
  # BC_17 : vortex_generator
  # BC_18 : nacelle_pylon
  # BC_19 : slat_attachments
  # BC_20 : flap_track_fairings_front
  # BC_21 : fuse_fairing
  # BC_22 : fuse_aft
  # BC_23 : fuse_forward
  # BC_24 : fuse_nose
  # BC_25 : htail_te
  # BC_26 : htail_le
  # BC_27 : htail_upper
  # BC_28 : htail_lower
  # BC_29 : htail_tip
  # BC_30 : vtail_te
  # BC_31 : vtail_le
  # BC_32 : vtail
  # BC_33 : vtail_tip
  # BC_34 : zMin
  # BC_35 : yMin
  # BC_36 : xMin
  # BC_37 : yMax
  # BC_38 : xMax
  # BC_39 : zMax

  BC_1  {type = wallfunction}
  BC_2  {type = wallfunction}
  BC_3  {type = wallfunction}
  BC_4  {type = wallfunction}
  BC_5  {type = wallfunction}
  BC_6  {type = wallfunction}
  BC_7  {type = wallfunction}
  BC_8  {type = wallfunction}
  BC_9  {type = wallfunction}
  BC_10 {type = wallfunction}
  BC_11 {type = wallfunction}
  BC_12 {type = wallfunction}
  BC_13 {type = wallfunction}
  BC_14 {type = wallfunction}
  BC_15 {type = wallfunction}
  BC_16 {type = wallfunction}
  BC_17 {type = wallfunction}
  BC_18 {type = wallfunction}
  BC_19 {type = wallfunction}
  BC_20 {type = wallfunction}
  BC_21 {type = wallfunction}
  BC_22 {type = wallfunction}
  BC_23 {type = wallfunction}
  BC_24 {type = wallfunction}
  BC_25 {type = wallfunction}
  BC_26 {type = wallfunction}
  BC_27 {type = wallfunction}
  BC_28 {type = wallfunction}
  BC_29 {type = wallfunction}
  BC_30 {type = wallfunction}
  BC_31 {type = wallfunction}
  BC_32 {type = wallfunction}
  BC_33 {type = wallfunction}
  BC_34 {type = farfield}
  BC_35 {type = symmetry}
  BC_36 {type = farfield}
  BC_37 {type = farfield}
  BC_38 {type = farfield}
  BC_39 {type = farfield}

}