2D viscosity implementation for fargo3D
From @elena in #187 (closed)
I have prepared the routines for the viscosity, namely using the same structure of viscosity.cpp of D3 I have written the
ViscosityTensor::ViscosityTensor(GasDisk::Velocity& velocity,
PolarGrid const& rhoGrid,
PolarGrid const& viscosityGrid)
for the 2D case,
it will work if in diskGrid.cpp we define phi = pi/2 and phiMed=pi/2 if dim=2
now I would like to insert it in the 3D code but I do not know how:
will there be ViscosityTensor::ViscosityTensor3D ViscosityTensor::updateVelocities3D
and ViscosityTensor::ViscosityTensor2D ViscosityTensor::updateVelocities2D
???
Whatever is the choice tell me when the structure is ready I am ready to insert the viscosity routines