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    CB215821  24/04/12    21:17:23     11897          
  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.       N2=1
  60.       SEGINI MCHAML
  61.       ICHAML(IPART)=MCHAML
  62.       NOMCHE(1)='TMPS'
  63.       TYPCHE(1)='REAL*8'
  64.       N1PTEL=2
  65.       N1EL=NBELEM
  66.       N2PTEL=0
  67.       N2EL=0
  68.       SEGINI MELVAL
  69.       IELVAL(1)=MELVAL
  70.       DO 30 I=1,NBELEM
  71.             VELCHE(1,I)=TPAR(I)
  72.             VELCHE(2,I)=TPAR(I+1)
  73.    30 CONTINUE
  74.  
  75. C     Création du modèle correspondant
  76. C     INFMOD et MATMOD sont pipos, CONMOD est vide, élément fini=2 (SEG2)
  77.       MN3=1
  78.       NFOR=1
  79.       NMAT=1
  80.       NPARMO=0
  81.       nobmod=0
  82.       SEGINI IMODEL
  83.       KMODEL(IPART)=IMODEL
  84.       IMAMOD=MELEME
  85.       NEFMOD=2
  86.       CONMOD=' '
  87.       FORMOD(1)(1:12)='TRAJECTOIRE'
  88.       MATMOD=' '
  89.  
  90.       SEGDES MELVAL
  91.       SEGDES IMODEL
  92.       SEGDES MCHAML
  93.       SEGDES MELEME
  94.       RETURN
  95.       END
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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