Add a taper on (gamma - gasAcceleration()) in disk.cpp
In Disk with a giant planet + no viscosity a vortex is generated at the gap's outer edge.
This perturbs the total torque (it oscillates according to vortex motion)
When releasing the planet we observe a rapid increase of the planet's eccentricity
We would like to introduce the disk action smoothly, with a taper in disk.cpp, namely:
Triplet vel = planet.velocity() + dt * (gamma - gasAcceleration());
should be replaced by:
diskTaper = (physicalTime-disk.taper.start*2*PI)/((disk.taper.final-disk.taper.start)*2*PI);
diskTaper = diskTaper > 1 ? 1 : std::pow(std::sin(diskTaper*PI/2),2);
Triplet vel = planet.velocity() + dt * diskTaper*(gamma - gasAcceleration());
@alainm I am not sure that disk.taper is the good way to do,
also this is may be a test that we do not wont to keep for long,
so I don't want to spend a lot of time in deciding where to put such parameters, having problem of compatibility with Hdf5 ... briefly I would like not to create all the troubles I did with the parameter "width" for the cavity
is there an easy way to add parameters without changing to much the code behaviour ?? especially for parameters that may be really temporary, like in this case, if the test is bad, we will look for another solution , may adding another piece of code with other parameters elsewhere ...