set_index1a3.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: SET_INDEX1A3(I, J, K, IJK, IMJK, IPJK, IJMK, IJPK,      C
!                            IJKM, IJKP, IJKW, IJKE, IJKS, IJKN,       C
!                            IJKB, IJKT)                               C
!  Purpose: Set the indices of the neighbors of cell ijk (brute force) C
!                                                                      C
!  Author: M. Syamlal                                 Date: 21-JAN-92  C
!  Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti,  Date: 29-JAN-92  C
!           W. Rogers                                                  C
!                                                                      C
!  Revision Number: 1                                                  C
!  Purpose: Modify index computations for K for setting periodic       C
!           boundary conditions in a cylindrical geometry where z goes C
!           from 0 to 2 pi                                             C
!  Author: M. Syamlal                                 Date: 10-MAR-92  C
!  Revision Number: 2                                                  C
!  Purpose:  Calculate only the nearest neighbor indices.( for code    C
!            optimization)                                             C
!  Author: M. Syamlal                                 Date: 23-SEP-92  C
!  Reviewer: M. Syamlal                               Date: 11-DEC-92  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced: I, J, K, IJK                                  C
!                                                                      C
!  Variables modified: IJKM, IJMK, IMJK, IPJK, IJPK, IJKP, IJKW, IJKE, C
!                      IJKS, IJKN, IJKB, IJKT                          C
!                                                                      C
!  Local variables: None                                               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE set_index1a3(I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, &
         ijkw, ijke, ijks, ijkn, ijkb, ijkt)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE physprop
      USE geometry
      USE compar
      USE fldvar
      USE indices
      USE functions
      USE function3
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
      INTEGER I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, IJKW, IJKE, &
         IJKS, IJKN, IJKB, IJKT
!-----------------------------------------------
!   L o c a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------

      imjk = undefined_i
      ipjk = undefined_i
      ijmk = undefined_i
      ijpk = undefined_i
      ijkm = undefined_i
      ijkp = undefined_i


      IF(im1_3(i).NE.undefined_i) THEN
        imjk = bound_funijk3(im1_3(i),j,k)
      ENDIF

      IF(ip1_3(i).NE.undefined_i) THEN
        ipjk = bound_funijk3(ip1_3(i),j,k)
      ENDIF

      IF(jm1_3(j).NE.undefined_i) THEN
        ijmk = bound_funijk3(i,jm1_3(j),k)
      ENDIF

      IF(jp1_3(j).NE.undefined_i) THEN
        ijpk = bound_funijk3(i,jp1_3(j),k)
      ENDIF

      IF(km1_3(k).NE.undefined_i) THEN
        ijkm = bound_funijk3(i,j,km1_3(k))
      ENDIF

      IF(kp1_3(k).NE.undefined_i) THEN
        ijkp = bound_funijk3(i,j,kp1_3(k))
      ENDIF
!
      RETURN
      END SUBROUTINE set_index1a3

!// Comments on the modifications for DMP version implementation
!// Modified calls to BOUND_FUNIJK to have a self consistent formulation