C CAPA5     SOURCE    CB215821  17/01/16    21:15:07     9279           

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=                3 en THERMIQUE RHO CP EPAI                           =
C=                2 en DIFFUSION CDIF EPAI                             =
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(NMATT)

C ===
C  1 - Recuperation des caracteristiques materielles
C      Ces composantes sont obligatoires donc IVAL(i) n'est pas nul !
C ===
      MPTVAL = IVAMAT
      DO i = 1, NMATT
        VACOMP(i)=XZero
        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)

C     CAS THERMIQUE on fait RHO.CP
      VACOMP(1) = VALRHO * VACOMP(2)
      RETURN
      END
      

C ===
C  2 - Ajout a la capacite calorifique du terme de chaleur latente lie
C      a un changement de phase
C      CB215821 - 2017 - Fait dans PROPHA et plus ici !!!
C ===
C      IF (IVAPHA.NE.0) THEN
C        MPTVAL = IVAPHA
Cc* NPHAT = 1 (normalement !)
C        MELVAL = IVAL(1)
C        IGMN = MIN(iGau,VELCHE(/1))
C        IEMN = MIN(iElt,VELCHE(/2))
C        VACOMP(1) = VACOMP(1) + VALRHO * VELCHE(iGau,iElt)
C      ENDIF
C
C      RETURN
C      END