This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: Kirkwood gaps code • prog5.f as a starting point is in your home directory • Includes variable declarations, format statements… a start IMPLICIT NONE INTEGER n,nsteps PARAMETER(NSTEPS=10000) DOUBLE PRECISION vxj,vyj,vxa,vya ! velocities of jupiter, asteroid DOUBLE PRECISION vxji,vyji,vxai,vyai ! initial velocities of jupiter, asteroid DOUBLE PRECISION xj,yj,xa,ya ! positions of jupiter, asteroid DOUBLE PRECISION xja,yja ! x,y separation between asteroid and jupiter DOUBLE PRECISION xji,yji,xai,yai ! initial positions of jupiter, asteroid DOUBLE PRECISION alpha,alpha_j,alpha_a,dt,time DOUBLE PRECISION rj,ra,rja ! separation between jupiter-sun, asteroid-sun, and DOUBLE PRECISION pi DOUBLE PRECISION ms,mj,ma,ekin,epot,etot,amom Kirkwood gaps code, more declarations PARAMETER(dt=0.05d0) ! time step in units of years PARAMETER(ms=1.0d0,mj=9.55d-4,ma=4.39d-10) ! asteroid mass is that o PARAMETER(xji=0.0d0,yji=5.200d0) ! initial position of jupiter in AU PARAMETER(vxji=2.75535903d0,vyji=0.0d0) ! initial velocity of jupiter in AUPARAMETER(vxji=2....
View Full Document
This note was uploaded on 08/08/2011 for the course PHZ 5156 taught by Professor Johnson,m during the Fall '08 term at University of Central Florida.
- Fall '08