{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

lecture5 - Kirkwood gaps code prog5.f as a starting point...

This preview shows pages 1–3. Sign up to view the full content.

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

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
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 AU c PARAMETER(xai=0.0d0,yai=3.000d0) ! initial position of asteroid in AU
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}