| 1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465 |
- c*************************************************************************
- c DRIFT.F
- c*************************************************************************
- c This subroutine loops thorugh the particles and calls the danby routine
- c
- c Input:
- c nbod ==> number of massive bodies (int scalar)
- c mass ==> mass of bodies (real array)
- c xj,yj,zj ==> initial position in jacobi coord
- c (real arrays)
- c vxj,vyj,vzj ==> initial position in jacobi coord
- c (real arrays)
- c dt ==> time step
- c Output:
- c xj,yj,zj ==> final position in jacobi coord
- c (real arrays)
- c vxj,vyj,vzj ==> final position in jacobi coord
- c (real arrays)
- c
- c Authors: Hal Levison
- c Date: 2/12/93
- c Last revision: 9/5/94
- subroutine drift(nbod,mass,xj,yj,zj,vxj,vyj,vzj,dt)
- include '../../swift.inc'
- c... Inputs Only:
- integer nbod
- real*8 mass(nbod),dt
- c... Inputs and Outputs:
- real*8 xj(nbod),yj(nbod),zj(nbod)
- real*8 vxj(nbod),vyj(nbod),vzj(nbod)
- c... Internals:
- real*8 etajm1,etaj,mu
- integer j,iflg
- c----
- c... Executable code
- c Take a drift forward dth
- etajm1 = mass(1)
- do j = 2,nbod
- etaj = etajm1 + mass(j)
- mu = mass(1)*etaj/etajm1
- call drift_one(mu,xj(j),yj(j),zj(j),
- & vxj(j),vyj(j),vzj(j),dt,iflg)
- if(iflg.ne.0) then
- write(*,*) ' Planet ',j,' is lost !!!!!!!!!'
- write(*,*) mu,dt
- write(*,*) xj(j),yj(j),zj(j)
- write(*,*) vxj(j),vyj(j),vzj(j)
- write(*,*) ' STOPPING '
- call util_exit(1)
- endif
- etajm1 = etaj
- enddo
- return
- end
- c--------------------------------------------------------------------------
|
