Sources Esope | Cast3M

modpoi
C MODPOI    SOURCE    FD218221  23/02/13    21:15:11     11591          
      SUBROUTINE MODPOI (R1,R2,R3,DENS,IPOIN)
************************************************************************
*
*                             M O D P O I
*                             -----------
*
* FONCTION:
* ---------
*
*     MODIFICATION DE LA DENSITE ET DES COORDONNEES D'UN POINT
*
*
*
* MODE D'APPEL:
* -------------
*
*     CALL MODPOI (R1,R2,R3,DENS,IPOIN)
*
* PARAMETRES:   (E)=ENTREE   (S)=SORTIE
* -----------
*
*     R1      REEL DP   (E)  1-ERE COORDONNEE DU POINT.
*     R2      REEL DP   (E)  2-EME COORDONNEE DU POINT.
*     R3      REEL DP   (E)  3-EME COORDONNEE DU POINT (QUELCONQUE SI
*                            PROBLEME BIDIMENSIONNEL).
*     DENS    REEL DP   (E)  DENSITE DU POINT
*     IPOIN   ENTIER    (E)  NUMERO DU POINT
*
      IMPLICIT INTEGER(I-N)
      REAL*8 R1,R2,R3,DENS
*
* MODULES UTILISES:
* -----------------
 
-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
*
* COMMENTAIRE
* -----------
*
*     ON NE VERIFIE PAS DANS CE SOUS-PROGRAMME QUE LE POINT EXISTE
*
* AUTEUR, DATE DE CREATION:
* -------------------------
*
*     PASCAL MANIGOT     23 NOVEMBRE 1987
*
* LANGAGE:
* --------
*
*     ESOPE + FORTRAN77
*
************************************************************************
*
      IF (IDIM .LT. 2 .OR. IDIM .GT. 3) THEN
         NUMERR = 219
         CALL ERREUR (NUMERR)
         RETURN
      END IF
*
      SEGACT,MCOORD*MOD
      NBPT1=IPOIN
      XCOOR((NBPT1-1)*(IDIM+1)+1) = R1
      XCOOR((NBPT1-1)*(IDIM+1)+2) = R2
      IF (IDIM .EQ. 3) THEN
         XCOOR((NBPT1-1)*(IDIM+1)+3) = R3
      END IF
      XCOOR(NBPT1*(IDIM+1)) = DENS
*
      END