C
C SUBROUTINE MEMORISANT LA TRAJECTOIRE DES PARTICULES SELECTIONNEES	
C	
	SUBROUTINE  TIMPA(NOM,NUMPA,NUMTIM,NK,ELX,
     &                        ,ERY,ERZ,BRY,BRZ,TRAJ)
	
	INCLUDE "PARAM"
               PARAMETER (FIELD1=FIELD3+2)
	

     	COMMON/PARTI/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)
     	LXL2=LXL-1
     	LXR2=LXR+1
    	 
    	EXL=ELX(LX)+DX*(ELX(LXR)-ELX(LXL))
        EYT=ERY(LX)+DX*(ERY(LXR)-ERY(LXL))
        EZT=ERZ(LX)+DX*(ERZ(LXR)-ERZ(LXL))
        BYT=BRY(LX)+DX*(BRY(LXR)-BRY(LXL))
        BZT=BRZ(LX)+DX*(BRZ(LXR)-BRZ(LXL))
     	
        TRAJ(I,NK,6)=EXL
        TRAJ(I,NK,7)=EYT
        TRAJ(I,NK,8)=EZT   	
        TRAJ(I,NK,9)=BYT
        TRAJ(I,NK,10)=BZT   	
     	
C
C   inclusion of all spacial gradients
C

C------ GRADIENT OF ELX ALONG X
         ELXX=(ELX(LXR)-ELX(LXL))
     &     +DX*(ELX(LXR2)-ELX(LX))
     &     -DX*(ELX(LX)-ELX(LXL2))
         ELXX=0.5*ELXX
C------ GRADIENT OF ETY ALONG X
         ETYX=(ERY(LXR)-ERY(LXL))
     &     +DX*(ERY(LXR2)-ERY(LX))
     &     -DX*(ERY(LX)-ERY(LXL2))
         ETYX=0.5*ETYX
C------ GRADIENT OF ETZ ALONG X        
         ETZX=(ERZ(LXR)-ERZ(LXL))
     &     +DX*(ERZ(LXR2)-ERZ(LX))
     &     -DX*(ERZ(LX)-ERZ(LXL2))
         ETZX=0.5*ETZX
C------ GRADIENT OF BTY ALONG X 
         BTYX=(BRY(LXR)-BRY(LXL))
     &     +DX*(BRY(LXR2)-BRY(LX))
     &     -DX*(BRY(LX)-BRY(LXL2))
         BTYX=0.5*BTYX
C------ GRADIENT OF BTZ ALONG X          
         BTZX=(BRZ(LXR)-BRZ(LXL))
     &     +DX*(BRZ(LXR2)-BRZ(LX))
     &     -DX*(BRZ(LX)-BRZ(LXL2))
         BTZX=0.5*BTZX 
   
        TRAJ(I,NK,14)=ELXX
        TRAJ(I,NK,15)=ETYX 	
     	TRAJ(I,NK,16)=ETZX
     	TRAJ(I,NK,17)=BTXX 	
     	TRAJ(I,NK,18)=BTYX	
     	TRAJ(I,NK,19)=BTZX  
        
     	
   10   CONTINUE
     	
     	RETURN
     	END