More comments added and csc notation and safe value

aart_func/ellipi_f.py

@@ -1,9 +1,10 @@
 from aart_func import *
 from params import * 
 
-def ellipi(n,phi,m): 
-    csc=1/np.sin(phi)**2
-    return 1/3*n*elliprj(csc-1,csc-m,csc,csc-n)+ellipf(phi,m)
+#Legendre’s Integrals as Symmetric Integrals https://dlmf.nist.gov/19.25#i
+def ellipi(n,phi,m):
+    csc2=1/(np.sin(phi)**2+1e-16)
+    return 1/3*n*elliprj(csc2-1,csc2-m,csc2,csc2-n)+ellipf(phi,m)
 
 def ellippi_lim(n,phi,m):
     #take arrays n, m, phi can either be a float or an array

params.py

@@ -7,9 +7,6 @@
 spin_case=0.94
 #Observer's inclination  
 i_case=17
-#Disk's inclination       
-i_disk=90    
-
 
 # Distance to M87 in meters
 dM87=5.214795112e23  
@@ -58,15 +55,16 @@
 #Number of snapshots in that range    
 snapshots=12
 
-# SU's parameters for the envelope 
-# Just used for the profiles computed within AART
 isco = rms(spin_case)
-
+# SU's parameters for the envelope 
+# Just used for the analytical profiles
 gammap=-3/2
 mup=1-sqrt(1-spin_case**2)
 sigmap=1/2 
-speed_p=1
+# Sarting point of the profile's cutoff
 cutoff_p=limits-1
+#Speed/slope of the cutoff
+speed_p=1
 
 #The power of the redshift factor
 gfactor=3
@@ -80,7 +78,11 @@
 #For the parallel generation of images (movies)
 nthreads=4
 
+#Observer's inclination in radians
 thetao=i_case*np.pi/180
+#Disk's inclination  
+#Current version is just implemented for equatorial models   
+i_disk=90    
 thetad=i_disk*np.pi/180
 
 Gc=6.67e-11 # G constant [m^3 kg^-1 s^-2]