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 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):