1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: SET_INDEX1(IJK, I, J, K, IMJK, IPJK, IJMK, IJPK, C 4 ! IJKM, IJKP, IJKW, IJKE, IJKS, IJKN, C 5 ! IJKB, IJKT, IM, JM, KM) C 6 ! Purpose: Set the indices of the first neighbors of cell ijk C 7 ! This version adds 'increments' stored in STORE_INCREMENTS C 8 ! to IJK to find indices of the neighbors. C 9 ! C 10 ! Author: M. Syamlal, W. A. Rogers Date: 17-Dec-91 C 11 ! Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti, Date: 29-JAN-92 C 12 ! W. Rogers C 13 ! C 14 ! Revision Number: 1 C 15 ! Purpose: Second index of store_increments changed to PARAMETER's. C 16 ! Use IM_OF etc. instead of BOUND_FUNIJK for computing C 17 ! IMJK etc. C 18 ! Author: M. Syamlal Date: 18-FEB-92 C 19 ! Revision Number:2 C 20 ! Purpose: change STORE_INCREMENTS to INCREMENT_FOR_xx. Do only C 21 ! calculation of the nearest neighbors. Remove MIN and MAX C 22 ! from IM, IP etc. calculations. C 23 ! Author: M. Syamlal Date: 18-SEP-92 C 24 ! Reviewer: M. Syamlal Date: 11-DEC-92 C 25 ! C 26 ! Literature/Document References: C 27 ! C 28 ! Variables referenced: CELL_CLASS, INCREMENTS_FOR_xx, IJK, I, J, K C 29 ! C 30 ! Variables modified: IJKN,IJKS,IJKE,IJKW,IJKT,IJKB, IJKM, IJMK, IMJK,C 31 ! IPJK, IJPK, IJKP C 32 ! C 33 ! Local variables: C 34 ! C 35 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 36 ! 37 SUBROUTINE SET_INDEX1(IJK, I, J, K, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, & 38 IJKW, IJKE, IJKS, IJKN, IJKB, IJKT, IM, JM, KM) 39 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 40 !...Switches: -xf 41 ! 42 ! Include param.inc file to specify parameter values 43 ! 44 !----------------------------------------------- 45 ! M o d u l e s 46 !----------------------------------------------- 47 USE param 48 USE param1 49 USE physprop 50 USE fldvar 51 USE geometry 52 USE constant 53 USE indices 54 USE compar 55 USE functions 56 IMPLICIT NONE 57 !----------------------------------------------- 58 ! G l o b a l P a r a m e t e r s 59 !----------------------------------------------- 60 !----------------------------------------------- 61 ! D u m m y A r g u m e n t s 62 !----------------------------------------------- 63 ! 64 ! Indices 65 INTEGER I, J, K, IJK, IMJK, IPJK, IJMK, IJPK, IJKM, IJKP, & 66 IJKW, IJKE, IJKS, IJKN, IJKB, IJKT, & 67 IM, JM, KM 68 ! 69 !----------------------------------------------- 70 71 I = I_OF(IJK) 72 J = J_OF(IJK) 73 K = K_OF(IJK) 74 ! 75 IM = IM1(I) 76 JM = JM1(J) 77 KM = KM1(K) 78 ! 79 ! 80 ! Determine the true indices of neighboring cells 81 ! 82 ! 83 IJKW = WEST_OF(IJK) 84 IJKE = EAST_OF(IJK) 85 IJKS = SOUTH_OF(IJK) 86 IJKN = NORTH_OF(IJK) 87 IJKB = BOTTOM_OF(IJK) 88 IJKT = TOP_OF(IJK) 89 IMJK = IM_OF(IJK) 90 IPJK = IP_OF(IJK) 91 IJMK = JM_OF(IJK) 92 IJPK = JP_OF(IJK) 93 IJKM = KM_OF(IJK) 94 IJKP = KP_OF(IJK) 95 ! 96 RETURN 97 END SUBROUTINE SET_INDEX1 98 99 !// Comments on the modifications for DMP version implementation 100 !// 001 Include header file and common declarations for parallelization 101