Peano
Loading...
Searching...
No Matches
para_plot.py File Reference

Go to the source code of this file.

Namespaces

namespace  para_plot
 

Functions

 para_plot.file_len (fname)
 

Variables

 para_plot.pi = np.pi
 
float para_plot.G = 6.67408*1e-11
 
float para_plot.M = 0.15
 
int para_plot.rho_c_ini = 1
 
float para_plot.rho_c0_value = 1.8788*10**(-26)
 
float para_plot.a_i = 0.001
 
 para_plot.parser = argparse.ArgumentParser(description='Tracer_conversion')
 
 para_plot.dest
 
 para_plot.required
 
 para_plot.type
 
 para_plot.str
 
 para_plot.help
 
 para_plot.default
 
 para_plot.choices
 
 para_plot.action
 
 para_plot.False
 
 para_plot.args = parser.parse_args()
 
 para_plot.IfXunitInCode = args.IfXunitInCode
 
list para_plot.snapshots = [10,30,50]
 main code

 
int para_plot.snapshot = 0
 
list para_plot.rho = [[]];
 
list para_plot.m_1 = [[]];
 
list para_plot.E_t = [[]]
 
list para_plot.coor = [[]]
 
 para_plot.f = open("para_"+str(snap)+".csv")
 
 para_plot.dat = f.readlines()[1:]
 
int para_plot.count = 0
 
 para_plot.tem = list(map(float,line.split(',')))
 
 para_plot.sample_point = len(rho[0])
 
int para_plot.origin_point = int(sample_point/2)+1
 
 para_plot.q_id = args.pq
 
 para_plot.coor_data = np.array(coor[snapshot])
 
 para_plot.e_data = np.array(E_t[snapshot])
 
 para_plot.m1_data = np.array(m_1[snapshot])
 
 para_plot.m2_data = np.array(m_2[snapshot])
 
 para_plot.m3_data = np.array(m_3[snapshot])
 
 para_plot.rho_data = np.array(rho[snapshot])
 
 para_plot.q_data = np.array(rho[snapshot])
 
list para_plot.t_code = snapshots[snapshot]
 
tuple para_plot.t_real = (0.1694*a_i**(-0.5)-t_code/11.8)**(-3)*10**15
 
float para_plot.a_scale = 0.0287*(0.1694*a_i**(-0.5)-t_code/11.8)**(-2)
 
int para_plot.coor_data_physical = 100*a_scale*coor_data
 
str para_plot.q_id_show = "$"+q_id+"/(km/s)$"
 
int para_plot.v_hubble = 100*(2.057*10**17/t_real)*coor_data_physical
 
float para_plot.r_ta = a_i*100*3.086*1e19*((3*M)/(4*pi*rho_c_ini))**(1/3)*(4/3/pi)**(8/9)*(t_real/6.505/1e12)**(8/9)
 
int para_plot.rho_scaling = 1/(6*pi*G*t_real**2)
 
float para_plot.v_scaling = r_ta/t_real
 
int para_plot.P_scaling = rho_scaling*r_ta**2*1e6/(t_real**2)
 
 para_plot.label