Télécharger la plaquette Cast3M


fr en

Télécharger trjcha.eso

Retour à la liste

Numérotation des lignes : 

trjcha

  1. C TRJCHA    SOURCE    OF166741  24/10/03    21:15:41     12022          
  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 PPARAM
  16. -INC SMCHAML
  17. -INC SMMODEL
  18.       SEGMENT  IPARPO
  19.         INTEGER NAPAR(NPOS),NUMP(NPOS)
  20.         REAL*8 CREF(NDIM,NPOS),TPAR(NPOS)
  21.       ENDSEGMENT
  22. C     GENERATION DU MAILLAGE
  23.       NBELEM=TPAR(/1)-1
  24.       NDIM=CREF(/1)
  25.       NBSOUS=0
  26.       NBREF=0
  27.       NBNN=2
  28.       SEGINI MELEME
  29.       ITYPEL=2
  30.       ZBB=0.D0
  31.       XBB=CREF(1,1)
  32.       YBB=CREF(2,1)
  33.       IF(NDIM.EQ.3)ZBB=CREF(3,1)
  34.       CALL CREPO1(XBB,YBB,ZBB,IP1)
  35.       IF(NDIM.EQ.2)THEN
  36.            DO 10 I=1,NBELEM
  37.               XBB=CREF(1,I+1)
  38.               YBB=CREF(2,I+1)
  39.               CALL CREPO1(XBB,YBB,ZBB,IP2)
  40.               NUM(1,I)=IP1
  41.               NUM(2,I)=IP2
  42.               IP1=IP2
  43.    10      CONTINUE
  44.       ELSEIF(NDIM.EQ.3)THEN
  45.            DO 20 I=1,NBELEM
  46.               XBB=CREF(1,I+1)
  47.               YBB=CREF(2,I+1)
  48.               ZBB=CREF(3,I+1)
  49.               CALL CREPO1(XBB,YBB,ZBB,IP2)
  50.               NUM(1,I)=IP1
  51.               NUM(2,I)=IP2
  52.               IP1=IP2
  53.    20      CONTINUE
  54.       ENDIF
  55.  
  56. C     Sauvegarde dans le CHAMELEM
  57.       IMACHE(IPART)=MELEME
  58.       CONCHE(IPART)=' '
  59.       INFCHE(IPART,6)=1
  60.       N2=1
  61.       SEGINI MCHAML
  62.       ICHAML(IPART)=MCHAML
  63.       NOMCHE(1)='TMPS'
  64.       TYPCHE(1)='REAL*8'
  65.       N1PTEL=2
  66.       N1EL=NBELEM
  67.       N2PTEL=0
  68.       N2EL=0
  69.       SEGINI MELVAL
  70.       IELVAL(1)=MELVAL
  71.       DO 30 I=1,NBELEM
  72.             VELCHE(1,I)=TPAR(I)
  73.             VELCHE(2,I)=TPAR(I+1)
  74.    30 CONTINUE
  75.  
  76. C     Création du modèle correspondant
  77. C     INFMOD et MATMOD sont pipos, CONMOD est vide, élément fini=2 (SEG2)
  78.       MN3=1
  79.       NFOR=1
  80.       NMAT=1
  81.       NPARMO=0
  82.       nobmod=0
  83.       SEGINI IMODEL
  84.       KMODEL(IPART)=IMODEL
  85.       IMAMOD=MELEME
  86.       NEFMOD=2
  87.       CONMOD=' '
  88.       FORMOD(1)(1:12)='TRAJECTOIRE'
  89.       MATMOD=' '
  90.  
  91.       SEGDES MELVAL
  92.       SEGDES IMODEL
  93.       SEGDES MCHAML
  94.       SEGDES MELEME
  95.       RETURN
  96.       END
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales