1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: SET_INDEX1A3(I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, C 4 ! IJKM, IJKP, IJKW, IJKE, IJKS, IJKN, C 5 ! IJKB, IJKT) C 6 ! Purpose: Set the indices of the neighbors of cell ijk (brute force) C 7 ! C 8 ! Author: M. Syamlal Date: 21-JAN-92 C 9 ! Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti, Date: 29-JAN-92 C 10 ! W. Rogers C 11 ! C 12 ! Revision Number: 1 C 13 ! Purpose: Modify index computations for K for setting periodic C 14 ! boundary conditions in a cylindrical geometry where z goes C 15 ! from 0 to 2 pi C 16 ! Author: M. Syamlal Date: 10-MAR-92 C 17 ! Revision Number: 2 C 18 ! Purpose: Calculate only the nearest neighbor indices.( for code C 19 ! optimization) C 20 ! Author: M. Syamlal Date: 23-SEP-92 C 21 ! Reviewer: M. Syamlal Date: 11-DEC-92 C 22 ! C 23 ! Literature/Document References: C 24 ! C 25 ! Variables referenced: I, J, K, IJK C 26 ! C 27 ! Variables modified: IJKM, IJMK, IMJK, IPJK, IJPK, IJKP, IJKW, IJKE, C 28 ! IJKS, IJKN, IJKB, IJKT C 29 ! C 30 ! Local variables: None C 31 ! C 32 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 33 ! 34 SUBROUTINE SET_INDEX1A3(I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, & 35 IJKW, IJKE, IJKS, IJKN, IJKB, IJKT) 36 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 37 !...Switches: -xf 38 ! 39 ! Include param.inc file to specify parameter values 40 ! 41 !----------------------------------------------- 42 ! M o d u l e s 43 !----------------------------------------------- 44 USE param 45 USE param1 46 USE physprop 47 USE geometry 48 USE compar 49 USE fldvar 50 USE indices 51 USE boundfunijk3 52 USE functions 53 IMPLICIT NONE 54 !----------------------------------------------- 55 ! G l o b a l P a r a m e t e r s 56 !----------------------------------------------- 57 !----------------------------------------------- 58 ! D u m m y A r g u m e n t s 59 !----------------------------------------------- 60 INTEGER I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, IJKW, IJKE, & 61 IJKS, IJKN, IJKB, IJKT 62 !----------------------------------------------- 63 ! L o c a l P a r a m e t e r s 64 !----------------------------------------------- 65 !----------------------------------------------- 66 67 IMJK = UNDEFINED_I 68 IPJK = UNDEFINED_I 69 IJMK = UNDEFINED_I 70 IJPK = UNDEFINED_I 71 IJKM = UNDEFINED_I 72 IJKP = UNDEFINED_I 73 74 75 IF(IM1_3(I).NE.UNDEFINED_I) THEN 76 IMJK = BOUND_FUNIJK3(IM1_3(I),J,K) 77 ENDIF 78 79 IF(IP1_3(I).NE.UNDEFINED_I) THEN 80 IPJK = BOUND_FUNIJK3(IP1_3(I),J,K) 81 ENDIF 82 83 IF(JM1_3(J).NE.UNDEFINED_I) THEN 84 IJMK = BOUND_FUNIJK3(I,JM1_3(J),K) 85 ENDIF 86 87 IF(JP1_3(J).NE.UNDEFINED_I) THEN 88 IJPK = BOUND_FUNIJK3(I,JP1_3(J),K) 89 ENDIF 90 91 IF(KM1_3(K).NE.UNDEFINED_I) THEN 92 IJKM = BOUND_FUNIJK3(I,J,KM1_3(K)) 93 ENDIF 94 95 IF(KP1_3(K).NE.UNDEFINED_I) THEN 96 IJKP = BOUND_FUNIJK3(I,J,KP1_3(K)) 97 ENDIF 98 ! 99 RETURN 100 END SUBROUTINE SET_INDEX1A3 101 102 !// Comments on the modifications for DMP version implementation 103 !// Modified calls to BOUND_FUNIJK to have a self consistent formulation 104