Commit 4c9bf063 authored by Maxime Rey's avatar Maxime Rey
Browse files

Improvement: added log option for 'rad_prof'.

parent a39224d9
......@@ -564,7 +564,7 @@ def half_mass_radius(genpath, folders, labels, timestep, var_str='gas'):
ax2.set_xlabel('Redshift')
def rad_prof(genpath, folders, labels, timesteps, var_str='H', overwhat='V', r_max=150, factor=2, nbins=100, split_flow=False, xlims=None, ylims=None, saveinfile=False, savefig=False):
def rad_prof(genpath, folders, labels, timesteps, var_str='H', overwhat='V', r_max=150, factor=2, nbins=100, split_flow=False, xlims=None, ylims=None, xlog=False, saveinfile=False, savefig=False):
"""
Plots the average density of the CGM and it's average mass-weighted temperature through radial chunks. It is also possible to remove satellites from the selection.
Stacks timesteps up to the rvir of the first timestep (the smaller one).
......@@ -725,6 +725,8 @@ def rad_prof(genpath, folders, labels, timesteps, var_str='H', overwhat='V', r_m
ax.legend(frameon=False)
ax.set_yscale('log')
if xlog:
ax.set_xscale('log')
ax.set_xlabel('$r$ [kpc]')
ax.set_ylabel(ylabel)
ax.set_xlim(xlims)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment