#~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------MESH PARAMETERS-----#
#length scale
reference_length: 1.0

# 2d or 3d mesh
2d: true

2d_plane_output:
  - active: false
    triangulate: false
    normalVector: [0, 1, 0]
    pointInPlane: [0,0.5,0]
#~~~~~~~~~~~~~~~~~~~~~~~~~~#

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
#------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: anisotropic_seeds_final.csv
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

near_wall_cell_type: isotropic
# near_wall_cell_type: anisotropic
#near_wall_cell_type: elongated

num_wall_pnt_lloyd_iters: 10

max_ray_angle: 10

merge_restricted_triangles_in_final_mesh: true

# 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: true

#partition point cloud file name
partitionPCfile: partitionPC.csv

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

#convergence criterion for load balancing
load_balancing_rms_threshold: 0.05

#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: container
    type: cube
    distance: [200, 200, 1.044e-3]
    seed: hexprism
    spacing_rel: [1e1, 1e1, 1e1]
    center: [0.0, 0.0, 0.0]

  - name: circle
    type: wall
    seed_toward_normal: true
    # Create a cylinder object
    trifile: cylinder
    geom_parameters:
    - center: [0, 0, 0]
      axis: [0, 0, 1]
      radius: 0.5
      height: 2.5
      angle_resolution: 2
    isDomainBoundary: true
    smoothing_layers: 3
    project_layers: 2
    face_merge_angle: 0
    refinements:
        - name: wall_telescope
          seed: hexprism
          smoothing_layers: 2
          level_points_per_thickness: 8
          level_spacing_inner_rel: 2e-3
          level_spacing_outer_rel: 10.0
          level_spacing_growth_rate: 1.5
          components: 1
        - name: wall_inner
          seed: hexprism
          smoothing_layers: 2
          level_points_per_thickness: 4
          num_levels: 1
          level_spacing_inner_rel: 1.044e-3
          level_spacing_outer_rel: 1.044e-3
          level_spacing_growth_rate: 1.
          components: 1

  - name: wake-near
    type: telescoping
    seed: hexprism
    smoothing_layers: 8
    points_per_bbox: [8, 8, 8]
    frustum:
      - axis_start: [0, 0, 0]
        axis_end: [5, 0, 0]
        radius_start: 3
        radius_end: 3
        inner_radius_start: 0.0
        inner_radius_end: 0.0
        spacing_start_rel: [4e-2, 4e-2, 4e-2]
        spacing_end_rel: [4e-2, 4e-2, 4e-2]
        spacing_growth: 1.
        num_points_per_layer: 20
    radial_refinement:
      - max_spacing_rel: 10.
        num_points_per_layer: 8
        spacing_growth: 2.

  - name: wake-far
    type: telescoping
    seed: hexprism
    smoothing_layers: 8
    points_per_bbox: [8, 8, 8]
    frustum:
      - axis_start: [5, 0, 0]
        axis_end: [100, 0, 0]
        radius_start: 3
        radius_end: 6
        inner_radius_start: 0.0
        inner_radius_end: 0.0
        spacing_start_rel: [4e-2, 4e-2, 4e-2]
        spacing_end_rel: [10, 10, 4e-2]
        spacing_growth: 1.1
        num_points_per_layer: 8
    radial_refinement:
      - max_spacing_rel: 10.
        num_points_per_layer: 8
        spacing_growth: 2.

  - name: cylinder_near_field
    type: telescoping
    seed: hexprism
    smoothing_layers: 8
    points_per_bbox: [8, 8, 8]
    frustum:
      - axis_start: [0, 0, -1]
        axis_end: [0, 0, 1]
        radius_start: 3
        radius_end: 3
        inner_radius_start: 0.0
        inner_radius_end: 0.0
        spacing_start_rel: [4e-2, 4e-2, 4e-2]
        spacing_end_rel: [4e-2, 4e-2, 4e-2]
        spacing_growth: 1.
        num_points_per_layer: 50
    radial_refinement:
      - max_spacing_rel: 10.
        num_points_per_layer: 8
        spacing_growth: 2.