pic_bc_mod.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Module name: pic_bc_mod.f                                           !
!                                                                      !
!  Purpose: Common elements needed for the pic mass inflow boundary    !
!  condition.                                                          !
!                                                                      !
!  Author: R. Garg                                   Date: 11-Jun-14   !
!                                                                      !
!  Comments:                                                           !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      MODULE pic_bc

      use param, only: dimension_bc
!......................................................................!

      INTEGER :: pic_bcmi
      INTEGER :: pic_bcmo

      LOGICAL :: pic_mio  ! either inlet or outlet exists

! Map between PIC MI/MO IDs and the user input BC index.
      INTEGER :: pic_bcmi_map(dimension_bc)
      INTEGER :: pic_bcmo_map(dimension_bc)


      INTEGER, ALLOCATABLE :: pic_bcmo_ijkstart(:)
      INTEGER, ALLOCATABLE :: pic_bcmo_ijkend(:)

      INTEGER, ALLOCATABLE :: pic_bcmo_ijk(:)


      INTEGER, ALLOCATABLE :: pic_bcmi_ijkstart(:)
      INTEGER, ALLOCATABLE :: pic_bcmi_ijkend(:)

      INTEGER, ALLOCATABLE :: pic_bcmi_ijk(:)

! Direction of the MI BC plane
      INTEGER, ALLOCATABLE :: pic_bcmi_normdir(:,:)
! BC planes will be calculated from the fluid node as
! XE-dx, YN-dy, and ZT-dz.
! For BC-plane 'S', YN-dy will provide the wrong plane, as
! it shud be just yn. The array below is used to provide this correction.
! BC_plane coords will be calculated as x_i - offset dx_i and offset will
! have default value of 1. It will be set to 0 for a case described above.
      INTEGER, ALLOCATABLE :: pic_bcmi_offset(:,:)


! This will store cumulative number of computational parcels
      Double precision, Allocatable :: pic_bcmi_cnp(:, :)
! This will store cumulative number of implied real particles
      Double precision, Allocatable :: pic_bcmi_rnp(:, :)

! To include or not to include cut-cells in the MI BC

      LOGICAL, Allocatable :: pic_bcmi_incl_cutcell(:)

! Logicals to print seeding and deletion of particles based on user needs
      LOGICAL :: pic_report_seeding_stats, pic_report_deletion_stats

! Mininum velocity for parcles at rest of offset gravitational forces
      DOUBLE PRECISION :: minvel(3), minvel_mag, oominvel_mag

      END MODULE pic_bc