diff --git a/usr/dehnen/falcON/src/public/lib/profile.cc b/usr/dehnen/falcON/src/public/lib/profile.cc index 46f29ad13..4f616c5f3 100644 --- a/usr/dehnen/falcON/src/public/lib/profile.cc +++ b/usr/dehnen/falcON/src/public/lib/profile.cc @@ -335,7 +335,7 @@ mr(0), rr(0), sd(0), vl(0), vr(0), sl(0), ba(0), ph(0), al(0) Mvr += xi * (mi*vp); add_outer_product2(Mxx,xi,mi); } - sd[i] = M / (Pi * i? square(R[ir[i+1]]) - square(R[ir[i-1]]) + sd[i] = M / Pi / (i? square(R[ir[i+1]]) - square(R[ir[i-1]]) : square(R[ir[i+1]]) ); vl[i] = Mvl/M; vr[i] = abs(Mvr)/M; @@ -361,7 +361,7 @@ mr(0), rr(0), sd(0), vl(0), vr(0), sl(0), ba(0), ph(0), al(0) M += mi; add_outer_product2(Mxx,xi,mi); } - sd[i] = M / (Pi * i? square(R[ir[i+1]]) - square(R[ir[i-1]]) + sd[i] = M / Pi / (i? square(R[ir[i+1]]) - square(R[ir[i-1]]) : square(R[ir[i+1]]) ); symmetrize2(Mxx); double ab = Mxx[0][0]-Mxx[1][1]; diff --git a/usr/dehnen/falcON/src/public/manip/projprof.cc b/usr/dehnen/falcON/src/public/manip/projprof.cc index f65eb0cd8..92878d196 100644 --- a/usr/dehnen/falcON/src/public/manip/projprof.cc +++ b/usr/dehnen/falcON/src/public/manip/projprof.cc @@ -127,7 +127,7 @@ namespace falcON { namespace Manipulate { }; ////////////////////////////////////////////////////////////////////////////// inline void projprof::print_line(bool V) const { - OUT <<"#--------------------------------------------------"; + OUT <<"#-----------------------------------------------------------------"; if(V) OUT<<"----------------------------------------------------"; OUT <<'\n'; @@ -171,6 +171,9 @@ namespace falcON { namespace Manipulate { else OUT<<'\n'; if(X0) OUT<<"# xcen = "<<(*X0)<<'\n'; if(V0) OUT<<"# vcen = "<<(*V0)<<'\n'; + OUT << "# Nb = " << PP.Nb() << '\n'; + OUT << "# N = " << PP.N() << '\n'; + OUT << "# Mtot = " << PP.Mtot() << '\n'; OUT <<"#\n" <<"# radius Sigma"; if(PP.has_vels())