cutcell_mod.f

Go to the documentation of this file.

      MODULE cutcell

      ! Use param
      ! Use param1
      USE progress_bar

!     CUT_CELL.LOG unit number
      INTEGER  unit_cut_cell_log
      parameter(unit_cut_cell_log = 111)

!     Flag to activate Cartesian grid

      LOGICAL :: cartesian_grid

!     Flag to activate cell re-indexing
      LOGICAL :: re_indexing

!     Flag to activate adjustment of each processor domain size
      LOGICAL :: adjust_proc_domain_size

!     Flag to report best domain size on each processor and quit
      LOGICAL :: report_best_domain_size

!     maximum number of cut cells
      INTEGER :: dimension_max_cut_cell

!     Factor used to allocate cut cells arrays
      DOUBLE PRECISION :: fac_dim_max_cut_cell

!     Flag to identify interior cells.
!     This flag is true for any cell within the range
!     istart1 <= i <= iend1 and so on for j and k.
!     Caution: it is local to each processor and this array
!     is not shared accross processor boundaries.
!     For example, consider 1-D decomposition over two procs.
!     on each proc, INTERIOR_CELL_AT will be true for
!     istart1 <= i <= iend1. On processor 1, it will be false
!     at i=iend2 even though i=iend2 is an interior cell in the
!     global sense!
      LOGICAL, DIMENSION(:), ALLOCATABLE :: interior_cell_at

!     One-Dimensional Arrays for East, North, Top location of
!     original (uncut) scalar cells
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::xg_e
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::yg_n
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::zg_t

!     location of U-momentum nodes
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_u
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_u
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_u

!     location of V-momentum nodes
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_v
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_v
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_v

!     location of W-momentum nodes
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_w
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_w
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_w

!     Intersection flags
      LOGICAL, DIMENSION(:), ALLOCATABLE :: intersect_x
      LOGICAL, DIMENSION(:), ALLOCATABLE :: intersect_y
      LOGICAL, DIMENSION(:), ALLOCATABLE :: intersect_z

!     Location of intersections
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_int

!     Location of original (uncut) corner cell nodes
      DOUBLE PRECISION, DIMENSION(0:15):: x_node
      DOUBLE PRECISION, DIMENSION(0:15):: y_node
      DOUBLE PRECISION, DIMENSION(0:15):: z_node
      DOUBLE PRECISION, DIMENSION(0:15):: f_node
      INTEGER, DIMENSION(0:15) :: ijk_of_node

!     Location of new (along intersecting edges) nodes
      INTEGER :: number_of_new_points
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_new_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_new_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_new_point

!     Location of new (along intersecting edges) nodes
      INTEGER :: number_of_new_u_points
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_new_u_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_new_u_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_new_u_point

!     Location of new (along intersecting edges) nodes
      INTEGER :: number_of_new_v_points
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_new_v_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_new_v_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_new_v_point

!     Location of new (along intersecting edges) nodes
      INTEGER :: number_of_new_w_points
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_new_w_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_new_w_point
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_new_w_point

!     Number of nodes
      INTEGER, DIMENSION(:), ALLOCATABLE ::  number_of_nodes
      INTEGER, DIMENSION(:), ALLOCATABLE ::  number_of_u_nodes
      INTEGER, DIMENSION(:), ALLOCATABLE ::  number_of_v_nodes
      INTEGER, DIMENSION(:), ALLOCATABLE ::  number_of_w_nodes

!     Connectivity
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: connectivity
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: connectivity_u
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: connectivity_v
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: connectivity_w

!     Processor assign to cell IJK
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: partition

!     Normal Vector Defining cut face in Scalar Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  normal_s

!     Reference point Defining cut face in Scalar Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  refp_s

!     Flags for Wall momentum cells
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  wall_u_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  wall_v_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  wall_w_at

!     Areas of cut faces
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  area_cut
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  area_u_cut
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  area_v_cut
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  area_w_cut

!     Distances from cell center to face center
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_ue
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_uw
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_un
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_us
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_ut
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_ub

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_ve
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_vw
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_vn
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_vs
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_vt
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_vb

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_we
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delx_ww
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_wn
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dely_ws
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_wt
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delz_wb

!     Location of face centers
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_u_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_u_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_u_ec

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_u_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_u_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_u_nc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_u_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_u_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_u_tc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_v_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_v_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_v_ec

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_v_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_v_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_v_nc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_v_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_v_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_v_tc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_w_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_w_ec
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_w_ec

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_w_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_w_nc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_w_nc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  x_w_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  y_w_tc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  z_w_tc

!     Distance to cut face in Scalar Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delh_scalar


!     Distance to cut face in U-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delh_u

!     Normal Vector Defining cut face in U-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  normal_u

!     Reference point Defining cut face in W-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  refp_u


!     Correction factors for U-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ue
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ue_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_u_ne
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_u_nw

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_u_te
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_u_tw

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_ue_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_u_e

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_un
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_un_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_un_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_u_n

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ut
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ut_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_ut_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_u_t

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_upg_e
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_upg_w

!     Distance to cut face in V-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delh_v

!     Normal Vector Defining cut face in V-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  normal_v

!     Reference point Defining cut face in V-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  refp_v


!     Correction factors for V-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_v_ne
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_v_se

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_vn
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_vn_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_v_nt
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_v_st

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ve
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_ve_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_ve_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_v_e

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_vn_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_v_n

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_vt
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_vt_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_vt_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_v_t

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_vpg_n
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_vpg_s


