ADC input for fargo3D
We need to scan each field and find out where something specific need to be done for 2D.
- velocity
- radial: same
- theta: same
- phi: set to 0 in 2D
- soundSpeed: same
- pressure: same
- viscosity: same
- energy: for 2D:
real
Fenergy( DiskDesc const& disk, int i, real phi) { // remove phi and propagate \pi/2
DiskGrid const& grid = disk.grid();
auto const& Rmed = grid.rMed();
auto const& aspectRatio = disk.aspectRatio();
assert(disk.physic().adiabaticIndex != 1.0); // if only for rounding....
real rcil = Rmed[i]*sin(phi); // rcil == rmed[i]
real energy0 = (std::pow(aspectRatio[i],2)*DensityTheo(disk, Rmed[i],phi)// 2D DensityTheo
*std::pow(rcil,-1.0)/(disk.physic().adiabaticIndex-1.0));
return energy0;
}
- density: rely on
DensityTheo
overloads. - deadDensity: does not exist in 2D, no deadDensity field in GasDisk for 2D.
- coolingTime: same