C FIMAX2    SOURCE    OF166741  24/12/13    21:15:44     12097          
      SUBROUTINE FIMAX2(ICEN,IX2D,IY2D,IS2D,IVOL,IPN,IGPN,ILGPN,IRES)
C************************************************************************
C
C PROJET            :  CASTEM 2000
C
C NOM               :  FIMAX2
C
C DESCRIPTION       :  Voir FIMAXI
C
C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
C
C AUTEUR            :  A. BECCANTINI, LTMF
C
C************************************************************************
C
C
C APPELES (Outils)  :  AUCUN
C
C APPELES (Calcul)  :  AUCUN
C
C
C************************************************************************
C
C ENTREES
C
C
C     1) Pointeurs de MELEMEs et de CHPOINTs "CENTRE" de la table DOMAINE
C
C     ICEN      : MELEME de 'POI1' SPG des CENTRES
C
C     IX2D,IY2D : coordonnés des centres 2D
C
C     IS2D      : surface 2D
C     IVOL      : volume
C
C     2) CHPOINTs "CENTRE"
C
C     IPN,IGPN,ILGPN : pression, gradient, limiteur
C
C
C SORTIES
C
C     IRES    : contribution au residu
C
C
C************************************************************************
C
C HISTORIQUE (Anomalies et modifications éventuelles)
C
C HISTORIQUE :  Créée le 25.2.2004
C
C
C
C************************************************************************
C
      IMPLICIT INTEGER(I-N)
C
      INTEGER ICEN,IX2D,IY2D,IS2D,IVOL,IPN,IGPN,ILGPN,IRES
     &     , NBEL, IELEM, NP, INDP, IGEOM
      REAL*8 DX, DY, GX, GY, GEOM, AL, P
      CHARACTER*8 TYPE
C
C**** Les Includes
C
-INC SMCOORD

-INC PPARAM
-INC CCOPTIO
-INC SMCHPOI
      POINTEUR MP.MPOVAL,     MGRP.MPOVAL,
     &         MALP.MPOVAL,   MPXC2D.MPOVAL,
     &         MPYC2D.MPOVAL, MPSU2D.MPOVAL,
     &         MPVOLU.MPOVAL, MPRES.MPOVAL
-INC SMELEME
C
      CALL LICHT(IPN,MP,TYPE,IGEOM)
      CALL LICHT(IGPN,MGRP,TYPE,IGEOM)
      CALL LICHT(ILGPN,MALP,TYPE,IGEOM)
      CALL LICHT(IX2D,MPXC2D,TYPE,IGEOM)
      CALL LICHT(IY2D,MPYC2D,TYPE,IGEOM)
      CALL LICHT(IS2D,MPSU2D,TYPE,IGEOM)
      CALL LICHT(IVOL,MPVOLU,TYPE,IGEOM)
      CALL LICHT(IRES,MPRES,TYPE,IGEOM)
C
C     SEGACT   MP
C     SEGACT   MGRP
C     SEGACT   MALP
C     SEGACT   MPXC2D
C     SEGACT   MPYC2D
C     SEGACT   MPSU2D
C     SEGACT   MPVOLU
C     SEGACT   MPRES
C
      MELEME=ICEN
      SEGACT MELEME
      NBEL=MELEME.NUM(/2)
      DO IELEM=1,NBEL,1
         NP=MELEME.NUM(1,IELEM)
         INDP = (NP-1)*3
         DX=MPXC2D.VPOCHA(IELEM,1) - XCOOR(INDP+1)
         DY=MPYC2D.VPOCHA(IELEM,1) - XCOOR(INDP+2)
         GX=MGRP.VPOCHA(IELEM,1)
         GY=MGRP.VPOCHA(IELEM,2)
         GEOM=MPSU2D.VPOCHA(IELEM,1)/MPVOLU.VPOCHA(IELEM,1)
         AL=MALP.VPOCHA(IELEM,1)
         P=MP.VPOCHA(IELEM,1)
         MPRES.VPOCHA(IELEM,1)=(P+(AL*((GX*DX)+(GY*DY))))*GEOM
      ENDDO
C
      SEGDES  MELEME
      SEGDES  MP
      SEGDES  MGRP
      SEGDES  MALP
      SEGDES  MPXC2D
      SEGDES  MPYC2D
      SEGDES  MPSU2D
      SEGDES  MPVOLU
      SEGDES  MPRES
C
 9999 CONTINUE
C
      RETURN
      END