February 28, 2008

mercury6.2jh is John Chamber's mercury6_2.for code, but tweaked slightly
to suit my purposes. Those tweaks include additional forces to
drive orbit decay and/or planet-migration. Go to http://www.arm.ac.uk/~jec/
to download John Chamber's original Mercury6.2 code, and consult
his manual for instructions on how to use his code. And to take advantage
of my tweaks to his code, have a look at the first few lines of 
mercury6.2jh.for.

Also included in this package are a number of IDL codes that I use
to read and inspect the output that is generated by mercury6.2jh. 
Those codes are:

g77_compile		this file contains the unix command that I use to 
			compile mercury6.2jh with the g77 compiler on my linux desktop PC

g95_compile             another compilation using the g95 compiler 

demo.pro		this idl batch file illustrates the use of the following 
			IDL subroutines. First run mercury6.2jh using the included
			*.in input files, then run demo.pro in IDL, via IDL> @demo 

read_mercury6.pro	read the output of mercury6.2jh, as well as Chamber's
			original mercury6_2 code

xv2el.pro		convert heliocentric cartesian coordinates to 
			orbit elements

bary2_coords.pro	convert heliocentric coordinates to barycentric 
			coordinates

energy2.pro		calculate the system's total energy

angular_momentum.pro	calculate the system's total angular momentum

jacobi2.pro		calculate the particles' Jacobi integrals


Questions, comments, or bug reports for any of my tweaks to Mercury can
be sent to me at the address below.

Joe Hahn
jhahn@spacescience.org