def read(snapshot_number): _, path_to_catalogue = find_files() path_to_catalogue.replace(f'{xlargs.snapshot_number:04d}', f'{snapshot_number:04d}') vr_handle = velociraptor.load(path_to_catalogue) r500 = vr_handle.spherical_overdensities.r_500_rhocrit[0].to('Mpc').value xcminpot = vr_handle.positions.xcminpot[0].to('Mpc').value ycminpot = vr_handle.positions.ycminpot[0].to('Mpc').value zcminpot = vr_handle.positions.zcminpot[0].to('Mpc').value return r500, xcminpot, ycminpot, zcminpot
import sys import os.path from matplotlib import pyplot as plt import numpy as np from unyt import Solar_Mass import traceback sys.path.append("..") from scaling_relations import EntropyFgasSpace, EntropyProfiles from register import find_files, set_mnras_stylesheet, xlargs, default_output_directory from literature import Sun2009, Pratt2010, Croston2008, Cosmology snap, cat = find_files() kwargs = dict(path_to_snap=snap, path_to_catalogue=cat) set_mnras_stylesheet() try: cosmology = Cosmology() gas_profile_obj = EntropyFgasSpace(max_radius_r500=1.) radial_bin_centres, cumulative_gas_mass_profile, cumulative_mass_profile, m500fb = gas_profile_obj.process_single_halo( **kwargs) entropy_profile_obj = EntropyProfiles(max_radius_r500=1) _, entropy_profile, K500 = entropy_profile_obj.process_single_halo( **kwargs) entropy_profile /= K500 gas_fraction_enclosed = cumulative_gas_mass_profile / m500fb
""" import sys import copy import os import numpy as np from matplotlib import pyplot as plt from matplotlib.colors import LogNorm sys.path.append("..") from scaling_relations import SliceGas from register import default_output_directory, find_files, xlargs, set_mnras_stylesheet from literature import Cosmology s, c = find_files() set_mnras_stylesheet() centres = np.load( 'map_centre_L0300N0564_VR18_-8res_MinimumDistance_fixedAGNdT8.5_Nheat1_alpha1p0.npy' ) import matplotlib.offsetbox from matplotlib.lines import Line2D class AnchoredHScaleBar(matplotlib.offsetbox.AnchoredOffsetbox): """ size: length of bar in data units extent : height of bar ends in axes units """