# Read in rdisc.1 file
#
hasharr_app, n_accr, r_d_kep, r_d_sigALMA, \
m_s, m_mri_d, m_d_kep, m_d_piv, m_d_sigALMA, = \
r1_r.read(dat_dir_tmp, plotdir, ea_run[i], \
snaparr, v_K[k], inclin[j])
#
# DE05.sink file(s) read (Planet sink data vs. radius vs. time)
#
pmass, pradius, timearr = \
sink_r.read(arch_dir, dat_dir_tmp, ic_dir, plotdir, ea_run[i], snaparr, pcAU)
#
# pdisc read (Disc parameters vs. radius for defined time)
#
r, vkep, col_arr = pdisc_r.read(arch_dir, dat_dir_tmp, plotdir, ea_run[i], \
snaparr, snaparr_IA, snapcore, EA_timeref, EA_lenref, pmass, pradius, \
hasharr_app, n_accr, r_limit, r_d_kep, r_d_sigALMA,
timearr, v_K[k], inclin[j])
#
# rdisc read (Disc parameters vs. time for defined radius)
#
rdisc_r.read(dat_dir_tmp, plotdir, ea_run[i], hasharr_app, \
n_accr, r_inspec, v_K[k], inclin[j])
#
if ((pv == "TRUE")): # and (inclin[j] != "0")):
#
# Position-Velocity diagram plot
#
pv_mass, kep_fit, raw_fit = pv_diag.pv(dat_dir_tmp, plotdir, \
ea_run[i], snaparr, snapcore, v_K[k], inclin[j], r, vkep, \
pmass, EA_lenref, EA_timeref, pcAU, AUm, G, Msol_kg)
#
# Read in rdisc.1 file # hasharr_app, n_accr, r_d_kep, r_d_sigALMA, \ m_s, m_mri_d, m_d_kep, m_d_sigALMA, m_d_piv = \ r1_r.read(arch_dir_tmp, plotdir, ea_run[i], \ snaparr, v_K[k], inclin[j]) # # DE05.sink file(s) read (Planet sink data vs. radius vs. time) # pmass, pradius, timearr = \ sink_r.read(arch_dir_tmp, plotdir, ea_run[i], snaparr, pcAU) # # pdisc read (Disc parameters vs. radius for defined time) # r, vkep = pdisc_r.read(arch_dir_tmp, plotdir, ea_run[i], \ snaparr, EA_timeref, EA_lenref, pmass, pradius, \ hasharr_app, n_accr, r_limit, spline, r_d_kep, r_d_sigALMA, timearr, v_K[k]) # # rdisc read (Disc parameters vs. time for defined radius) # rdisc_r.read(arch_dir_tmp, plotdir, ea_run[i], hasharr_app, \ n_accr, r_inspec, v_K[k], inclin[j]) # if ((pv == "TRUE") and (inclin[j] != "0")): # # Position-Velocity diagram plot # pv_mass, raw_fit = pv_diag.pv(arch_dir_tmp, plotdir, ea_run[i], \ snaparr, v_K[k], inclin[j], r, vkep, EA_lenref, EA_timeref, \ pcAU, AUm, G, Msol_kg) #
# DE05.sink file(s) read (Planet sink data vs. radius vs. time)
#
pmass, pradius, timearr = sink_r.read(arch_dir, dat_dir_tmp, plotdir, ea_run[i], snaparr, pcAU)
#
# pdisc read (Disc parameters vs. radius for defined time)
#
r, vkep, col_arr = pdisc_r.read(
dat_dir_tmp,
plotdir,
ea_run[i],
snaparr,
snapcore,
EA_timeref,
EA_lenref,
pmass,
pradius,
hasharr_app,
n_accr,
r_limit,
r_d_kep,
r_d_sigALMA,
timearr,
v_K[k],
inclin[j],
)
#
# rdisc read (Disc parameters vs. time for defined radius)
#
rdisc_r.read(dat_dir_tmp, plotdir, ea_run[i], hasharr_app, n_accr, r_inspec, v_K[k], inclin[j])
#
if (pv == "TRUE") and (inclin[j] != "0"):