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