C TRJCHA    SOURCE    OF166741  24/10/03    21:15:41     12022          
      SUBROUTINE TRJCHA(IPARPO,MCHELM,MMODEL,IPART)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C     Construction d'un CHAMELEM representant les trajectoires
C     et du modèle correspondant
C     Les coordonnées des points sont dans un maillage et le temps
C     dans un segment MELVAL
C     MCHELM est supposé actif et suffisament dimensionné
C     on crée un segment MCHAML pour chaque trajectoire
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
-INC SMELEME
-INC PPARAM
-INC SMCHAML
-INC SMMODEL
      SEGMENT  IPARPO
        INTEGER NAPAR(NPOS),NUMP(NPOS)
        REAL*8 CREF(NDIM,NPOS),TPAR(NPOS)
      ENDSEGMENT
C     GENERATION DU MAILLAGE
      NBELEM=TPAR(/1)-1
      NDIM=CREF(/1)
      NBSOUS=0
      NBREF=0
      NBNN=2
      SEGINI MELEME
      ITYPEL=2
      ZBB=0.D0
      XBB=CREF(1,1)
      YBB=CREF(2,1)
      IF(NDIM.EQ.3)ZBB=CREF(3,1)
      CALL CREPO1(XBB,YBB,ZBB,IP1)
      IF(NDIM.EQ.2)THEN
           DO 10 I=1,NBELEM
              XBB=CREF(1,I+1)
              YBB=CREF(2,I+1)
              CALL CREPO1(XBB,YBB,ZBB,IP2)
              NUM(1,I)=IP1
              NUM(2,I)=IP2
              IP1=IP2
   10      CONTINUE
      ELSEIF(NDIM.EQ.3)THEN
           DO 20 I=1,NBELEM
              XBB=CREF(1,I+1)
              YBB=CREF(2,I+1)
              ZBB=CREF(3,I+1)
              CALL CREPO1(XBB,YBB,ZBB,IP2)
              NUM(1,I)=IP1
              NUM(2,I)=IP2
              IP1=IP2
   20      CONTINUE
      ENDIF

C     Sauvegarde dans le CHAMELEM
      IMACHE(IPART)=MELEME
      CONCHE(IPART)=' '
      INFCHE(IPART,6)=1
      N2=1
      SEGINI MCHAML
      ICHAML(IPART)=MCHAML
      NOMCHE(1)='TMPS'
      TYPCHE(1)='REAL*8'
      N1PTEL=2
      N1EL=NBELEM
      N2PTEL=0
      N2EL=0
      SEGINI MELVAL
      IELVAL(1)=MELVAL
      DO 30 I=1,NBELEM
            VELCHE(1,I)=TPAR(I)
            VELCHE(2,I)=TPAR(I+1)
   30 CONTINUE

C     Création du modèle correspondant
C     INFMOD et MATMOD sont pipos, CONMOD est vide, élément fini=2 (SEG2)
      MN3=1
      NFOR=1
      NMAT=1
      NPARMO=0
      nobmod=0
      SEGINI IMODEL
      KMODEL(IPART)=IMODEL
      IMAMOD=MELEME
      NEFMOD=2
      CONMOD=' '
      FORMOD(1)(1:12)='TRAJECTOIRE'
      MATMOD=' '

      SEGDES MELVAL
      SEGDES IMODEL
      SEGDES MCHAML
      SEGDES MELEME
      RETURN
      END












 
 
 
