subroutine trajectory(t) 

	structure /traj/
	   real*8 epoch
	   real xgse,ygse,zgse
	   real xgsm,ygsm,zgsm
	   real xgei,ygei,zgei
           real r_km,r_re,mlt,mlat,lshell
	   real ra_sat,dec_sat,ra_vsat,dec_vsat,vsat ! km/sec
           real gei_gse(3,3),gse_gsm(3,3)
	end structure
        record /traj/t

        REAL*8 TRJ,CD,tod
C**********************************************************************
        DIMENSION ID(324),TRJ(77),CD(21)
C**********************************************************************
        EQUIVALENCE (ID(130),CD(1))
        EQUIVALENCE (ID(180),TRJ(1))
	common/trecord/id
C
C*************************************************************************
C       TRAJECTORY INFORMATION
C*************************************************************************
C
c	trj(4)  longitude in gse
c	trj(5)  latitude in gse
c	trj(6)  longitude in gmc converted to Magnetic local time
c       amlt =  mod( ( (trj(6)-trj(28)+180.)/15.+24.0 ),24.d0)
c	trj(7)  latitude in gmc
c	trj(9)  radius in RE
c	trj(10) gse-x
c	trj(11) gse-y
c	trj(12) gse-z
c	trj(13) gsm-x
c	trj(14) gsm-y
c	trj(15) gsm-z
	tod =(id(5)*60d0+id(6))*60d0+id(7)/100d0 
	t.epoch = epoch0(id(3),id(4),tod)
        t.mlat = trj(7)
	t.mlt =  mod( ( (trj(6)-trj(28)+180.)/15.+24.0 ),24.d0)
	t.r_re = trj(8)
        t.r_km = trj(9)
	t.xgse = trj(10)
 	t.ygse = trj(11)
 	t.zgse = trj(12)
	t.xgsm = trj(13)
 	t.ygsm = trj(14)
 	t.zgsm = trj(15)
	t.xgei = trj(22)
 	t.ygei = trj(23)
 	t.zgei = trj(24)
	do i=1,3
           t.gse_gsm(1,i)=trj(31+i)
           t.gse_gsm(2,i)=trj(34+i)
           t.gse_gsm(3,i)=trj(37+i)
           t.gei_gse(1,i)=trj(40+i)
           t.gei_gse(2,i)=trj(43+i)
           t.gei_gse(3,i)=trj(46+i)
	enddo
	t.ra_sat = trj(50)
	t.dec_sat = trj(51)
	t.ra_vsat = trj(52)
	t.dec_vsat = trj(53)
	t.vsat = trj(54)
	t.lshell = trj(55)
c	t. = trj()
c        WRITE (*,503) TRJ(9),TRJ(4),TRJ(5)
c        WRITE (*,8)
c  503   FORMAT (23X,4HR = ,F7.3,3X,6HLON = ,F7.3,3X,6HLAT = ,F6.3)
c        WRITE (*,8)
c        WRITE (*,510)
c        WRITE (*,8)
c        WRITE (*,504) CD(10),CD(11),CD(12),CD(13),CD(14)
c  504   FORMAT (13X,5F10.3)
c        WRITE (*,8)
c        WRITE (*,505)
c  505   FORMAT (23X,'SATELLITE POSITION AND VELOCITY IN CI')
c        WRITE (*,8)
c        WRITE (*,506) TRJ(50),TRJ(51),TRJ(52),TRJ(53)
c  506   FORMAT(20X,'RAS =',F7.3,'DEC =',F7.3,'REV =',F7.3,'DEC =',F7.3)
c        WRITE (*,505)
c        WRITE (*,8)
c        WRITE (*,507)
c  507   FORMAT (35X,'GSE TO GSM')
c        WRITE (*,8)
c        WRITE (*,508) TRJ(32),TRJ(33),TRJ(34)
c        WRITE (*,508) TRJ(35),TRJ(36),TRJ(37)
c        WRITE (*,508) TRJ(38),TRJ(39),TRJ(40)
c  508   FORMAT (14X,3F15.6)
c        WRITE (*,8)
c        WRITE (*,509)
c  509   FORMAT (35X,'GEI TO GSE')
c        WRITE (*,8)
c        WRITE (*,508) TRJ(41),TRJ(42),TRJ(43)
c        WRITE (*,508) TRJ(44),TRJ(45),TRJ(46)
c        WRITE (*,508) TRJ(47),TRJ(48),TRJ(49)
c        WRITE (*,8)
c        WRITE (*,8)
c  510   FORMAT (31X,'ATTITUDE DATA')
	return
	end