# ---
# jupyter:
#   jupytext:
#     text_representation:
#       extension: .py
#       format_name: light
#       format_version: '1.5'
#       jupytext_version: 1.14.4
#   kernelspec:
#     display_name: Python 3
#     language: python
#     name: python3
# ---

# # $M \propto R\sigma^2$
#
# Start from GalaxyMaker save files.

import load

pwd

gal = load.rd_GM.rd_gal(187, 10, base='./29172/')

import tree.halomodule as hmo

t = hmo.Halo(base='./29176/', nout=187, is_gal = True, load=True)

t.data.dtype


import matplotlib.pyplot as plt

t.data['sig']

plt.scatter(t.data['m'], t.data['sigbulge']**2 * t.data['radius'][0])
plt.show()



StarcoderdataJupyterScriptsDedupFiltered