def plot_chi_uds_x(dirnames, ax):
chis, uds, uds_err = utils.chi_uds_x(dirnames)
i_sort = np.argsort(chis)
chis, uds = chis[i_sort], uds[i_sort]
ax.plot(chis, uds)
ax.set_xlim(-0.02, 1.0)
ax.set_ylim(0.0, 1.1)
use_pgf = True
ejm_rcparams.set_pretty_plots(use_latex, use_pgf)
fig = plt.figure(figsize=(12, 12 * ejm_rcparams.golden_ratio))
ax = fig.add_subplot(111)
ejm_rcparams.prettify_axes(ax)
dirnames_1D_p_1side_S = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_1D,dt=0.01,seed=1,n=5000,align=(0),origin=(1),v=20,p=1,chi=*,side=1,type=S,obs=NoObs,c=NoC')
dirnames_1D_p_2side_S = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_1D,dt=0.01,seed=1,n=5000,align=(0),origin=(1),v=20,p=1,chi=*,side=2,type=S,obs=NoObs,c=NoC')
dirnames_1D_p_1side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_1D,dt=0.01,seed=1,n=5000,align=(0),origin=(1),v=20,p=1,chi=*,side=1,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
dirnames_1D_p_2side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_1D,dt=0.01,seed=1,n=5000,align=(0),origin=(1),v=20,p=1,chi=*,side=2,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
chis_th = np.linspace(0.0, 0.99, 100)
chis, uds, uds_err = utils.chi_uds_x(dirnames_1D_p_2side_S)
i_sort = np.argsort(chis)
chis, uds, uds_err = chis[i_sort], uds[i_sort], uds_err[i_sort]
ax.errorbar(chis, uds, yerr=uds_err,
label='Spatial, two-sided', c=ejm_rcparams.set2[0])
ax.plot(chis_th, chis_th, label=r'$u_x = \chi$',
ls='dashed', c=ejm_rcparams.set2[0])
chis, uds, uds_err = utils.chi_uds_x(dirnames_1D_p_1side_S)
i_sort = np.argsort(chis)
chis, uds, uds_err = chis[i_sort], uds[i_sort], uds_err[i_sort]
ax.errorbar(chis, uds, yerr=uds_err,
label='Spatial, one-sided', c=ejm_rcparams.set2[1])
ax.plot(chis_th, chis_th / (2.0 - chis_th), label=r'$u_x = \chi / (2 - \chi)$',
ls='dashed', c=ejm_rcparams.set2[1])
use_latex = save_flag
use_pgf = True
ejm_rcparams.set_pretty_plots(use_latex, use_pgf)
fig = plt.figure(figsize=(12, 12 * ejm_rcparams.golden_ratio))
ax = fig.add_subplot(111)
ejm_rcparams.prettify_axes(ax)
dirnames_2D_p_1side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_2D,dt=0.01,seed=1,n=5000,align=(0,0),origin=(1,1),v=20,p=1,chi=*,side=1,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
dirnames_2D_p_2side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_2D,dt=0.01,seed=1,n=5000,align=(0,0),origin=(1,1),v=20,p=1,chi=*,side=2,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
dirnames_2D_Dr_1side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_2D,dt=0.01,seed=1,n=5000,align=(0,0),origin=(1,1),v=20,Dr=1,chi=*,side=1,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
dirnames_2D_Dr_2side_T = glob('/Users/ewj/Desktop/porous/ahoy_data/chi_scan/ahoy_2D,dt=0.01,seed=1,n=5000,align=(0,0),origin=(1,1),v=20,Dr=1,chi=*,side=2,type=T,dtMem=0.1,tMem=5,obs=NoObs,c=NoC')
chis, uds, uds_err = utils.chi_uds_x(dirnames_2D_p_2side_T)
i_sort = np.argsort(chis)
chis, uds, uds_err = chis[i_sort], uds[i_sort], uds_err[i_sort]
ax.errorbar(chis, uds, yerr=uds_err,
label='Tumbling, two-sided', c=ejm_rcparams.set2[0])
chis, uds, uds_err = utils.chi_uds_x(dirnames_2D_p_1side_T)
i_sort = np.argsort(chis)
chis, uds, uds_err = chis[i_sort], uds[i_sort], uds_err[i_sort]
ax.errorbar(chis, uds, yerr=uds_err,
label='Tumbling, one-sided', c=ejm_rcparams.set2[1])
chis, uds, uds_err = utils.chi_uds_x(dirnames_2D_Dr_2side_T)
i_sort = np.argsort(chis)
chis, uds, uds_err = chis[i_sort], uds[i_sort], uds_err[i_sort]
ax.errorbar(chis, uds, yerr=uds_err,
