display_resid.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Module name: DISPLAY_RESID(NIT, IER)                                !
!  Author: M. Syamlal                                 Date: 8-JUL-96   !
!                                                                      !
!  Purpose: Display residuals                                          !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE display_resid(NIT)

      USE residual, only: group_resid

      IMPLICIT NONE

! iteration number
      INTEGER, INTENT(IN) :: NIT

! Print Location of Max_Resid
!      LOGICAL,PARAMETER:: Print_ijk=.FALSE.


      IF(group_resid) THEN
         CALL display_group_resid
      ELSE
         CALL display_field_resid
      ENDIF

!     IF(PRINT_IJK) 





!
!
!     Display maximum values of residuals
!     IF(PRINT_IJK) WRITE(*,'(A, G12.3, 3I6, A, G12.3, 3I6, A, G12.3)') &
!     & " Max Res/IJK: P_g: ", MAX_RESID(RESID_P, 0), &
!     & I_OF_G(IJK_RESID(RESID_P, 0)), &
!     & J_OF_G(IJK_RESID(RESID_P, 0)), &
!     & K_OF_G(IJK_RESID(RESID_P, 0)), &
!     & " P_s: ", MAX_RESID(RESID_p, 1), &
!     & I_OF_G(IJK_RESID(RESID_p, 1)), &
!     & J_OF_G(IJK_RESID(RESID_p, 1)), &
!     & K_OF_G(IJK_RESID(RESID_p, 1)), &
!     & " P_star=",  P_star(IJK_RESID(RESID_p, 1))
!
      RETURN

      contains


!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Module name: DISPLAY_FIELD_RESID(NIT, IER)                          !
!  Author: M. Syamlal                                 Date: 8-JUL-96   !
!                                                                      !
!  Purpose: Display residuals for each field variable.                 !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE display_field_resid

      use param1, only: undefined_i
      use residual, only: resid_string, resid_index, resid

      use error_manager

      IMPLICIT NONE

      INTEGER :: LL, LC, LS, LE

      IF(nit == 1) THEN
         WRITE(err_msg(1)(1:5),'("  Nit")')
         lc=1
         DO ll = 1, 8
            IF (resid_index(ll,1) /= undefined_i) THEN
               ls= 6+10*(lc-1)
               le= 5+10*(lc)
               WRITE(err_msg(1)(ls:le),'(5X,A4)') resid_string(ll)
               lc=lc+1
            ENDIF
         END DO
         IF(resid_index(8,1) == undefined_i) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),'(2X,A7)') 'Max res'
         ENDIF
         CALL flush_err_msg(header=.false., footer=.false., log=.false.)
      ENDIF


      WRITE(err_msg(1)(1:5),'(I5)') nit
      lc=1
      DO ll = 1, 8
         IF(resid_index(ll,1) /= undefined_i) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),'(2X,1PG8.1)') &
               resid(resid_index(ll,1),resid_index(ll,2))
            lc=lc+1
         ENDIF
      ENDDO
      IF(resid_index(8,1) == undefined_i) THEN
         ls= 6+10*(lc-1)
         le= 3+10*(lc)
         WRITE(err_msg(1)(ls:le),'(4X,A4)') resid_string(8)
      ENDIF
      CALL flush_err_msg(header=.false., footer=.false., log=.false.)

      RETURN

      END SUBROUTINE display_field_resid



!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Module name: DISPLAY_GROUP_RESID(NIT, IER)                          !
!  Author: M. Syamlal                                 Date: 8-JUL-96   !
!                                                                      !
!  Purpose: Display residuals grouped by equation type.                !
!                                                                      !
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
      SUBROUTINE display_group_resid

      use residual, only: resid_string
      use residual, only: resid_grp, resid_grp_string
      use residual, only: hydro_grp, theta_grp, energy_grp
      use residual, only: species_grp, scalar_grp, ke_grp
      use run, only: granular_energy, energy_eq
      use run, only: any_species_eq, k_epsilon
      use scalars, only : nscalar

      use error_manager

      IMPLICIT NONE


      INTEGER :: LC, LS, LE

      IF (nit == 1) THEN
         WRITE(err_msg(1)(1:5),'("  Nit")')
         lc=1

         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1000) resid_grp_string(hydro_grp)
         lc=lc+1

         IF(granular_energy) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),1000) resid_grp_string(theta_grp)
            lc=lc+1
         ENDIF

         IF(energy_eq) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),1000) resid_grp_string(energy_grp)
            lc=lc+1
         ENDIF

         IF(any_species_eq) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),1000) resid_grp_string(species_grp)
            lc=lc+1
         ENDIF

         IF(nscalar > 0) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),1000) resid_grp_string(scalar_grp)
            lc=lc+1
         ENDIF

         IF(k_epsilon) THEN
            ls= 6+10*(lc-1)
            le= 5+10*(lc)
            WRITE(err_msg(1)(ls:le),1000) resid_grp_string(ke_grp)
            lc=lc+1
         ENDIF

         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1000) 'Max res '

         CALL flush_err_msg(header=.false., footer=.false., log=.false.)
      ENDIF

 1000 FORMAT(3x,a7)



     WRITE(err_msg(1)(1:5),'(I5)') nit
      lc=1

      ls= 6+10*(lc-1)
      le= 5+10*(lc)
      WRITE(err_msg(1)(ls:le),1100) resid_grp(hydro_grp)
      lc=lc+1

      IF(granular_energy) THEN
         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1100) resid_grp(theta_grp)
         lc=lc+1
      ENDIF

      IF(energy_eq) THEN
         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1100) resid_grp(energy_grp)
         lc=lc+1
      ENDIF

      IF(any_species_eq) THEN
         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1100) resid_grp(species_grp)
         lc=lc+1
      ENDIF

      IF(nscalar > 0) THEN
         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1100) resid_grp(scalar_grp)
         lc=lc+1
      ENDIF

      IF(k_epsilon) THEN
         ls= 6+10*(lc-1)
         le= 5+10*(lc)
         WRITE(err_msg(1)(ls:le),1100) resid_grp(ke_grp)
         lc=lc+1
      ENDIF

      ls= 6+10*(lc-1)
      le= 3+10*(lc)
      WRITE(err_msg(1)(ls:le),'(4X,A4)') resid_string(8)

 1100 FORMAT(2x,1pg8.1)

      CALL flush_err_msg(header=.false., footer=.false., log=.false.)

      RETURN
      END SUBROUTINE display_group_resid

      END SUBROUTINE display_resid