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)
#
# Mass comparison of simulation to PV data
#
if ((raw_fit == "FALSE") and (mcomp_tmp == "TRUE")):
mass_comp.comp(arch_dir_tmp, plotdir, ea_run[i], \
snaparr, timearr, v_K, inclin, m_s, m_mri_d, pmass, \
m_d_kep, m_d_piv, m_d_sigALMA, pv_mass, EA_lenref)
#
# Convert all images to .eps format, suitable for publication
#
os.system('for f in '+plotdir+'*.pdf; do pdftops -eps $f; done')
os.system('rm -r '+plotdir+'*.pdf')
#
#
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)
#
# Mass comparison of simulation to PV data
#
if ((raw_fit == "FALSE") and (kep_fit == "TRUE")
and (mcomp_tmp == "TRUE")):
mass_comp.comp(dat_dir_tmp, plotdir, ea_run[i], snaparr, \
snapcore, timearr, v_K, inclin, m_s, m_mri_d, pmass, \
m_d_kep, m_d_piv, m_d_sigALMA, pv_mass, EA_lenref)
#
# Convert all images to .eps format, suitable for publication
#
os.system('for f in ' + plotdir +
'*.pdf; do pdftops -eps $f; done')
os.system('rm -r ' + plotdir + '*.pdf')
#
# DE05.sink file(s) read (Planet sink data vs. radius vs. time)
#
pmass, pradius, timearr = sink.read(arch_dir, ic_dir, plotdir, \
ea_run[i], snaparr, pcAU)
#
# pdisc read (Disc parameters vs. radius for defined time)
#
r, vkep = pdisc.read(arch_dir, dat_dir_tmp, snaparr)
#
# Position-Velocity diagram plot
#
if ((pv == "TRUE")):
pv_mass60, pv_mass90, kep_fit, raw_fit = pv_diag.pv(dat_dir_tmp, \
plotdir, ea_run[i], snaparr, v_K[k], inclin[j], r, vkep, pmass, \
EA_lenref, EA_timeref, pcAU, AUm, G, Msol_kg, pv_MASK60, \
pv_MASK90, pcm, mH, plot_form, kb, c, h, cgsSIopac, mu, Jerg, \
cm2m2, part_mass )
#
if (inclin == "60"):
pv_MASK60 = pv_mass60
#
# Mass comparison of simulation to PV data (both 60 and 90 deg.)
#
if ((raw_fit == "FALSE") and (kep_fit == "TRUE")
and (mcomp_tmp == "TRUE")):
mass_comp.comp(dat_dir_tmp, plotdir, ea_run[i], snaparr, \
timearr, v_K, inclin, m_s, m_mri_d, pmass, m_d_kep, \
m_d_piv, m_d_sigALMA, pv_mass60, pv_mass90, EA_lenref, plot_form)
#
# Convert all images to .eps format, suitable for publication
#
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,
)
#
# Mass comparison of simulation to PV data
#
if (raw_fit == "FALSE") and (kep_fit == "TRUE") and (mcomp_tmp == "TRUE"):
mass_comp.comp(
dat_dir_tmp,
plotdir,