copy_a.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: COPY_A_g(A_VEL, A_m, IER)                              C
!  Purpose: Copy A_VEL_g to A_m                                        C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-DEC-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 copy_a_g(A_VEL, A_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 geometry
      USE indices
      USE compar
      USE functions
      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
!-----------------------------------------------
!
!                      Stored coefficients
      DOUBLE PRECISION A_VEL(dimension_3, -3:3)
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Indices
      INTEGER          IJK
!-----------------------------------------------

      ijk = 1
      IF (ijkmax2 > 0) THEN
         a_m(:ijkmax2,west,0) = a_vel(:ijkmax2,west)
         a_m(:ijkmax2,east,0) = a_vel(:ijkmax2,east)
         a_m(:ijkmax2,south,0) = a_vel(:ijkmax2,south)
         a_m(:ijkmax2,north,0) = a_vel(:ijkmax2,north)
         a_m(:ijkmax2,bottom,0) = a_vel(:ijkmax2,bottom)
         a_m(:ijkmax2,top,0) = a_vel(:ijkmax2,top)
         ijk = ijkmax2 + 1
      ENDIF
      RETURN
      END SUBROUTINE copy_a_g
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: COPY_A_s(A_VEL, A_m, M, IER)                           C
!  Purpose: Copy A_VEL_s to A_m                                        C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-DEC-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 copy_a_s(A_VEL, A_M, 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 geometry
      USE indices
      USE physprop
      USE compar
      USE functions
      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
!-----------------------------------------------
!
!                      Stored coefficients
      DOUBLE PRECISION A_VEL(dimension_3, -3:3, dimension_m)
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Indices
      INTEGER          IJK, M
!-----------------------------------------------

      ijk = 1
      IF (ijkmax2 > 0) THEN
         a_m(:ijkmax2,west,m) = a_vel(:ijkmax2,west,m)
         a_m(:ijkmax2,east,m) = a_vel(:ijkmax2,east,m)
         a_m(:ijkmax2,south,m) = a_vel(:ijkmax2,south,m)
         a_m(:ijkmax2,north,m) = a_vel(:ijkmax2,north,m)
         a_m(:ijkmax2,bottom,m) = a_vel(:ijkmax2,bottom,m)
         a_m(:ijkmax2,top,m) = a_vel(:ijkmax2,top,m)
         ijk = ijkmax2 + 1
      ENDIF
      RETURN
      END SUBROUTINE copy_a_s

!// Comments on the modifications for DMP version implementation
!// 001 Include header file and common declarations for parallelization