RODEO: a ROe solver for Disc-Embedded Objects
As is well-known, 99 % of the visible matter in the Universe consists of gas, which makes gas dynamics (or hydrodynamics) a very important subject in astrophysics. Numerical methods become ever more important in this field. RODEO is a C++ implementation of the Roe solver, specifically aimed at gas disc simulations (Paardekooper & Mellema 2006). The version provided here does 2D (r,phi) hydrodynamics (isothermal as well as non-isothermal), uses the FARGO algorithm (Masset 2000) to allow for large time steps, and is parallelised using MPI. It is very much a 'vanilla' version, but can serve as a starting point for doing disc simulations. The aim was to keep this version short and simple.
Quick start
Download rodeo.tar.gz. Unpack using
gzip -d rodeo.tar.gz and tar -xf
rodeo.tar. This will create a directory rodeo. Go
to the source directory rodeo/src and enter
make. This compiles the code using g++. Once
compiling has finished, go to rodeo/out, where a default
input file Disc.in can be found. Enter
../bin/rodeo and the code should start to run the problem
described in de Val-Borro et
al. (2006): an inviscid isothermal simulation of a Jupiter mass
planet. For a complete reference see rodeo.pdf.Comments or questions? Feel free to contact me.
Disclaimer
RODEO is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.