Télécharger trjsto.eso

Retour à la liste

Numérotation des lignes :

  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