#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------MESH PARAMETERS-----#
# Model diameter ~ 1.477m
#reference_length: 0.02954   # x-coarse (50 points per diameter)
#reference_length: 0.01477   # coarse (100 points per diameter)
reference_length: 0.007385  # medium (200 points per diameter)
#reference_length: 0.0036925 # fine (400 points per diameter)

# 2d or 3d mesh
2d: false
#~~~~~~~~~~~~~~~~~~~~~~~~~~#

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------POINT CLOUD PARAMETERS-----#
# if reading point cloud from a csv file specify true
read_points_from_file: false

# point cloud csv file name
# You may specify *.csv file type or as *.bin file type which are produced by the point cloud export functionality.
point_cloud_file: smoothed_final.csv

# when saving the unsmoothed point cloud and the final smoothed point cloud, you may specify the file format
# which is saved to the file system. If this variable is not set as "csv" or "bin" then we default to csv.
point_cloud_export_type: csv
#point_cloud_export_type: bin

# if marking wall points and smoothing points specify true
# User may want to specify false if a smoothed point cloud which was previously labeled is used.
do_pointCloud_labeling: true

# total number of lloyd iterations
lloyd_iter: 10
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------PARTITION PARAMETERS-----#
#number of partitions used while generating the mesh
num_partitions: 1

#specify true if user wants to restart from a previously defined partition.
read_partition_from_file: false

#partition point cloud file name
partitionPCfile: partitionPC.csv

#convergence criterion for load balancing
load_balancing_rms_threshold: 0.05

#do lloyd smoothing without recomputing true local graph for N iterations, then recompute it
recompute_local_graph_every: 1

#blending between CVT and VP algorithms in load balacing
load_balancing_blending: 0.8

#max number of iterations in load balancing if convergence is not met
load_balancing_max_iter: 2000

#distance multiplier when trying to find ghost cells.
#value needs to be >= 1.
#An increased value means more ghosts but also means more safety with respect to the connectivity of ghost split cells.
dist_coeff: 1.5

#specify true if user wants to check if all faces have neighbors.
face_connectivity_check: false
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------OUTPUT PARAMETERS-----#
# Specify if scalars have to be written into the
# CGNS mesh
write_cgns_scalars: false
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~#

regions:
  # First, specify the container.
  # Container layers are zero by default
  # and angles do not work for containers
  - name: farfield
    type: cube
    distance: [32, 24, 24]
    seed: bodycentered
    spacing_rel: [24, 24, 24]
    center: [10, 0, 0]

  #--------------------------------------------
  # Near body refinements
  #--------------------------------------------
  # PID: 1   name: front
  # PID: 2   name: shoulder
  # PID: 3   name: middle
  # PID: 4   name: back_fillet
  # PID: 5   name: back

  - name: body
    type: wall
    smoothing_layers: 4
    project_layers: 2
    tangential_plane_pert_factor: 0.1
    seed_toward_normal: true
    trifile: genesis_capsule_surface_tri.cgns
    face_merge_angle: 90

    refinements:
      - name: telescope
        seed: bodycentered
        smoothing_layers: 2 # how many smoothing layers per side
        num_levels: 4
        # Absolute level thickness specification to make sure we capture the bow shock
        level_spacings_rel: [2  , 4  , 8  , 16 ]
        level_thicknesses:  [0.4, 0.6, 1.0, 1.8]
        components: 1-5
      - name: wall_inner
        seed: bodycentered
        smoothing_layers: 2 # how many smoothing layers per side
        level_points_per_thickness: 4
        level_spacing_inner_rel: 1
        level_spacing_outer_rel: 1
        level_spacing_growth_rate: 2
        components: 1-3

  - name: wake
    type: frustum
    axis_start: [0.36, 0, 0]
    axis_end: [6.3, 0, 0]
    radius_start: 1.14
    radius_end: 1.14
    seed: bodycentered
    smoothing_layers: 2
    spacing_rel: [2, 2, 2]
    telescoping_refinement:
      - spacing_growth: 1.5
        num_points_per_layer: 8
        max_spacing_rel: 24