1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: COPY_A_g(A_VEL, A_m, IER) C 4 ! Purpose: Copy A_VEL_g to A_m C 5 ! C 6 ! Author: M. Syamlal Date: 24-DEC-96 C 7 ! Reviewer: Date: C 8 ! C 9 ! C 10 ! Literature/Document References: C 11 ! C 12 ! Variables referenced: C 13 ! Variables modified: C 14 ! C 15 ! Local variables: C 16 ! C 17 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 18 ! 19 SUBROUTINE COPY_A_G(A_VEL, A_M) 20 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 21 !...Switches: -xf 22 ! 23 ! Include param.inc file to specify parameter values 24 ! 25 !----------------------------------------------- 26 ! M o d u l e s 27 !----------------------------------------------- 28 USE param 29 USE param1 30 USE parallel 31 USE geometry 32 USE indices 33 USE compar 34 USE functions 35 IMPLICIT NONE 36 !----------------------------------------------- 37 ! G l o b a l P a r a m e t e r s 38 !----------------------------------------------- 39 !----------------------------------------------- 40 ! D u m m y A r g u m e n t s 41 !----------------------------------------------- 42 ! 43 ! Stored coefficients 44 DOUBLE PRECISION A_VEL(DIMENSION_3, -3:3) 45 ! 46 ! Septadiagonal matrix A_m 47 DOUBLE PRECISION A_m(DIMENSION_3, -3:3, 0:DIMENSION_M) 48 ! 49 ! Indices 50 INTEGER IJK 51 !----------------------------------------------- 52 53 IJK = 1 54 IF (IJKMAX2 > 0) THEN 55 A_M(:IJKMAX2,west,0) = A_VEL(:IJKMAX2,west) 56 A_M(:IJKMAX2,east,0) = A_VEL(:IJKMAX2,east) 57 A_M(:IJKMAX2,south,0) = A_VEL(:IJKMAX2,south) 58 A_M(:IJKMAX2,north,0) = A_VEL(:IJKMAX2,north) 59 A_M(:IJKMAX2,bottom,0) = A_VEL(:IJKMAX2,bottom) 60 A_M(:IJKMAX2,top,0) = A_VEL(:IJKMAX2,top) 61 IJK = IJKMAX2 + 1 62 ENDIF 63 RETURN 64 END SUBROUTINE COPY_A_G 65 ! 66 ! 67 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 68 ! C 69 ! Module name: COPY_A_s(A_VEL, A_m, M, IER) C 70 ! Purpose: Copy A_VEL_s to A_m C 71 ! C 72 ! Author: M. Syamlal Date: 24-DEC-96 C 73 ! Reviewer: Date: C 74 ! C 75 ! C 76 ! Literature/Document References: C 77 ! C 78 ! Variables referenced: C 79 ! Variables modified: C 80 ! C 81 ! Local variables: C 82 ! C 83 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 84 ! 85 SUBROUTINE COPY_A_S(A_VEL, A_M, M) 86 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 87 !...Switches: -xf 88 ! 89 ! Include param.inc file to specify parameter values 90 ! 91 !----------------------------------------------- 92 ! M o d u l e s 93 !----------------------------------------------- 94 USE param 95 USE param1 96 USE parallel 97 USE geometry 98 USE indices 99 USE physprop 100 USE compar 101 USE functions 102 IMPLICIT NONE 103 !----------------------------------------------- 104 ! G l o b a l P a r a m e t e r s 105 !----------------------------------------------- 106 !----------------------------------------------- 107 ! D u m m y A r g u m e n t s 108 !----------------------------------------------- 109 ! 110 ! Stored coefficients 111 DOUBLE PRECISION A_VEL(DIMENSION_3, -3:3, DIMENSION_M) 112 ! 113 ! Septadiagonal matrix A_m 114 DOUBLE PRECISION A_m(DIMENSION_3, -3:3, 0:DIMENSION_M) 115 ! 116 ! Indices 117 INTEGER IJK, M 118 !----------------------------------------------- 119 120 IJK = 1 121 IF (IJKMAX2 > 0) THEN 122 A_M(:IJKMAX2,west,M) = A_VEL(:IJKMAX2,west,M) 123 A_M(:IJKMAX2,east,M) = A_VEL(:IJKMAX2,east,M) 124 A_M(:IJKMAX2,south,M) = A_VEL(:IJKMAX2,south,M) 125 A_M(:IJKMAX2,north,M) = A_VEL(:IJKMAX2,north,M) 126 A_M(:IJKMAX2,bottom,M) = A_VEL(:IJKMAX2,bottom,M) 127 A_M(:IJKMAX2,top,M) = A_VEL(:IJKMAX2,top,M) 128 IJK = IJKMAX2 + 1 129 ENDIF 130 RETURN 131 END SUBROUTINE COPY_A_S 132 133 !// Comments on the modifications for DMP version implementation 134 !// 001 Include header file and common declarations for parallelization 135