indices_mod.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Module: indices                                                     !
!  Author: M. Syamlal                                 Date: dd-mmm-yy  !
!                                                                      !
!  Purpose: Global arrays for index computations.                      !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      MODULE indices

      Use param1, only: max_class

! Increments used for index computation
      INTEGER :: increment_for_n (max_class)
      INTEGER :: increment_for_s (max_class)
      INTEGER :: increment_for_e (max_class)
      INTEGER :: increment_for_w (max_class)
      INTEGER :: increment_for_t (max_class)
      INTEGER :: increment_for_b (max_class)
      INTEGER :: increment_for_im(max_class)
      INTEGER :: increment_for_ip(max_class)
      INTEGER :: increment_for_jm(max_class)
      INTEGER :: increment_for_jp(max_class)
      INTEGER :: increment_for_km(max_class)
      INTEGER :: increment_for_kp(max_class)

      INTEGER :: increment_for_mp(6,max_class) = 0
      INTEGER :: increment_for_nb(6,max_class) = 0

! Increments used for index computation of 3rd layer
      INTEGER :: increment3_for_im(max_class)
      INTEGER :: increment3_for_ip(max_class)
      INTEGER :: increment3_for_jm(max_class)
      INTEGER :: increment3_for_jp(max_class)
      INTEGER :: increment3_for_km(max_class)
      INTEGER :: increment3_for_kp(max_class)

! Store LM index values
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: store_lm

! Identification of the cell class
      INTEGER, DIMENSION(:), ALLOCATABLE :: cell_class

! Maps between IJK and the base I/J/K indices
      INTEGER, DIMENSION(:), ALLOCATABLE :: i_of ! IJK --> I
      INTEGER, DIMENSION(:), ALLOCATABLE :: j_of ! IJK --> J
      INTEGER, DIMENSION(:), ALLOCATABLE :: k_of ! IJK --> K

! +/- increments, shifted for cyclic BCs and DMP partitions.
      INTEGER, DIMENSION(:), ALLOCATABLE :: im1, ip1 ! I-1, I+1
      INTEGER, DIMENSION(:), ALLOCATABLE :: jm1, jp1 ! J-1, J+1
      INTEGER, DIMENSION(:), ALLOCATABLE :: km1, kp1 ! K-1, K+1

! Identification of the cell class for higher order scheme
      INTEGER, DIMENSION(:), ALLOCATABLE :: cell_class3

      INTEGER, DIMENSION(:), ALLOCATABLE :: i3_of
      INTEGER, DIMENSION(:), ALLOCATABLE :: j3_of
      INTEGER, DIMENSION(:), ALLOCATABLE :: k3_of

      INTEGER, DIMENSION(:), ALLOCATABLE :: im1_3, ip1_3
      INTEGER, DIMENSION(:), ALLOCATABLE :: jm1_3, jp1_3
      INTEGER, DIMENSION(:), ALLOCATABLE :: km1_3, kp1_3

! Save original IJK value of Background grid (new to old mapping)
      INTEGER, DIMENSION(:), ALLOCATABLE :: background_ijk_of
! Save new IJK value of Background grid (old to new mapping)
      INTEGER, DIMENSION(:), ALLOCATABLE :: ijk_of_background
! Save original IJKEND3 value of Background grid
      INTEGER :: background_ijkend3

      END MODULE indices