|
|
|
|
|
## Presentation
|
|
|
|
|
|
The code FARGOCA is a 3D version of FARGO, developed in 2011 in its isothermal version,
|
|
|
has evolved since then with the addition of the treatment of energy equation with energy provided by the viscous and compressional heating and by the energy irradiated from the star.
|
|
|
|
|
|
![fig50EMd](/uploads/5fefec1df065b698bdaaadae26014518/fig50EMd.png)
|
|
|
|
|
|
## Publications
|
|
|
|
|
|
The code '''FARGOCA''' has been used in all or part of the simulations of the papers:
|
|
|
|
|
|
1. Bitsch B., A. Morbidelli, E. Lega, K. Kretke, A. Crida, A&A, 2014, 570, 75.
|
|
|
2. Bitsch B., Morbidelli A., Lega E., Crida A., 2014, A&A, 564, A135.
|
|
|
3. Lega E., Crida A., Bitsch B., Morbidelli A., 2014, MNRAS, 440, 683.
|
|
|
4. Morbidelli A., Szulagyi J., Crida A., Lega E., Bitsch B., Tanigawa T., 2014, Icarus, 232, 266.
|
|
|
5. Lambrecths, M., Jhoansen, A., Morbidelli, A., 2014, A&A, 572.
|
|
|
6. Bitsch, B., Johansen A., Lamberechts M., Morbidelli, A., A&A 2015, 575, A28.
|
|
|
7. Lega E., Morbidelli A., Bitsch B., Crida A., Szulagyi J., 2015, MNRAS, 452, 1717-1726.
|
|
|
8. Szulagyi J., Masset, F., Lega E., Morbidelli, A., Guillot T. and Crida, A., MNRAS 2016, 460, 2853-2861
|
|
|
9. Fung, J., Masset, F., Lega E. and Velasco, D., Astron. Journal, 2017, 153
|
|
|
10. Lamberechts, M., Lega E., A&A, 2017, under revision
|
|
|
|
|
|
|
|
|
## Getting and Building the code
|
|
|
Follow [this link](Building)
|
|
|
|
|
|
## SVN content
|
|
|
This needs a rewrite.
|
|
|
|
|
|
* **fargo** original 2D version by Masset
|
|
|
* **fargo2DSymba** a version of fargo 2D where the Swift integrator replaces the 5th order Runge-Kutta for the N-body part
|
|
|
* **fargoADC** a fargo2D based version with the addition of an adiabatic equation of state, the treatment of gas-dust interaction (C. Robert master) and the modelisation of accretion disk (constant Mdot) under study
|
|
|
* **fargo3D** -> the 3D version of the code with energy, stellar heating on the whole disc with uniform and non uniform grid or on a sector with both geometries of the grid (OpenMP/MPI)
|
|
|
* **fargo3DSymba** -> the equivalent of fargo2DSymba in 3D (isothermal equation of state only)
|
|
|
* **FARGOCUDA1.2** -> fargo2D -3D (Benitez-Masset) for MPI/GPU
|
|
|
|
|
|
In order to get the code use the svn command '''chekout ''':
|
|
|
```
|
|
|
svn co https://forge.oca.eu/codes/DISK/fargo
|
|
|
svn co https://forge.oca.eu/codes/DISK/fargo3D
|
|
|
svn co https://forge.oca.eu/codes/DISK/fargoADC
|
|
|
etc ...
|
|
|
```
|
|
|
checkout of the projekct DISK (a directory DISK is created
|
|
|
and the files copied inside)
|
|
|
|
|
|
In order to compile the '''code''' 3D in ./DISK/fargo3D/src:
|
|
|
```
|
|
|
./gmake
|
|
|
```
|
|
|
|
|
|
Then ''' run''':
|
|
|
```
|
|
|
./fargo in/template.par
|
|
|
```
|
|
|
|
|
|
In order to '''compile the program for the plot''' of the gas quantities:
|
|
|
```
|
|
|
./PostProces/make
|
|
|
```
|
|
|
|
|
|
|
|
|
In order to '''make the plot of gas density and velocities ''':
|
|
|
```
|
|
|
./fargo3Dplot -h
|
|
|
./streamlines -h for streamlines
|
|
|
./average-cylindric for the analysis of flow in a cylindric frame centred on the planet
|
|
|
etc ... (see README in fargo3D/PostProces)
|
|
|
```
|
|
|
|
|
|
for the help on using fargo3Dplot. fargo3Dplot prepares different files for plot, it must be followed ./gnuexe in order to plot them with gnuplot
|
|
|
|
|
|
The output file containing the orbital parameters at each DT is orbit''i''.dat where ''i'' is the planet number.
|
|
|
The columns contain respectively:
|
|
|
```
|
|
|
the date
|
|
|
the eccentricity
|
|
|
the semi-major axis
|
|
|
the inclination
|
|
|
the planet mean anomaly
|
|
|
the argument of pericenter
|
|
|
the longitude of the node
|
|
|
```
|
|
|
|
|
|
![Cress](/uploads/1ce4238ed6d32b6464053c4bfb50fad0/Cress.jpeg)
|
|
|
|
|
|
|
|
|
Revision e [71] in [source:/fargo3D/src@71 ] corresponds to the full 3D isothermal model with gas + planets.
|
|
|
The figure shows the evolution of semi-major axis eccentricity and inclination for a 20 Earth mass planets.
|
|
|
The parameter are those in Cresswell et al. 2007, and the results are in good agreement (see Fig.16 of Cresswell et al. paper).
|
|
|
|
|
|
<!---
|
|
|
In revision [82] in [source:/fargo3DAD/src@82 ] the adiabatic part is completed.
|
|
|
In revision [88] in [source:/fargo3DAD/src@82 ] Fargo3DSymba has been tested.
|
|
|
In revision [180] in [source:/fargo3D/src@180] the radiative part has been tested
|
|
|
In revision [215] in [source:/fargo3D/src@215] the stellar radiation routine has been added
|
|
|
In revision [280] in [source:/fargo3D/src@215] Fargo-Sector tested
|
|
|
In revision [390] in [source:/fargo3D/src@390] a geometry of non-uniform grid has been added
|
|
|
--->
|
|
|
|
|
|
|
|
|
|
|
|
|