vtk_mod.f

Go to the documentation of this file.

      MODULE vtk

      Use param1
! Maximum number of solids phases.
      use param, only: dim_m
! Maximum number of gas phase species
      use param, only: dim_n_g
! Maximum number of solids phase species
      use param, only: dim_n_s
! Maximum number of scalar equations
      use param, only: dim_scalar
! Maximum number of DEM solids phase species
use param, only: dimension_n_s



      INTEGER number_of_cells
      INTEGER number_of_cut_cells
      INTEGER number_of_blocked_cells
      INTEGER number_of_standard_cells
      INTEGER number_of_vtk_cells
      INTEGER number_of_vtk_parts

      LOGICAL             :: write_vtk_files
      LOGICAL             :: time_dependent_filename
      LOGICAL             :: reset_frame_at_time_zero=.true.
      CHARACTER (LEN=255)  :: vtu_dir
      CHARACTER (LEN=255)  :: vtk_filename,frame_char,vtu_filename,pvd_filename,pvtu_filename
      CHARACTER (LEN=255)  :: vtu_frame_filename='VTU_FRAME_INDEX.TXT'
      CHARACTER (LEN=512) :: buffer

      CHARACTER (LEN=1), PARAMETER  :: end_rec = char(10)

      INTEGER :: boundary_unit=122
      INTEGER :: vtk_unit=123,vtu_unit=124,pvd_unit=125
      INTEGER :: pvtu_unit=126,vtu_frame_unit=127
      INTEGER, PARAMETER :: dim_vtk_var = 20
      INTEGER, DIMENSION(DIM_VTK_VAR) :: vtk_var


      INTEGER :: poly_counter,number_of_points

      integer, allocatable :: cleaned_connectivity(:,:)
      REAL, allocatable :: coords_of_points(:,:)

      integer, allocatable :: global_i_of(:)
      integer, allocatable :: global_j_of(:)
      integer, allocatable :: global_k_of(:)
      integer, allocatable :: global_connectivity(:,:)
      integer, allocatable :: global_cleaned_connectivity(:,:)
      integer, allocatable :: global_number_of_nodes(:)

      REAL, allocatable :: global_coords_of_points(:,:)

      LOGICAL, allocatable :: global_interior_cell_at(:)
      LOGICAL, allocatable :: global_blocked_cell_at(:)
      LOGICAL, allocatable :: global_standard_cell_at(:)
      LOGICAL, allocatable :: global_cut_cell_at(:)
      LOGICAL, allocatable :: global_snap(:)
      DOUBLE PRECISION, allocatable :: global_f_at(:)

      double precision, allocatable :: global_x_new_point(:)
      double precision, allocatable :: global_y_new_point(:)
      double precision, allocatable :: global_z_new_point(:)

      INTEGER :: global_number_of_new_points


      LOGICAL :: global_var_allocated

      LOGICAL :: grid_info_printed_on_screen

      LOGICAL :: write_ani_cutcell


      INTEGER :: vtu_offset

      LOGICAL, allocatable :: belongs_to_vtk_subdomain(:)
      ! LOGICAL, allocatable :: PART_BELONGS_TO_VTK_SUBDOMAIN(:)

      INTEGER, PARAMETER :: dimension_vtk = 100
! Max size for VTK_nRR
      INTEGER, PARAMETER :: vtk_nrrmax = 100

! Current VTK region
      INTEGER :: vtk_region

! Time interval at which vtk files are saved
      DOUBLE PRECISION :: vtk_dt(dimension_vtk)

! Flag to distiguish between regular file (VTK_DBG_FILE=.FALSE., default value)
!                        and debug file (VTK_DBG_FILE=.TRUE.)
      LOGICAL :: vtk_dbg_file(dimension_vtk)

! Current vtk time
      DOUBLE PRECISION :: vtk_time(dimension_vtk)

! Type of data in vtk region: 'C':cell data, 'P': particle data
      CHARACTER(LEN=1) :: vtk_data(dimension_vtk)

! FRAME index of vtk file
      INTEGER :: frame(dimension_vtk)

! PVD file initialization flag
      LOGICAL :: pvd_file_initialized(dimension_vtk)=.false.

! Logical variable to determine whether an vtk region is defined
      LOGICAL :: vtk_defined (dimension_vtk)

! VTK region West face, X-coordinate
      DOUBLE PRECISION :: vtk_x_w (dimension_vtk)

! VTK region East face, X-coordinate
      DOUBLE PRECISION :: vtk_x_e (dimension_vtk)

! VTK region South face, Y-coordinate
      DOUBLE PRECISION :: vtk_y_s (dimension_vtk)

! VTK region North face, Y-coordinate
      DOUBLE PRECISION :: vtk_y_n (dimension_vtk)

! VTK region Bottom face, Z-coordinate
      DOUBLE PRECISION :: vtk_z_b (dimension_vtk)

