C CAPA5     SOURCE    FANDEUR   11/07/19    21:15:21     7042           

C=======================================================================
C=                            C A P A 5                                =
C=                            ---------                                =
C=                                                                     =
C=  Fonction :                                                         =
C=  ----------                                                         =
C=  Calcul, en un point de Gauss, du terme Rho.Cp intervenant dans la  =
C=  contribution de ce point de Gauss a la matrice de CAPACITE d'un    =
C=  element fini. Ce terme prend egalement en compte la presence d'un  =
C=  eventuel changement de phase.                                      =
C=                                                                     =
C=  Parametres :  (E)=Entree  (S)=Sortie                               =
C=  ------------                                                       =
C=   iGau    (E)  Numero du point de Gauss de l'element considere      =
C=   iElt    (E)  Numero de l'ELEMENT fini considere                   =
C=   IVAMAT  (E)  Segment MPTVAL                               (ACTIF) =
C=   NMATT   (E)  Nombre de composantes recherchees                    =
C=   IVAPHA  (E)  Segment MPTVAL pour changement de phase      (ACTIF) =
C=   NPHAT   (E)  Nombre de composantes recherchees                    =
C=   VACOMP  (S)  Valeur de Rho.Cp au point de Gauss...                =
C=                                                                     =
C=  Denis ROBERT, le 18 fevrier 1988.                                  =
C=======================================================================

      SUBROUTINE CAPA5 (iGau,iElt,IVAMAT,NMATT,IVAPHA,NPHAT, VACOMP)

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)

-INC CCREEL
-INC SMCHAML

      SEGMENT MPTVAL
        INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU)
        CHARACTER*16 TYVAL(NCOSOU)
      ENDSEGMENT

      DIMENSION VACOMP(*)

C ===
C  1 - Recuperation des caracteristiques materielles
C      Ces composantes sont obligatoires donc IVAL(i) n'est pas nul !
C ===
      VACOMP(1) = XZero
      VACOMP(2) = XZero
      VACOMP(3) = XZero
      MPTVAL = IVAMAT
      DO i = 1, NMATT
        MELVAL = IVAL(i)
c*      IF (IVAL(i).NE.0) THEN
          IGMN = MIN(iGau,VELCHE(/1))
          IEMN = MIN(iElt,VELCHE(/2))
          VACOMP(i) = VELCHE(IGMN,IEMN)
c*      ENDIF
      ENDDO
      VALRHO = VACOMP(1)
      VACOMP(1) = VALRHO * VACOMP(2)
C ===
C  2 - Ajout a la capacite calorifique du terme de chaleur latente lie
C      a un changement de phase
C ===
      IF (IVAPHA.NE.0) THEN
        MPTVAL = IVAPHA
c* NPHAT = 1 (normalement !)
        MELVAL = IVAL(1)
        IGMN = MIN(iGau,VELCHE(/1))
        IEMN = MIN(iElt,VELCHE(/2))
        VACOMP(1) = VACOMP(1) + VALRHO * VELCHE(iGau,iElt)
      ENDIF

      RETURN
      END