Esope source files

   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 :