Open Morisset opened 4 years ago
You are right! That's because I made a mistake in calling the Service_profile_disk function:
Sersic_profiledisk(A, Re, B, n, r, r_ * np.cos(theta))
which should actually be:
Sersic_profiledisk(A, Re, B, n, np.sqrt((r * np.cos(phi)) * 2 + (r_ np.sin(phi) / q) * 2), r_ np.cos(theta)) where q is the axis ratio!
Thanks and I will update it!
Best, Zeren
Hi, I do not understand why you are running a loop on phi, which does not contain any parameter varying with phi? Why not using a single theta variation, the use the axy-symmetry of the situation to obtain your disk? Christophe