FARGO features

  • Solves Navier-Stokes and continuity equations for a Keplerian disk subject to the gravity of the central object and that of embedded protoplanets.
  • Isothermal equation of state, with arbitrary temperature (or sound speed) radial profile.
 A beta version including an energy equation (adiabatic flow) can be found at fargo.in2p3.fr
  • 2D polar mesh hydrocode.
 The mesh spacing in azimuth is constant, but the radial spacing can be arbitrary.
  • uses a van Leer upwind algorithm on a staggered mesh.
  • Arbitrary number of planets.
 The planets are point-like masses. Their trajectory is determined using a fifth order Runger-Kutta integrator.
  • Includes the FARGO algorithm (Fast Advection in Rotating Gaseous Objects).
 The code can be run using either the FARGO algorithm or a standard advection technique, for comparison purposes.
  • Frame is non-inertial and centered on the primary.
  • The full viscous stress tensor of Navier-Stokes equation is implemented.
  • Parallellized both with OpenMP and MPI.
 FARGO can be run using OpenMP and MPI at the same time, on machines which have nodes of processors with shared memory. Compile time options produce either a sequential or parallel executable from the same source files.
  • Frame is either non-rotating, rotating at a fixed angular velocity, or corotating with a planet or its guiding center.
 The additional forces that arise in the rotating case are conservatively implemented.
  • A run can be stopped (or killed) and restarted.
 Some parameters can be adjusted before the restart, such as the time sampling of outputs, etc. A parallel run can be restarted as sequential, or with a different number of processors.
  • The disk is non self-gravitating.
 A beta version with the disk's self gravity implemented can be found at fargo.in2p3.fr