blockPointers.F90

Go to the documentation of this file.

module blockpointers
    !
    !       This module contains the pointers for all variables inside a
    !       block. The pointers are set via the subroutine setPointers,
    !       which can be found in the utils directory. In this way the
    !       code becomes much more readable. The relation to the original
    !       multiblock grid is not copied, because it does not affect the
    !       computation.
    !       See the module block for the meaning of the variables.
    !       Note that the dimensions are not pointers, but integers.
    !       Consequently changing dimensions of a block must be done only
    !       with the variables of floDoms.
    !
    use constants, only: inttype, realtype, portype
    use block, only: fringetype, bcdatatype, viscsubfacetype, flowdoms, ndom
#ifndef USE_TAPENADE
    use block, only: flowdomsd
#endif
    implicit none
    !
    !       Additional info, such that it is known to which block the data
    !       inside this module belongs.
    !
    ! sectionID:   the section to which this block belongs.
    ! nbkLocal :   local block number.
    ! nbkGlobal:   global block number in the original cgns grid.
    ! mgLevel:     the multigrid level.
    ! spectralSol: the spectral solution index of this block.
 
    integer(kind=intType) :: sectionid
    integer(kind=intType) :: nbklocal, nbkglobal, mglevel
    integer(kind=intType) :: spectralsol
    !
    !       Variables, which are either copied or the pointer is set to
    !       the correct variable in the block. See the module block for
    !       meaning of the variables.
    !
    integer(kind=intType) :: nx, ny, nz, il, jl, kl
    integer(kind=intType) :: ie, je, ke, ib, jb, kb
    integer(kind=intType) :: maxdim, imaxdim, jmaxdim
 
    logical :: righthanded
 
    integer(kind=intType) :: ibegor, iendor, jbegor, jendor
    integer(kind=intType) :: kbegor, kendor
 
    integer(kind=intType) :: nsubface, n1to1, nbocos, nviscbocos
 
    integer(kind=intType), dimension(:), pointer :: bctype
    integer(kind=intType), dimension(:), pointer :: bcfaceid
 
    integer(kind=intType), dimension(:), pointer :: cgnssubface
 
    integer(kind=intType), dimension(:), pointer :: inbeg, inend
    integer(kind=intType), dimension(:), pointer :: jnbeg, jnend
    integer(kind=intType), dimension(:), pointer :: knbeg, knend
 
    integer(kind=intType), dimension(:), pointer :: dinbeg, dinend
    integer(kind=intType), dimension(:), pointer :: djnbeg, djnend
    integer(kind=intType), dimension(:), pointer :: dknbeg, dknend
 
    integer(kind=intType), dimension(:), pointer :: icbeg, icend
    integer(kind=intType), dimension(:), pointer :: jcbeg, jcend
    integer(kind=intType), dimension(:), pointer :: kcbeg, kcend
 
    integer(kind=intType), dimension(:), pointer :: neighblock
    integer(kind=intType), dimension(:), pointer :: neighproc
    integer(kind=intType), dimension(:), pointer :: l1, l2, l3
    integer(kind=intType), dimension(:), pointer :: groupnum
 
    integer(kind=intType), dimension(:, :, :), pointer :: iblank
    integer(kind=intType), dimension(:, :, :), pointer :: status
    integer(kind=intType), dimension(:, :, :), pointer :: forcedrecv
    type(fringetype), dimension(:), pointer :: fringes
    integer(kind=intType), pointer :: ndonors
    integer(kind=intType), dimension(:, :, :, :), pointer :: fringeptr
    integer(kind=intType), dimension(:, :, :, :), pointer :: gind
    integer(kind=intType), dimension(:, :), pointer :: orphans
    integer(kind=intType) :: norphans
 
    integer(kind=intType), dimension(:), pointer :: neighblockover
    integer(kind=intType), dimension(:), pointer :: neighprocover
 
    type(bcdatatype), dimension(:), pointer :: bcdata
    type(viscsubfacetype), dimension(:), pointer :: viscsubface
 
    integer(kind=intType), dimension(:, :), pointer :: visciminpointer
    integer(kind=intType), dimension(:, :), pointer :: viscimaxpointer
    integer(kind=intType), dimension(:, :), pointer :: viscjminpointer
    integer(kind=intType), dimension(:, :), pointer :: viscjmaxpointer
    integer(kind=intType), dimension(:, :), pointer :: visckminpointer
    integer(kind=intType), dimension(:, :), pointer :: visckmaxpointer
 
    real(kind=realtype), dimension(:, :, :, :), pointer :: x
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: xold
    real(kind=realtype), dimension(:, :, :, :), pointer :: si, sj, sk
    real(kind=realtype), dimension(:, :, :), pointer :: vol
    real(kind=realtype), dimension(:, :, :), pointer :: volref
    real(kind=realtype), dimension(:, :, :, :), pointer :: volold
    real(kind=realtype), dimension(:, :, :), pointer :: skew
    real(kind=realtype), dimension(:, :, :, :), pointer :: dadidata
 
    integer(kind=porType), dimension(:, :, :), pointer :: pori, porj, pork
 
    integer(kind=intType), dimension(:, :, :), pointer :: indfamilyi
    integer(kind=intType), dimension(:, :, :), pointer :: indfamilyj
    integer(kind=intType), dimension(:, :, :), pointer :: indfamilyk
 
    integer(kind=intType), dimension(:, :, :), pointer :: factfamilyi
    integer(kind=intType), dimension(:, :, :), pointer :: factfamilyj
    integer(kind=intType), dimension(:, :, :), pointer :: factfamilyk
 
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: rotmatrixi
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: rotmatrixj
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: rotmatrixk
 
    logical :: blockismoving, addgridvelocities
 
    real(kind=realtype), dimension(:, :, :), pointer :: sfacei, sfacej, sfacek
    real(kind=realtype), dimension(:, :, :, :), pointer :: w
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: wold
 
    real(kind=realtype), dimension(:, :, :), pointer :: p, gamma, aa
    real(kind=realtype), dimension(:, :, :), pointer :: shocksensor
    real(kind=realtype), dimension(:, :, :), pointer :: rlv, rev
    real(kind=realtype), dimension(:, :, :, :), pointer :: s
    real(kind=realtype), dimension(:, :, :), pointer :: p1
    real(kind=realtype), dimension(:, :, :, :), pointer :: dw, fw
    real(kind=realtype), dimension(:, :, :, :), pointer :: scratch
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: dwoldrk
    real(kind=realtype), dimension(:, :, :, :), pointer :: w1, wr
    real(kind=realtype), dimension(:, :, :), pointer :: ux, uy, uz
    real(kind=realtype), dimension(:, :, :), pointer :: vx, vy, vz
    real(kind=realtype), dimension(:, :, :), pointer :: wx, wy, wz
    real(kind=realtype), dimension(:, :, :), pointer :: qx, qy, qz
 
    integer(kind=intType), dimension(:, :), pointer :: mgifine
    integer(kind=intType), dimension(:, :), pointer :: mgjfine
    integer(kind=intType), dimension(:, :), pointer :: mgkfine
 
    real(kind=realtype), dimension(:), pointer :: mgiweight
    real(kind=realtype), dimension(:), pointer :: mgjweight
    real(kind=realtype), dimension(:), pointer :: mgkweight
 
    integer(kind=intType), dimension(:, :), pointer :: mgicoarse
    integer(kind=intType), dimension(:, :), pointer :: mgjcoarse
    integer(kind=intType), dimension(:, :), pointer :: mgkcoarse
 
    real(kind=realtype), dimension(:, :, :, :), pointer :: wn
    real(kind=realtype), dimension(:, :, :), pointer :: pn
    real(kind=realtype), dimension(:, :, :), pointer :: dtl
    real(kind=realtype), dimension(:, :, :), pointer :: radi, radj, radk
 
    real(kind=realtype), dimension(:, :, :), pointer :: d2wall
    real(kind=realtype), dimension(:, :, :), pointer :: intermittency
    real(kind=realtype), dimension(:, :, :), pointer :: filterdes  ! eran-des
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmti1
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmti2
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmtj1
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmtj2
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmtk1
    real(kind=realtype), dimension(:, :, :, :), pointer :: bmtk2
    real(kind=realtype), dimension(:, :, :), pointer :: bvti1, bvti2
    real(kind=realtype), dimension(:, :, :), pointer :: bvtj1, bvtj2
    real(kind=realtype), dimension(:, :, :), pointer :: bvtk1, bvtk2
 
    integer(kind=intType), dimension(:, :, :), pointer :: globalnode
    integer(kind=intType), dimension(:, :, :), pointer :: globalcell
    real(kind=realtype), dimension(:, :, :, :), pointer :: xseed
    integer(kind=intType), dimension(:, :, :), pointer :: wallind
 
    real(kind=realtype), dimension(:, :, :, :), pointer :: w_offtimeinstance
    real(kind=realtype), dimension(:, :, :), pointer :: vol_offtimeinstance
 
    ! Added by HDN
    real(kind=realtype), dimension(:, :, :, :), pointer :: xale
    real(kind=realtype), dimension(:, :, :, :), pointer :: sveloiale, svelojale, svelokale
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: siale, sjale, skale
    real(kind=realtype), dimension(:, :, :, :), pointer :: sfaceiale, sfacejale, sfacekale
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: dwale, fwale
 
