under_relax.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: UNDER_RELAX_S(Var, A_m, B_m, M, UR, IER)               C
!  Purpose: Underrelax scalar equations                                C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-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 under_relax_s(VAR, A_M, B_M, M, UR)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE indices
      USE compar
      USE sendrecv
      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
!-----------------------------------------------
!
!                      phase index
      INTEGER          M
!
!                      index
      INTEGER          IJK
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Vector b_m
      DOUBLE PRECISION B_m(dimension_3, 0:dimension_m)
!
!                      Variable
      DOUBLE PRECISION Var(dimension_3)
!
!                      Under-relaxation factor
      DOUBLE PRECISION UR
!
!                      Functions of under-relaxation factor
      DOUBLE PRECISION f1, f2
!
!                      center coefficient
      DOUBLE PRECISION ap

      f1 = one/ur
      f2 = f1 - one
!
!!!$omp    parallel do private(IJK,AP)
      DO ijk = ijkstart3, ijkend3
         IF (fluid_at(ijk)) THEN
            ap = a_m(ijk,0,m)
            IF (ap /= (-one)) THEN
               a_m(ijk,0,m) = ap*f1
               b_m(ijk,m) = b_m(ijk,m) + ap*var(ijk)*f2
            ENDIF
         ENDIF
      END DO
      RETURN
      END SUBROUTINE under_relax_s
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: UNDER_RELAX_U(Var, A_m, B_m, M, UR, IER)               C
!  Purpose: Underrelax u-momentum equations                            C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-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 under_relax_u(VAR, A_M, B_M, M, UR)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE indices
      USE compar
      USE sendrecv
      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
!-----------------------------------------------
!
!                      phase index
      INTEGER          M
!
!                      index
      INTEGER          IJK
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Vector b_m
      DOUBLE PRECISION B_m(dimension_3, 0:dimension_m)
!
!                      Variable
      DOUBLE PRECISION Var(dimension_3)
!
!                      Under-relaxation factor
      DOUBLE PRECISION UR
!
!                      Functions of under-relaxation factor
      DOUBLE PRECISION f1, f2
!
!                      center coefficient
      DOUBLE PRECISION ap

      f1 = one/ur
      f2 = f1 - one
!
!!!$omp    parallel do private(IJK,AP)
      DO ijk = ijkstart3, ijkend3
         IF (flow_at_e(ijk)) THEN
            ap = a_m(ijk,0,m)
            IF (ap /= (-one)) THEN
               a_m(ijk,0,m) = ap*f1
               b_m(ijk,m) = b_m(ijk,m) + ap*var(ijk)*f2
            ENDIF
         ENDIF
      END DO
      RETURN
      END SUBROUTINE under_relax_u
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: UNDER_RELAX_V(Var, A_m, B_m, M, UR, IER)               C
!  Purpose: Underrelax v-momentum equations                            C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-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 under_relax_v(VAR, A_M, B_M, M, UR)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE indices
      USE compar
      USE sendrecv
      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
!-----------------------------------------------

!                      phase index
      INTEGER          M
!
!                      index
      INTEGER          IJK
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Vector b_m
      DOUBLE PRECISION B_m(dimension_3, 0:dimension_m)
!
!                      Variable
      DOUBLE PRECISION Var(dimension_3)
!
!                      Under-relaxation factor
      DOUBLE PRECISION UR
!
!                      Functions of under-relaxation factor
      DOUBLE PRECISION f1, f2
!
!                      center coefficient
      DOUBLE PRECISION ap

      f1 = one/ur
      f2 = f1 - one
!
!!!$omp    parallel do private(IJK,AP)
      DO ijk = ijkstart3, ijkend3
         IF (flow_at_n(ijk)) THEN
            ap = a_m(ijk,0,m)
            IF (ap /= (-one)) THEN
               a_m(ijk,0,m) = ap*f1
               b_m(ijk,m) = b_m(ijk,m) + ap*var(ijk)*f2
            ENDIF
         ENDIF
      END DO
      RETURN
      END SUBROUTINE under_relax_v
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: UNDER_RELAX_W(Var, A_m, B_m, M, UR, IER)               C
!  Purpose: Underrelax w-momentum equations                            C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 24-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 under_relax_w(VAR, A_M, B_M, M, UR)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE indices
      USE compar
      USE sendrecv
      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
!-----------------------------------------------

!                      phase index
      INTEGER          M
!
!                      index
      INTEGER          IJK
!
!                      Septadiagonal matrix A_m
      DOUBLE PRECISION A_m(dimension_3, -3:3, 0:dimension_m)
!
!                      Vector b_m
      DOUBLE PRECISION B_m(dimension_3, 0:dimension_m)
!
!                      Variable
      DOUBLE PRECISION Var(dimension_3)
!
!                      Under-relaxation factor
      DOUBLE PRECISION UR
!
!                      Functions of under-relaxation factor
      DOUBLE PRECISION f1, f2
!
!                      center coefficient
      DOUBLE PRECISION ap

      f1 = one/ur
      f2 = f1 - one
!
!!!$omp    parallel do private(IJK,AP)
      DO ijk = ijkstart3, ijkend3
         IF (flow_at_t(ijk)) THEN
            ap = a_m(ijk,0,m)
            IF (ap /= (-one)) THEN
               a_m(ijk,0,m) = ap*f1
               b_m(ijk,m) = b_m(ijk,m) + ap*var(ijk)*f2
            ENDIF
         ENDIF
      END DO
      RETURN
      END SUBROUTINE under_relax_w

!// Comments on the modifications for DMP version implementation
!// 001 Include header file and common declarations for parallelization
!// 350 Changed do loop limits: 1,ijkmax2-> ijkstart3, ijkend3