program subtour c real d1,d2,d3,d4,d5,d6,d7,d8,d9,d10,d11,day,x,y,z real min,x,y,z,r,theta,phi integer i open(10,file='T42traj',status='old') c open(10,file='enc61a.dat',status='old') open(20,file='oT42traj',status='new') do i=1,39 read(10,*) enddo c1 read(10,*,end=99)d1,d2,d3,d4,d5,d6,d7,d8,x,y,z,d9,d10,d11 c day = d2 + d5/24 + d6/(60*24) + d7/(3600*24) + d8/(1e6*3600*24) 1 read(10,*,end=99)min,x,y,z,r,theta,phi c write(*,*)min,x,y,z,r,theta write(20,*)x,y,z,min goto 1 99 continue end