stl_mod.f

Go to the documentation of this file.

      MODULE stl

      Use param1, only: zero

!     Maximum number of facets that can be read
      INTEGER, PARAMETER          :: dim_stl = 10000000   !10 Million
!     Number of facets
      INTEGER                     :: n_facets
      INTEGER, PARAMETER :: west_faceid = 9000000
      INTEGER, PARAMETER :: east_faceid = 9000001
      INTEGER, PARAMETER :: south_faceid = 9000002
      INTEGER, PARAMETER :: north_faceid = 9000003
      INTEGER, PARAMETER :: bottom_faceid = 9000004
      INTEGER, PARAMETER :: top_faceid = 9000005

!     Number of facets for des. This could be a diiferent number from
! N_FACETS if the outer boundary is triangulated here
      INTEGER                     :: n_facets_des
!     Vertex Coordinates X ,Y and Z
      DOUBLE PRECISION, DIMENSION(3,3,DIM_STL) :: vertex
!     Face normal vector (normalized)
      DOUBLE PRECISION, DIMENSION(3,DIM_STL) :: norm_face
!     TRANSLATION COMPONENTS
      DOUBLE PRECISION :: tx_stl,ty_stl,tz_stl
      DOUBLE PRECISION :: tx_msh,ty_msh,tz_msh
!     SCALING FACTOR
      DOUBLE PRECISION :: scale_stl
      DOUBLE PRECISION :: scale_msh
!     RANGE OF STL FILE:
      DOUBLE PRECISION :: xmin_stl,xmax_stl
      DOUBLE PRECISION :: ymin_stl,ymax_stl
      DOUBLE PRECISION :: zmin_stl,zmax_stl
      DOUBLE PRECISION :: xmin_msh,xmax_msh
      DOUBLE PRECISION :: ymin_msh,ymax_msh
      DOUBLE PRECISION :: zmin_msh,zmax_msh
!     VALUE OF F_STL OUTSIDE OF STL BOUNDING BOX
      DOUBLE PRECISION :: out_stl_value
      DOUBLE PRECISION :: out_msh_value
!     SMALLEST ANGLE FOR DETECTION OF SMALL TRIANGLES
      DOUBLE PRECISION :: stl_small_angle
      DOUBLE PRECISION :: msh_small_angle
!     DIRECTION OF RAY TRACED TO DETERMINE WHETHER A POINT IS INSIDE/OUTSIDE
      CHARACTER(LEN=3) :: ray_dir
!     Tolerance for polygone edge detection
      DOUBLE PRECISION :: tol_stl
      DOUBLE PRECISION :: tol_msh
      DOUBLE PRECISION :: tol_stl_dp
!     Boundary condition ID
      INTEGER :: stl_bc_id

      INTEGER, DIMENSION(DIM_STL) :: bc_id_stl_face

!     Maximum number of facets per cell. The arrays below are used
! to define cut-cells under the CG modules
      INTEGER          :: dim_facets_per_cell
      INTEGER, DIMENSION (:), ALLOCATABLE ::  n_facet_at
      INTEGER, DIMENSION (:,:), ALLOCATABLE ::  list_facet_at

! List of facets intersecting each DES grid cell
      TYPE facets_to_dg
         INTEGER :: count
         INTEGER, ALLOCATABLE :: id(:)
         INTEGER, ALLOCATABLE :: dir(:)
         DOUBLE PRECISION, ALLOCATABLE :: min(:)
         DOUBLE PRECISION, ALLOCATABLE :: max(:)
      END TYPE facets_to_dg

! Start/End positions for the different kinds of STLs
      INTEGER, PARAMETER :: all_stl = 0
      INTEGER, PARAMETER :: base_stl = 1
      INTEGER, PARAMETER :: bcwalls_stl = 2
      INTEGER, PARAMETER :: imprmbl_stl = 3
      INTEGER, PARAMETER :: default_stl = 4
      INTEGER :: stl_start(4), stl_end(4)

      type(facets_to_dg), ALLOCATABLE ::  facets_at_dg(:)
      CHARACTER(LEN=3) :: cad_propagate_order

      END MODULE stl