Esope source files

   1 : C TRJCHA    SOURCE    CHAT      11/03/16    21:32:43     6902           
   2 :       SUBROUTINE TRJCHA(IPARPO,MCHELM,MMODEL,IPART)
   3 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   4 : C     Construction d'un CHAMELEM representant les trajectoires
   5 : C     et du modèle correspondant
   6 : C     Les coordonnées des points sont dans un maillage et le temps
   7 : C     dans un segment MELVAL
   8 : C     MCHELM est supposé actif et suffisament dimensionné
   9 : C     on crée un segment MCHAML pour chaque trajectoire
  10 : C
  11 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8 (A-H,O-Z)
  14 : -INC SMELEME
  15 : -INC SMCHAML
  16 : -INC SMMODEL
  17 :       SEGMENT  IPARPO
  18 :         INTEGER NAPAR(NPOS),NUMP(NPOS)
  19 :         REAL*8 CREF(NDIM,NPOS),TPAR(NPOS)
  20 :       ENDSEGMENT
  21 : C     GENERATION DU MAILLAGE
  22 :       NBELEM=TPAR(/1)-1
  23 :       NDIM=CREF(/1)
  24 :       NBSOUS=0
  25 :       NBREF=0
  26 :       NBNN=2
  27 :       SEGINI MELEME
  28 :       ITYPEL=2
  29 :       ZBB=0.D0
  30 :       XBB=CREF(1,1)
  31 :       YBB=CREF(2,1)
  32 :       IF(NDIM.EQ.3)ZBB=CREF(3,1)
  32 :       CALL CREPO1(XBB,YBB,ZBB,IP1)
  33 :       IF(NDIM.EQ.2)THEN
  34 :            DO 10 I=1,NBELEM
  35 :               XBB=CREF(1,I+1)
  36 :               YBB=CREF(2,I+1)
  36 :               CALL CREPO1(XBB,YBB,ZBB,IP2)
  37 :               NUM(1,I)=IP1
  38 :               NUM(2,I)=IP2
  39 :               IP1=IP2
  40 :    10      CONTINUE
  41 :       ELSEIF(NDIM.EQ.3)THEN
  42 :            DO 20 I=1,NBELEM
  43 :               XBB=CREF(1,I+1)
  44 :               YBB=CREF(2,I+1)
  45 :               ZBB=CREF(3,I+1)
  45 :               CALL CREPO1(XBB,YBB,ZBB,IP2)
  46 :               NUM(1,I)=IP1
  47 :               NUM(2,I)=IP2
  48 :               IP1=IP2
  49 :    20      CONTINUE
  50 :       ENDIF
  51 : 
  52 : C     Sauvegarde dans le CHAMELEM
  53 :       IMACHE(IPART)=MELEME
  54 :       CONCHE(IPART)=' '
  55 :       N2=1
  56 :       SEGINI MCHAML
  57 :       ICHAML(IPART)=MCHAML
  58 :       NOMCHE(1)='TMPS'
  59 :       TYPCHE(1)='REAL*8'
  60 :       N1PTEL=2
  61 :       N1EL=NBELEM
  62 :       N2PTEL=0
  63 :       N2EL=0
  64 :       SEGINI MELVAL
  65 :       IELVAL(1)=MELVAL
  66 :       DO 30 I=1,NBELEM
  67 :             VELCHE(1,I)=TPAR(I)
  68 :             VELCHE(2,I)=TPAR(I+1)
  69 :    30 CONTINUE
  70 : 
  71 : C     Création du modèle correspondant
  72 : C     INFMOD et MATMOD sont pipos, CONMOD est vide, élément fini=2 (SEG2)
  73 :       MN3=1
  74 :       NFOR=1
  75 :       NMAT=1
  76 :       NPARMO=0
  77 :       nobmod=0
  78 :       SEGINI IMODEL
  79 :       KMODEL(IPART)=IMODEL
  80 :       IMAMOD=MELEME
  81 :       NEFMOD=2
  82 :       CONMOD=' '
  83 :       FORMOD(1)(1:12)='TRAJECTOIRE'
  84 :       MATMOD=' '
  85 : 
  86 :       SEGDES MELVAL
  87 :       SEGDES IMODEL
  88 :       SEGDES MCHAML
  89 :       SEGDES MELEME
  90 :       RETURN
  91 :       END
  92 : 
  93 : 
  94 : 
  95 : 
  96 : 
  97 : 
  98 : 
  99 : 
 100 :  
 101 :  
 102 :  
 103 :