parse_resid_string.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: Parse_RESID_string(IER)                                C
!  Purpose: Initialize residuals                                       C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date:12-MAY-97   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE parse_resid_string()
!...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 physprop
      USE residual
      USE funits
      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
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   V a r i a b l e s
!-----------------------------------------------

!
!                      local index
      INTEGER          L, L1
!
!                      error message
      CHARACTER(LEN=80)     LINE
!
!                      LOGICAL
      LOGICAL          STRING_DEFINED

!
!-----------------------------------------------
!
!
!
!
!  If user did not define any residual strings use default values
!

      IF(group_resid) THEN

         resid_grp_string(hydro_grp)   = 'HYDRO'
         resid_grp_string(theta_grp)   = 'THETA'
         resid_grp_string(energy_grp)  = 'ENERGY'
         resid_grp_string(species_grp) = 'SPECIES'
         resid_grp_string(scalar_grp)  = 'SCALAR'
         resid_grp_string(ke_grp)      = 'K-EPS.'

         resid_string = undefined_c
         resid_index(8,1) = undefined_i

         RETURN
      ENDIF

      string_defined = .false.
      DO l = 1, max_resid_index
         IF (resid_string(l) /= undefined_c) string_defined = .true.
      END DO
      IF (.NOT.string_defined) THEN
         resid_string(1) = 'P0'
         resid_string(2) = 'P1'
         resid_string(3) = 'U0'
         resid_string(4) = 'V0'
         resid_string(5) = 'U1'
         resid_string(6) = 'V1'
      ENDIF
!
!C
      DO l = 1, max_resid_index
         resid_index(l,1) = undefined_i
         string_defined = .false.
         DO l1 = 1, nprefix
!
            IF (resid_string(l)(1:1) == resid_prefix(l1)) THEN
               resid_index(l,1) = l1
               string_defined = .true.
               EXIT
            ENDIF
         END DO
         IF (string_defined) THEN
!
            resid_index(l,2) = ichar(resid_string(l)(2:2)) - 48
            !print *,"RESID_INDEX(L,2) = ",RESID_INDEX(L,2)
            IF (resid_index(l,2)<0 .OR. resid_index(l,2)>mmax) THEN
               WRITE (line, '(A, A1, A, A4, A)') 'Error: Phase index ', &
                  resid_string(l)(2:2), ' in RESID_STRING ', resid_string(l), &
                  ' out of bounds.'
               CALL write_error ('PARSE_RESID_STRING', line, 1)
            ENDIF
!
            IF (resid_string(l)(1:1)=='G' .AND. resid_index(l,2)==0) THEN
!
               WRITE (line, '(A, A1, A, A4, A)') 'Error: Phase index ', &
                  resid_string(l)(2:2), ' in RESID_STRING ', resid_string(l), &
                  ' should be > 0.'
               CALL write_error ('PARSE_RESID_STRING', line, 1)
            ENDIF
            IF (resid_string(l)(1:1) == 'X') resid_index(l,1) = resid_x + &
               (ichar(resid_string(l)(3:3)) - 48) * 10 &
              + ichar(resid_string(l)(4:4)) - 48 - 1
         ENDIF
      END DO
      RETURN
      END SUBROUTINE parse_resid_string

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