! VTK region Top face, Z-coordinate
      DOUBLE PRECISION :: vtk_z_t (dimension_vtk)

! VTK number of slices in x-direction
      INTEGER :: vtk_nxs(dimension_vtk)

! VTK number of slices in y-direction
      INTEGER :: vtk_nys(dimension_vtk)

! VTK number of slices in z-direction
      INTEGER :: vtk_nzs(dimension_vtk)

! VTK slice tolerance
      DOUBLE PRECISION :: vtk_slice_tol(dimension_vtk)

! Flag to write only cut cell data in VTK file
      LOGICAL :: vtk_cutcell_only(dimension_vtk)

! VTK filename base
      CHARACTER(LEN=255) :: vtk_filebase(dimension_vtk)

! Gas phase volume fraction
      LOGICAL :: vtk_ep_g (dimension_vtk)

! Gas pressure
      LOGICAL :: vtk_p_g (dimension_vtk)

! Solids pressure
      LOGICAL :: vtk_p_star(dimension_vtk)

! X-component of gas velocity
      LOGICAL :: vtk_u_g(dimension_vtk)

! X-component of solids phase velocity
      LOGICAL :: vtk_u_s(dimension_vtk, dim_m)

! Y-component of gas velocity
      LOGICAL :: vtk_v_g(dimension_vtk)

! Y-component of solids phase velocity
      LOGICAL :: vtk_v_s(dimension_vtk, dim_m)

! Z-component of gas velocity
      LOGICAL :: vtk_w_g(dimension_vtk)

! Z-component of solids phase velocity
      LOGICAL :: vtk_w_s(dimension_vtk, dim_m)

! Gas velocity vector
      LOGICAL :: vtk_vel_g(dimension_vtk)

! Solids velocity vector
      LOGICAL :: vtk_vel_s(dimension_vtk, dim_m)

! Macroscopic density of solids phases
      LOGICAL :: vtk_rop_s(dimension_vtk, dim_m)

! Solids phase volume fraction
      LOGICAL :: vtk_ep_s (dimension_vtk, dim_m)

! Gas temperature
      LOGICAL :: vtk_t_g(dimension_vtk)

! Solids temperature
      LOGICAL :: vtk_t_s(dimension_vtk, dim_m)

! Gas species mass fractions
      LOGICAL :: vtk_x_g(dimension_vtk, dim_n_g)

! Solids species mass fractions
      LOGICAL :: vtk_x_s(dimension_vtk, dim_m, dim_n_s)

! Granular temperature
      LOGICAL :: vtk_theta_m(dimension_vtk, dim_m)

! Scalar value
      LOGICAL :: vtk_scalar(dimension_vtk, dim_scalar)

! Reaction rates
      LOGICAL :: vtk_rrate(dimension_vtk,vtk_nrrmax)

! K & Epsilon values
      LOGICAL :: vtk_k_turb_g(dimension_vtk)
      LOGICAL :: vtk_e_turb_g(dimension_vtk)

! Vorticity magnitude
      LOGICAL :: vtk_vorticity(dimension_vtk)

! Lambda_2 (vortex core indicator)
      LOGICAL :: vtk_lambda_2(dimension_vtk)

! Grid partition
      LOGICAL :: vtk_partition(dimension_vtk)

! Boundary condition ID
      LOGICAL :: vtk_bc_id(dimension_vtk)

! Wall distance
      LOGICAL :: vtk_dwall(dimension_vtk)

! Facet count (DES)
      LOGICAL :: vtk_facet_count_des(dimension_vtk)

! Neighboring facets (DES)
      LOGICAL :: vtk_nb_facet_des(dimension_vtk)

! Cell IJK index
      LOGICAL :: vtk_ijk(dimension_vtk)

! Cut face normal vector
      LOGICAL :: vtk_normal(dimension_vtk)

! Debug variable
      LOGICAL :: vtk_debug(dimension_vtk,15)

      INTEGER, DIMENSION(DIMENSION_VTK,DIM_VTK_VAR) :: vtk_varlist

! VTK Particle selection: 'C': center, 'P': entire particle, 'I': particle intersect
      CHARACTER(LEN=1) :: vtk_select_mode(dimension_vtk)

! Particle radius
      LOGICAL :: vtk_part_diameter(dimension_vtk)

! Particle velocity
      LOGICAL :: vtk_part_vel(dimension_vtk)

! Particle angular velocity
      LOGICAL :: vtk_part_angular_vel(dimension_vtk)

! Particle orientation
      LOGICAL :: vtk_part_orientation(dimension_vtk)

! Particle user-defined variable
      LOGICAL :: vtk_part_usr_var(dimension_vtk,3)

! Particle temperature
      LOGICAL :: vtk_part_temp(dimension_vtk)

! Particle species mass fraction
      LOGICAL :: vtk_part_x_s(dimension_vtk,100)

! Particle cohesion
      LOGICAL :: vtk_part_cohesion(dimension_vtk)

      END MODULE vtk