C
C SUBROUTINE MEMORISANT LA TRAJECTOIRE DES PARTICULES SELECTIONNEES	
C	
	SUBROUTINE  TIMPA(NOM,NUMPA,NUMTIM,NK,ELX,ELY
     &                        ,ERX1,ERY1,ERZ,BRX,BRY,BRZ,TRAJ)
	
	INCLUDE "PARAM"
               PARAMETER (FIELD1=FIELD3+2)
	
     	COMMON/PART/X(PART),Y(PART)
     	COMMON/VITE/PX(PART),PY(PART),PZ(PART),PSA(PART)

     	COMMON/PERIOD/IR(FIELD3),IL(FIELD3),JR(FIELD2),JL(FIELD2)
      
     	DIMENSION NOM(NUMPA),TRAJ(NUMPA,NUMTIM,NUMFIC)

     	DIMENSION ELX(FIELD1,FIELD2),ELY(FIELD1,FIELD2)
     	DIMENSION ERX1(FIELD1,FIELD2),ERY1(FIELD1,FIELD2),
     &              ERZ(FIELD1,FIELD2)
     	DIMENSION BRX(FIELD1,FIELD2),BRY(FIELD1,FIELD2),BRZ(FIELD1,FIELD2)

     	DO 10 I=1,NUMPA
     	TRAJ(I,NK,1)=X(NOM(I))
     	TRAJ(I,NK,2)=Y(NOM(I))
     	TRAJ(I,NK,3)=PX(NOM(I))
     	TRAJ(I,NK,4)=PY(NOM(I))
     	TRAJ(I,NK,5)=PZ(NOM(I))
     	
     	
     	LX=X(NOM(I))+1.5
     	DX=(X(NOM(I))-LX+1.)*0.5
     	LXL=IL(LX)
     	LXR=IR(LX)
     	
     	LY=Y(NOM(I))+1.5
     	DY=(Y(NOM(I))-LY+1)*0.5
     	LYL=JL(LY)
    	LYR=JR(LY)
    	 
    	EXL=ELX(LX,LY)+DX*(ELX(LXR,LY)-ELX(LXL,LY))+DY*(ELX(LX,LYR)
     &            -ELX(LX,LYL))
     	EYL=ELY(LX,LY)+DX*(ELY(LXR,LY)-ELY(LXL,LY))+DY*(ELY(LX,LYR)
     &             -ELY(LX,LYL))
         
         
     	EXT=ERX1(LX,LY)+DX*(ERX1(LXR,LY)-ERX1(LXL,LY))+DY*(ERX1(LX,LYR)
     &          -ERX1(LX,LYL))
     	EYT=ERY1(LX,LY)+DX*(ERY1(LXR,LY)-ERY1(LXL,LY))+DY*(ERY1(LX,LYR)
     &           -ERY1(LX,LYL))
     	EZT=ERZ(LX,LY)+DX*(ERZ(LXR,LY)-ERZ(LXL,LY))+DY*(ERZ(LX,LYR)
     &           -ERZ(LX,LYL))
         
     	BXT=BRX(LX,LY)+DX*(BRX(LXR,LY)-BRX(LXL,LY))+DY*(BRX(LX,LYR)-BRX(LX
     &,LYL))
     	BYT=BRY(LX,LY)+DX*(BRY(LXR,LY)-BRY(LXL,LY))+DY*(BRY(LX,LYR)-BRY(LX
     &,LYL))
     	BZT=BRZ(LX,LY)+DX*(BRZ(LXR,LY)-BRZ(LXL,LY))+DY*(BRZ(LX,LYR)-BRZ(LX
     &,LYL))
     	
     	
     	TRAJ(I,NK,6)=EXT
     	TRAJ(I,NK,7)=EYT
     	TRAJ(I,NK,8)=EZT   	
     	TRAJ(I,NK,9)=BXT
     	TRAJ(I,NK,10)=BYT
     	TRAJ(I,NK,11)=BZT   	
     	TRAJ(I,NK,12)=EXL
     	TRAJ(I,NK,13)= EYL	
     	
   10   CONTINUE
     	
     	RETURN
     	END