Belt-Questions

Makemake2.1

Gas and Dust Density Defintions

src/Makemake2.1/Irradiation.c:Irradiation(Data *data, Grid *grid)

 rhogas  = data->Vc[RHO][kl][jl][il];
 Tgas    = data->GasTemperature[kl][jl][il];
 Trad    = IrradiationTemperature;
 #if IRRADIATION == YES
 rhodust = data->DustToGasMassRatio[kl][jl][il] * rhogas;
 Tdust   = data->DustTemperature[kl][jl][il];
 #endif

• Inertia of dust is neglected?