import matplotlib.pyplot as plt
import isochrones as iso
import numpy as np
log_age = 9.
zp01 = iso.read("zp01.dat")
zp1_temp = []
zp1_mbol = []
zp1_intimf = []
zp1_diff = [0]
for i in range(0, len(zp01['age'])):
if zp01['age'][i] == log_age:
zp1_temp.append(zp01['logte'][i])
zp1_mbol.append(zp01['mbol'][i])
zp1_intimf.append(zp01['intimf'][i])
for j in range(1, len(zp1_intimf)):
zp1_diff.append(10**(abs(zp1_intimf[j] - zp1_intimf[j - 1])))
fig, ax2 = plt.subplots()
number_bins = 30
xbin_width = (max(zp1_temp) - min(zp1_temp)) / float(number_bins)
ybin_width = (max(zp1_mbol) - min(zp1_temp)) / float(number_bins)
xedges = [min(zp1_temp)]
yedges = [min(zp1_mbol)]
for j in range(1, 30):
xedges.append(xedges[j - 1] + xbin_width)
import matplotlib.pyplot as plt
import isochrones as iso
import numpy as np
zp00 = iso.read("zp00.dat")
zp02 = iso.read("zp02.dat")
zp04 = iso.read("zp04.dat")
agetest00 = zp00['age']
agetest02 = zp02['age']
agetest04 = zp04['age']
log_age = 9.
'''
zp0 = {'temp' : zp00['logte'], 'grav' : zp00['logg'], 'mbol' : zp00['mbol'], 'z' : zp00['z'], 'age' : zp00['age']}
zp2 = {'temp' : zp02['logte'], 'grav' : zp02['logg'], 'mbol' : zp02['mbol'], 'z' : zp02['z'], 'age' : zp02['age']}
zp4 = {'temp' : zp04['logte'], 'grav' : zp04['logg'], 'mbol' : zp04['mbol'], 'z' : zp04['z'], 'age' : zp04['age']}
'''
zp0_temp = []
zp0_grav = []
zp0_mbol = []
zp2_temp = []
zp2_grav = []
zp2_mbol = []
zp4_temp = []
zp4_grav = []
import matplotlib.pyplot as plt
import isochrones as iso
import numpy as np
zp00 = iso.read("zp00.dat")
zp02 = iso.read("zp02.dat")
zp04 = iso.read("zp04.dat")
agetest00 = zp00['age']
agetest02 = zp02['age']
agetest04 = zp04['age']
log_age = 9.
'''
zp0 = {'temp' : zp00['logte'], 'grav' : zp00['logg'], 'mbol' : zp00['mbol'], 'z' : zp00['z'], 'age' : zp00['age']}
zp2 = {'temp' : zp02['logte'], 'grav' : zp02['logg'], 'mbol' : zp02['mbol'], 'z' : zp02['z'], 'age' : zp02['age']}
zp4 = {'temp' : zp04['logte'], 'grav' : zp04['logg'], 'mbol' : zp04['mbol'], 'z' : zp04['z'], 'age' : zp04['age']}
'''
zp0_temp = []
zp0_grav = []
zp0_mbol = []
zp2_temp = []
zp2_grav = []
zp2_mbol = []
zp4_temp = []
zp4_grav = []
import matplotlib.pyplot as plt
import isochrones as iso
import numpy as np
log_age = 9.
zp01 = iso.read("zp01.dat")
zp1_temp = []
zp1_mbol = []
zp1_intimf = []
zp1_diff = [0]
for i in range(0, len(zp01['age'])):
if zp01['age'][i] == log_age:
zp1_temp.append(zp01['logte'][i])
zp1_mbol.append(zp01['mbol'][i])
zp1_intimf.append(zp01['intimf'][i])
for j in range(1, len(zp1_intimf)):
zp1_diff.append(10 ** (abs(zp1_intimf[j] - zp1_intimf[j-1])))
fig, ax2 = plt.subplots()
number_bins = 30
xbin_width = (max(zp1_temp) - min(zp1_temp)) / float(number_bins)
ybin_width = (max(zp1_mbol) - min(zp1_temp)) / float(number_bins)
xedges = [min(zp1_temp)]
yedges = [min(zp1_mbol)]
for j in range(1, 30):