#ifndef USE_TAPENADE
    TYPE(viscsubfacetype), DIMENSION(:), POINTER :: viscsubfaced
 
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: xd
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: sid, sjd, skd
 
    real(kind=realtype), dimension(:, :, :), pointer :: vold
 
    REAL(kind=realtype), DIMENSION(:, :, :, :, :), POINTER :: rotmatrixid
    REAL(kind=realtype), DIMENSION(:, :, :, :, :), POINTER :: rotmatrixjd
    REAL(kind=realtype), DIMENSION(:, :, :, :, :), POINTER :: rotmatrixkd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: sfaceid
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: sfacejd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: sfacekd
 
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: wd
    REAL(kind=realtype), DIMENSION(:, :, :, :, :), POINTER :: woldd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: uxd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: uyd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: uzd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: vxd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: vyd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: vzd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: wxd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: wyd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: wzd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: qxd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: qyd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: qzd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: pd, gammad, aad
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: rlvd, revd
 
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: sd
 
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: dwd, fwd
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: w1d, wrd
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: scratchd
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: dtld
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: radid
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: radjd
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: radkd
 
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmti1d
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmti2d
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmtj1d
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmtj2d
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmtk1d
    REAL(kind=realtype), DIMENSION(:, :, :, :), POINTER :: bmtk2d
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: bvti1d, bvti2d
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: bvtj1d, bvtj2d
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: bvtk1d, bvtk2d
 
    REAL(kind=realtype), DIMENSION(:, :, :), POINTER :: d2walld
 
    real(kind=realtype), dimension(:, :, :, :), pointer :: w_offtimeinstanced
    real(kind=realtype), dimension(:, :, :), pointer :: vol_offtimeinstanced
 
    type(bcdatatype), dimension(:), pointer :: bcdatad
 
    real(kind=realtype), dimension(:, :, :, :, :), pointer :: pcmat
    real(kind=realtype), dimension(:, :, :, :), pointer :: pcvec1, pcvec2
 
    real(kind=realtype), dimension(:, :, :, :), pointer :: i_d_fact, j_d_fact, k_d_fact
    real(kind=realtype), dimension(:, :, :, :), pointer :: i_l_fact, j_l_fact, k_l_fact
    real(kind=realtype), dimension(:, :, :, :), pointer :: i_u_fact, j_u_fact, k_u_fact
    real(kind=realtype), dimension(:, :, :, :), pointer :: i_u2_fact, j_u2_fact, k_u2_fact
    integer(kind=intType), dimension(:, :, :, :), pointer :: i_ipiv, j_ipiv, k_ipiv
#endif
 
end module blockpointers