To start your first run, after compiling, enter the following command line in the parent directory of the
src/ directories :
This makes the simulation of a Jupiter mass planet initially on a circular orbit at r=1, in a disk extending from 0.04 to 20. The 2D grid extends from 0.35 to 2.5 (but the rings between 0.35 and 0.4 and between 2.45 and 2.5 are ghosts). The rest is covered by the 1D grid. The resolution is
The aspect ratio is H/r=0.04 and the viscosity is constant . The initial density profile is .
The length of the simulation is 200 orbits, and an output is produced every 10 orbits. The density profile evolution can be plotted, using the
gasdens.ascii.N.dat files. For instance with gnuplot :