!     Distance to cut face in W-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delh_w

!     Normal Vector Defining cut face in W-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  normal_w

!     Reference point Defining cut face in W-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::  refp_w


!     Correction factors for W-Momentum Cell
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_w_te
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_w_be

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_w_tn
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_w_bn

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_wt
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_wt_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_we
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_we_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_we_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_w_e

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_wn
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  theta_wn_bar

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_wn_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_w_n

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  alpha_wt_c
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  noc_w_t

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_wpg_t
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  a_wpg_b

!     1/dx, 1/dy, 1/dz
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodx_e_u
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneody_n_u
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodz_t_u


      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodx_e_v
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneody_n_v
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodz_t_v

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodx_e_w
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneody_n_w
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  oneodz_t_w



      LOGICAL, DIMENSION(:), ALLOCATABLE :: along_domain_boundary

!     Location of intersection points
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  xn_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  xn_u_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  xn_v_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  xn_w_int

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ye_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ye_u_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ye_v_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ye_w_int

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  zt_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  zt_u_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  zt_v_int
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  zt_w_int

!     Cut cell treatment flags
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_treatment_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_u_treatment_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_v_treatment_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_w_treatment_at

!     Various cell flags
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_u_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_v_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  cut_w_cell_at

      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  small_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  blocked_cell_at
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  small_cell_flag

      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  blocked_u_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  blocked_v_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  blocked_w_cell_at

      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  standard_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  standard_u_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  standard_v_cell_at
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  standard_w_cell_at

!     Tolerance for snapping procedure
      DOUBLE PRECISION, DIMENSION(3) :: tol_snap

!     Tolerances for wall distance
      DOUBLE PRECISION :: tol_delh

!     Tolerance for merging nodes
      DOUBLE PRECISION :: tol_merge

!     Tolerance for detecting small scalar cells
      DOUBLE PRECISION :: tol_small_cell
      DOUBLE PRECISION :: tol_small_area
!     Maximum value of ALPHA correction factor
      DOUBLE PRECISION :: alpha_max

!     Flags to include effect of cut cells
      LOGICAL :: noc,noc_ug,noc_vg,noc_wg,noc_us,noc_vs,noc_ws,noc_trdg,noc_trds
      LOGICAL :: cut_tau_ug,cut_tau_vg,cut_tau_wg,cut_tau_us,cut_tau_vs,cut_tau_ws

!     pressure gradient option flag
      INTEGER :: pg_option

!     Number of cells
      INTEGER number_of_u_cut_cells
      INTEGER number_of_v_cut_cells
      INTEGER number_of_w_cut_cells
      INTEGER number_of_small_cells

      INTEGER number_of_u_wall_cells
      INTEGER number_of_v_wall_cells
      INTEGER number_of_w_wall_cells

!     Vorticity and lambda2
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  vorticity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  lambda2

!     Re-ordering array
      INTEGER, DIMENSION(15) :: order

!     Snapping flag
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  snap

      INTEGER, DIMENSION(10) :: cg_safe_mode
      LOGICAL :: print_warnings
      LOGICAL :: set_corner_cells

!     Master cell of wall cell (FSW)
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  u_master_of
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  v_master_of
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  w_master_of

      INTEGER :: n_usr_def

      LOGICAL :: use_polygon

      LOGICAL :: use_stl

      LOGICAL :: use_msh

!     Boundary condition flag
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  bc_id
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  bc_u_id
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  bc_v_id
      INTEGER,  DIMENSION(:), ALLOCATABLE ::  bc_w_id

      INTEGER :: nsw_ghost_bc_id

!     Under-relaxation flag applied to cut cells
      DOUBLE PRECISION, DIMENSION(9):: cg_ur_fac

!     Debugging_variables
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::debug_cg

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::u_g_cc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::v_g_cc
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::w_g_cc

      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::u_s_cc
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::v_s_cc
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::w_s_cc

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  trd_g_out
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  pp_g_out
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  epp_out

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dudx_out
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dvdy_out
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  delv_out

      LOGICAL :: partial_check_03

!     Keep track of whether the CG Hader was printed on screen
      LOGICAL :: cg_header_was_printed = .false.

      LOGICAL, DIMENSION(:), ALLOCATABLE ::scalar_node_atwall

      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE ::scalar_node_xyz

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::ovol_around_node

!     Snapping flag
      LOGICAL,  DIMENSION(:), ALLOCATABLE ::  potential_cut_cell_at

      INTEGER, PARAMETER :: max_cp = 100

      DOUBLE PRECISION, DIMENSION(0:MAX_CP) ::  cpx,cpy,cpz  ! Control point location
      DOUBLE PRECISION, DIMENSION(MAX_CP)   ::  erx,ery,erz  ! Expansion Ratio
      INTEGER, DIMENSION(MAX_CP)             ::  ncx,ncy,ncz  ! Number of cell in a segment
      DOUBLE PRECISION, DIMENSION(MAX_CP)   ::  first_dx,last_dx! DX values at segment extremities
      DOUBLE PRECISION, DIMENSION(MAX_CP)   ::  first_dy,last_dy! DY values at segment extremities
      DOUBLE PRECISION, DIMENSION(MAX_CP)   ::  first_dz,last_dz! DZ values at segment extremities


      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  f_at

      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  dwall
! Grid partition for reporting best size
      INTEGER :: nodesi_report,nodesj_report,nodesk_report

      LOGICAL :: minimize_send_recv

      LOGICAL :: dwall_brute_force

      END MODULE cutcell