init_ab_m.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: Init_Ab_m(A_m, b_m, IJKMAX2, M, IER)                   C                     C
!  Purpose:Initialiize the sparse matrix coefficients and the          C
!           source vector.                                             C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 16-MAY-96  C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      SUBROUTINE init_ab_m(A_M, B_M, IJKMAX2A, M)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE parallel
      USE compar
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Phase index
      INTEGER          M
!
!                      Maximum dimension
      INTEGER          IJKMAX2A
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Source vector
      DOUBLE PRECISION b_m(dimension_3, 0:dimension_m)
!
!-----------------------------------------------
!
!      IJK = 1
      IF (ijkmax2a > 0) THEN
!$omp parallel
!$omp sections
         a_m(:,bottom,m) = zero
!$omp section
         a_m(:,south,m) = zero
!$omp section
         a_m(:,west,m) = zero
!$omp section
         a_m(:,0,m) = -one
!$omp section
         a_m(:,east,m) = zero
!$omp section
         a_m(:,north,m) = zero
!$omp section
         a_m(:,top,m) = zero
!$omp section
         b_m(:,m) = zero
!$omp end sections
!$omp end parallel
      ENDIF
      RETURN
    END SUBROUTINE init_ab_m