Télécharger la plaquette Cast3M


fr en

Télécharger trjsto.eso

Retour à la liste

Numérotation des lignes : 

trjsto

  1. C TRJSTO    SOURCE    CHAT      05/01/13    03:51:14     5004
  2.       SUBROUTINE TRJSTO(NDIM,IEL1,NPOS,ITER,IPARPO,TDEP,XDEP2)
  3. ****************************************************************************
  4. ***  Permet de sauvegarder pdt l'avancee de la particule les
  5. *** coord reelles de la particule ainsi que le tps ecoule associé.
  6. ***
  7. *** APPELES 1 = aucun
  8. *** APPELES 2 = aucun
  9. ***
  10. *** E = 'NDIM' dimension de l'espace
  11. ***     'IEL1' n° global element contenant particule
  12. ***
  13. *** E/S = 'NPOS' taille maximale des tableaux du segment 'IPARPO'
  14. ***       'ITER' indice des tableaux de 'IPARPO' pour la sauvegarde
  15. ***       'IPARPO' segment ou sont sauvegardés pts trajectoire particule
  16. ***       'TDEP' tps reel départ particule sauvegardé
  17. ***       'XDEP2' position réelle de départ particule sauvegardé
  18. ***
  19. ***************************************************************************
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.       SEGMENT  IPARPO
  24.          INTEGER NAPAR(NPOS),NUMP(NPOS)
  25.          REAL*8 CREF(NDIM,NPOS),TPAR(NPOS)
  26.       ENDSEGMENT
  27.       DIMENSION XDEP2(3)
  28. ***   allongement du segment 'IPARPO' si nécessaire
  29.       IF ((NPOS-ITER).LE.1) THEN
  30.          NPOS=NPOS+50
  31.          SEGADJ IPARPO
  32.       ENDIF
  33. ***   on passe à l'iteration suivante
  34.       ITER=ITER+1
  35. ***   affectation du n° element et coords réelles depart
  36.       NAPAR(ITER)=IEL1
  37.       DO 10 I=1,NDIM
  38.          CREF(I,ITER)=XDEP2(I)
  39.  10   CONTINUE
  40. ***   affectation du tps réel depart
  41.       TPAR(ITER)=TDEP
  42.  
  43.       RETURN
  44.       END
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  

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