corotating binary configuration with equal-mass
particles in both stars, you can use the
``setup2cm'' subroutine, which is called by setting the three
letter code in sph.init to ``2cm''.
The code lays down each star in an analogous manner to the equal-mass
single star routines, and the same instructions about the relaxation
timescale and shutoff time apply.
Note that the parameter ``N'' used to initialize the run refers
to the number of particles in EACH star, the total number of particles
will be 
.