1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: Init_Ab_m(A_m, b_m, IJKMAX2, M, IER) C C 4 ! Purpose:Initialiize the sparse matrix coefficients and the C 5 ! source vector. C 6 ! C 7 ! C 8 ! Author: M. Syamlal Date: 16-MAY-96 C 9 ! Reviewer: Date: C 10 ! C 11 ! C 12 ! Literature/Document References: C 13 ! C 14 ! Variables referenced: C 15 ! Variables modified: C 16 ! C 17 ! Local variables: C 18 ! C 19 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 20 SUBROUTINE INIT_AB_M(A_M, B_M, IJKMAX2A, M) 21 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 22 !...Switches: -xf 23 ! 24 ! Include param.inc file to specify parameter values 25 ! 26 !----------------------------------------------- 27 ! M o d u l e s 28 !----------------------------------------------- 29 USE param 30 USE param1 31 USE parallel 32 USE compar 33 IMPLICIT NONE 34 !----------------------------------------------- 35 ! G l o b a l P a r a m e t e r s 36 !----------------------------------------------- 37 !----------------------------------------------- 38 ! D u m m y A r g u m e n t s 39 !----------------------------------------------- 40 ! 41 ! Phase index 42 INTEGER M 43 ! 44 ! Maximum dimension 45 INTEGER IJKMAX2A 46 ! 47 ! Septadiagonal matrix A_m 48 DOUBLE PRECISION A_m(DIMENSION_3, -3:3, 0:DIMENSION_M) 49 ! 50 ! Source vector 51 DOUBLE PRECISION b_m(DIMENSION_3, 0:DIMENSION_M) 52 ! 53 !----------------------------------------------- 54 ! 55 ! IJK = 1 56 IF (IJKMAX2A > 0) THEN 57 !$omp parallel 58 !$omp sections 59 A_M(:,bottom,M) = ZERO 60 !$omp section 61 A_M(:,south,M) = ZERO 62 !$omp section 63 A_M(:,west,M) = ZERO 64 !$omp section 65 A_M(:,0,M) = -ONE 66 !$omp section 67 A_M(:,east,M) = ZERO 68 !$omp section 69 A_M(:,north,M) = ZERO 70 !$omp section 71 A_M(:,top,M) = ZERO 72 !$omp section 73 B_M(:,M) = ZERO 74 !$omp end sections 75 !$omp end parallel 76 ENDIF 77 RETURN 78 END SUBROUTINE INIT_AB_M 79