Mass Ratio q1.0: Equal-mass Particles

To construct a corotating binary configuration for stars of unequal mass, but with equal-mass SPH particles in both stars, you can use the ``setup2qm'' subroutine, which is called by setting the three letter code in sph.init to ``2qm''. 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. The mass ratio is determined by the parameter ${\tt QDAR}=M_2/M_1$.

Note that the parameter ``N'' used to initialize the run refers to the number of particles in the primary star, which we typically assume to be the larger star, although no attempt is made to enforce this in the code. The number of particles in the secondary is given by ${\tt N}_2={\rm INT}({\tt QDAR}\times{\tt N})$. As always make sure the parameter ``NMAX'' in the spha.h header file is as large or larger than the total number of particles needed!