import mppnp as mp
import utils
import mesa as ms

'''
logdir='/rpod2/athira/mesa/case1/LOGS'
p=ms.mesa_profile(logdir,100,num_type='model')
cycles_1=p.log_ind.keys()
cycles_1.sort()
m=mp.se('/rpod2/athira/mppnp/case1/H5_out')

#specs1a_l=['$^1H$','$^4He$','$^{12}C$','$^{13}C$','$^{14}N$','$^{15}N$','$^{16}O$','$^{18}O$','$^{19}F$','$^{22}Ne$']

specs1a_l=['$^1H$','$^4He$','$^{12}C$','$^{13}C$','$^{14}N$']
specs1b_l=['$^{15}N$','$^{16}O$','$^{18}O$','$^{19}F$','$^{22}Ne$']

lines=['b--','ro-','bs-','k-.','r+:','bd-','kx--','r>-','kp-','r.:','k^-']

matplotlib.rc('xtick',labelsize=13,color='k')
matplotlib.rc('ytick',labelsize=10)
matplotlib.rc('ytick.major',pad=10)
matplotlib.rc('axes',labelsize=12,labelcolor='k')
'''
ind=1

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a_l=['$^4He$','$^{12}C$','$^{14}N$','$^{16}O$','$^{22}Ne$']

specs1a=[he4,c12,n14,o16,ne22]

close(1)
ifig=figure(1,figsize=(9,5))

ax1=ifig.add_subplot(1,1,1)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])


for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.80,markevery=15)

ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
ax1.set_xlabel('$Mr/M_\mathrm{\odot}$')
ylim(1e-3,1.3)
xlim(0.0,0.90)
text(0.25,8e-3,"Core",rotation=90,color='k',fontsize=10)
text(0.462,4e-3,"Buffer",rotation=90,color='k',fontsize=10)
text(0.535,8e-3,"SOF",rotation=90,color='k',fontsize=10)
text(0.75,8e-3,"Envelope",rotation=90,color='k',fontsize=10)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)

'''
close(40)
ifig=figure(40,figsize=(10,12))
subplot(2,2,1)

ax1=ifig.add_subplot(2,2,1)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1

specs2_l=['$Na$','$Mg$','$Al$','$Si$','$S$','$Ca$','$Ti$']

na23=m.get(ind,'Na-23');
mg24=m.get(ind,'Mg-24');mg25=m.get(ind,'Mg-25');mg26=m.get(ind,'Mg-26')
al27=m.get(ind,'Al-27')
si28=m.get(ind,'Si-28');si29=m.get(ind,'Si-29');si30=m.get(ind,'Si-30')
s32=m.get(ind,'S-32');s33=m.get(ind,'S-33');s34=m.get(ind,'S-34');s36=m.get(ind,'S-36')
ca40=m.get(ind,'Ca-40');ca42=m.get(ind,'Ca-42');ca43=m.get(ind,'Ca-43');ca44=m.get(ind,'Ca-44');ca46=m.get(ind,'Ca-46');ca48=m.get(ind,'Ca-48')
ti46=m.get(ind,'Ti-46');ti47=m.get(ind,'Ti-47');ti48=m.get(ind,'Ti-48');ti49=m.get(ind,'Ti-49');ti50=m.get(ind,'Ti-50')

na=na23/23.0
mg=mg24/24.0+mg25/25.0+mg26/26.0
al=al27/27.0
si=si28/28.0+si29/29.0+si30/30.0
s=s32/32.0+s33/33.0+s34/34.0+s36/36.0
ca=ca40/40.0+ca42/42.0+ca43/43.0+ca44/44.0+ca46/46.0+ca48/48.0
ti=ti46/46.0+ti47/47.0+ti48/48.0+ti49/49.0+ti50/50.0

specs2=[na,mg,al,si,s,ca,ti]
mass=m.get(ind,'mass')

i=0
for a in range (len(specs2)):
    ax1.semilogy(mass,specs2[a],lines[i], label=specs2_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
setp(ax1.get_xticklabels(), visible=False)

ylim(1e-6,1e-3)
xlim(0.44,0.90)
text(0.35, 1e-3, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

subplot(2,2,2)
ax1=ifig.add_subplot(2,2,2)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1

specs3_l=['$Fe$','$Ni$','$Zn$','$Y$','$Zr$','$Ba$','$La$']

ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
setp(ax1.get_xticklabels(), visible=False)


fe56=m.get(ind,'Fe-56');fe57=m.get(ind,'Fe-57');fe58=m.get(ind,'Fe-58')
ni58=m.get(ind,'Ni-58');ni60=m.get(ind,'Ni-60');ni61=m.get(ind,'Ni-61');ni62=m.get(ind,'Ni-62');ni64=m.get(ind,'Ni-64')
zn64=m.get(ind,'Zn-64');zn66=m.get(ind,'Zn-66');zn67=m.get(ind,'Zn-67');zn68=m.get(ind,'Zn-68');zn70=m.get(ind,'Zn-70')
y89=m.get(ind,'Y-89')
zr90=m.get(ind,'Zr-90');zr91=m.get(ind,'Zr-91');zr92=m.get(ind,'Zr-92');zr94=m.get(ind,'Zr-94');zr96=m.get(ind,'Zr-96')
ba130=m.get(ind,'Ba-130');ba132=m.get(ind,'Ba-132');ba134=m.get(ind,'Ba-134');ba135=m.get(ind,'Ba-135');ba136=m.get(ind,'Ba-136');ba137=m.get(ind,'Ba-137');ba138=m.get(ind,'Ba-138')
la138=m.get(ind,'La-138');la139=m.get(ind,'La-139')

fe=fe56/56.0+fe57/57.0+fe58/58.0
ni=ni58/58.0+ni60/60.0+ni61/61.0+ni62/62.0+ni64/64.0
zn=zn64/64.0+zn66/66.0+zn67/67.0+zn68/68.0+zn70/70.0
y=y89/89.0
zr=zr90/90.0+zr91/91.0+zr92/92.0+zr94/94.0+zr96/96.0
ba=ba130/130.0+ba132/132.0+ba134/134.0+ba135/135.0+ba136/136.0+ba137/137.0+ba138/138.0
la=la138/138.0+la139/139.0

specs3=[fe,ni,zn,y,zr,ba,la]

for a in range (len(specs3)):
    ax1.semilogy(mass,specs3[a],lines[a], label=specs3_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

setp(ax1.get_xticklabels(), visible=False)
ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})

#ylim(1e-6,1e-3)
xlim(0.44,0.90)
text(0.95, 1e-4, 'b', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)


ax1=ifig.add_subplot(2,2,3)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2500

specs2_l=['$Na$','$Mg$','$Al$','$Si$','$S$','$Ca$','$Ti$']

na23=m.get(ind,'Na-23');
mg24=m.get(ind,'Mg-24');mg25=m.get(ind,'Mg-25');mg26=m.get(ind,'Mg-26')
al27=m.get(ind,'Al-27')
si28=m.get(ind,'Si-28');si29=m.get(ind,'Si-29');si30=m.get(ind,'Si-30')
s32=m.get(ind,'S-32');s33=m.get(ind,'S-33');s34=m.get(ind,'S-34');s36=m.get(ind,'S-36')
ca40=m.get(ind,'Ca-40');ca42=m.get(ind,'Ca-42');ca43=m.get(ind,'Ca-43');ca44=m.get(ind,'Ca-44');ca46=m.get(ind,'Ca-46');ca48=m.get(ind,'Ca-48')
ti46=m.get(ind,'Ti-46');ti47=m.get(ind,'Ti-47');ti48=m.get(ind,'Ti-48');ti49=m.get(ind,'Ti-49');ti50=m.get(ind,'Ti-50')

na=na23/23.0
mg=mg24/24.0+mg25/25.0+mg26/26.0
al=al27/27.0
si=si28/28.0+si29/29.0+si30/30.0
s=s32/32.0+s33/33.0+s34/34.0+s36/36.0
ca=ca40/40.0+ca42/42.0+ca43/43.0+ca44/44.0+ca46/46.0+ca48/48.0
ti=ti46/46.0+ti47/47.0+ti48/48.0+ti49/49.0+ti50/50.0

specs2=[na,mg,al,si,s,ca,ti]
mass=m.get(ind,'mass')

i=0
for a in range (len(specs2)):
    ax1.semilogy(mass,specs2[a],lines[i], label=specs2_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

#ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-6,1e-3)
text(0.35, 1e-3, 'c', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

neut_den_1=[]
neut=m.get(ind,'Neutron-1')
rho= m.get(ind,'rho')
mass=m.get(ind,'mass')

for a in range(len(neut)):
    res=(neut[a]*6*10**23*rho[a])
    neut_den_1.append(res)

ax2=ax1.twinx()
ax2.plot(mass,log10(neut_den_1),'orange',linestyle='-',linewidth=1.2)
ylim(0,max(log10(neut_den_1)+1))

#setp(ax1.get_xticklabels(), visible=False)

xlim(0.44,0.90)


subplot(2,2,4)
ax1=ifig.add_subplot(2,2,4)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2500

specs3_l=['$Fe$','$Ni$','$Zn$','$Y$','$Zr$','$Ba$','$La$']

ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
#setp(ax1.get_xticklabels(), visible=False)

xlim(0.44,0.90)

fe56=m.get(ind,'Fe-56');fe57=m.get(ind,'Fe-57');fe58=m.get(ind,'Fe-58')
ni58=m.get(ind,'Ni-58');ni60=m.get(ind,'Ni-60');ni61=m.get(ind,'Ni-61');ni62=m.get(ind,'Ni-62');ni64=m.get(ind,'Ni-64')
zn64=m.get(ind,'Zn-64');zn66=m.get(ind,'Zn-66');zn67=m.get(ind,'Zn-67');zn68=m.get(ind,'Zn-68');zn70=m.get(ind,'Zn-70')
y89=m.get(ind,'Y-89')
zr90=m.get(ind,'Zr-90');zr91=m.get(ind,'Zr-91');zr92=m.get(ind,'Zr-92');zr94=m.get(ind,'Zr-94');zr96=m.get(ind,'Zr-96')
ba130=m.get(ind,'Ba-130');ba132=m.get(ind,'Ba-132');ba134=m.get(ind,'Ba-134');ba135=m.get(ind,'Ba-135');ba136=m.get(ind,'Ba-136');ba137=m.get(ind,'Ba-137');ba138=m.get(ind,'Ba-138')
la138=m.get(ind,'La-138');la139=m.get(ind,'La-139')

fe=fe56/56.0+fe57/57.0+fe58/58.0
ni=ni58/58.0+ni60/60.0+ni61/61.0+ni62/62.0+ni64/64.0
zn=zn64/64.0+zn66/66.0+zn67/67.0+zn68/68.0+zn70/70.0
y=y89/89.0
zr=zr90/90.0+zr91/91.0+zr92/92.0+zr94/94.0+zr96/96.0
ba=ba130/130.0+ba132/132.0+ba134/134.0+ba135/135.0+ba136/136.0+ba137/137.0+ba138/138.0
la=la138/138.0+la139/139.0

specs3=[fe,ni,zn,y,zr,ba,la]

neut_den_1=[]
neut=m.get(ind,'Neutron-1')
rho= m.get(ind,'rho')
mass=m.get(ind,'mass')

for a in range (len(specs3)):
    ax1.semilogy(mass,specs3[a],lines[a], label=specs3_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})
ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
text(0.95, 1e-4, 'b', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

for a in range(len(neut)):
    res=(neut[a]*6*10**23*rho[a])
    neut_den_1.append(res)

ax2=ax1.twinx()
ax2.plot(mass,log10(neut_den_1),'orange',linestyle='-',linewidth=1.2)
#ax2.set_ylabel('Log Neutron Density (/cc)')
ylim(0,max(log10(neut_den_1)+1))

#setp(ax1.get_xticklabels(), visible=False)

xlim(0.44,0.90)
'''
'''
specs1a_l=['$^1H$','$^4He$','$^{12}C$','$^{13}C$','$^{14}N$']
specs1b_l=['$^{15}N$','$^{16}O$','$^{18}O$','$^{19}F$','$^{22}Ne$']

close(30)
ifig=figure(30,figsize=(10,12))

subplot(4,2,1)

ax1=ifig.add_subplot(4,2,1)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1

i=0
h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14]
specs1b=[n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[i], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
setp(ax1.get_xticklabels(), visible=False)

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

subplot(4,2,2)
ax1=ifig.add_subplot(4,2,2)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1

for a in range (len(specs1b)):
    ax1.semilogy(mass,specs1b[a],lines[i], label=specs1b_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

setp(ax1.get_xticklabels(), visible=False)
ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.95, 1.0, 'b', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)


ax1=ifig.add_subplot(4,2,3)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=700

i=0
h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14]
specs1b=[n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[i], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

#ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
setp(ax1.get_xticklabels(), visible=False)

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.35, 1.0, 'c', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

subplot(4,2,4)
ax1=ifig.add_subplot(4,2,4)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=700

for a in range (len(specs1b)):
    ax1.semilogy(mass,specs1b[a],lines[i], label=specs1b_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

setp(ax1.get_xticklabels(), visible=False)
#ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.95, 1.0, 'd', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

ax1=ifig.add_subplot(4,2,5)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1150

i=0
h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14]
specs1b=[n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[i], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

#ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
setp(ax1.get_xticklabels(), visible=False)

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.35, 1.0, 'e', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

subplot(4,2,6)
ax1=ifig.add_subplot(4,2,6)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1150

for a in range (len(specs1b)):
    ax1.semilogy(mass,specs1b[a],lines[i], label=specs1b_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

setp(ax1.get_xticklabels(), visible=False)
#ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.95, 1.0, 'f', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

ax1=ifig.add_subplot(4,2,7)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2500

i=0
h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14]
specs1b=[n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[i], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

#ax1.legend(loc='center', bbox_to_anchor=(-0.25, 0.36),prop={'size':9.0})
#setp(ax1.get_xticklabels(), visible=False)
ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')

ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.35, 1.0, 'g', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)

subplot(4,2,8)
ax1=ifig.add_subplot(4,2,8)
box = ax1.get_position()
#ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2500

for a in range (len(specs1b)):
    ax1.semilogy(mass,specs1b[a],lines[i], label=specs1b_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)
    i+=1

#setp(ax1.get_xticklabels(), visible=False)
#ax1.legend(loc='center', bbox_to_anchor=(1.17, 0.36),prop={'size':9.0})
ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
text(0.95, 1.0, 'h', ha='center', va='center', rotation='horizontal',color='k',fontsize=14)
'''
'''
ifig=figure(20,figsize=(10,12))

ax1=ifig.add_subplot(4,1,3)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1150

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('c')

subplot(4,1,4)
ax1=ifig.add_subplot(4,1,4)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2000

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('d')
'''

'''
close(10)
ifig=figure(10,figsize=(8.25,11.6))

subplot(3,1,1)
ax1=ifig.add_subplot(3,1,1)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.0,0.90)
text(0.25,4e-3,"Core (Hydro)",rotation=90,color='g',fontsize=9)
text(0.462,4e-3,"Buffer",rotation=90,color='g',fontsize=9)
text(0.535,1e-3,"SOF",rotation=90,color='g',fontsize=9)
text(0.75,9e-4,"Envelope",rotation=90,color='g',fontsize=9)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('a')

ifig=figure(10,figsize=(8.25,11.6))

subplot(3,2,3)
ax1=ifig.add_subplot(3,2,3)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=700

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('b')

subplot(3,2,4)
ax1=ifig.add_subplot(3,2,4)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=1150

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('c')

subplot(3,2,5)
ax1=ifig.add_subplot(3,2,5)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2000

h1=m.get(ind,'H-1');he4=m.get(ind,'He-4');c12=m.get(ind,'C-12');c13=m.get(ind,'C-13');n14=m.get(ind,'N-14');n15=m.get(ind,'N-15');o16=m.get(ind,'O-16');o18=m.get(ind,'O-18');f19=m.get(ind,'F-19');ne22=m.get(ind,'Ne-22')
mass=m.get(ind,'mass')

specs1a=[h1,he4,c12,c13,n14,n15,o16,o18,f19,ne22]

for a in range (len(specs1a)):
    ax1.semilogy(mass,specs1a[a],lines[a], label=specs1a_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

#ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Mass fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
ylim(1e-8,1.3)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('d')

subplot(3,2,6)
ax1=ifig.add_subplot(3,2,6)
box = ax1.get_position()
ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height])
ind=2000

specs3_l=['$Fe$','$Ni$','$Zn$','$Y$','$Zr$','$Ba$','$La$']

neut_den_1=[]
neut=m.get(ind,'Neutron-1')
rho= m.get(ind,'rho')
mass=m.get(ind,'mass')

for a in range(len(neut)):
    res=(neut[a]*6*10**23*rho[a])
    neut_den_1.append(res)

fe56=m.get(ind,'Fe-56');fe57=m.get(ind,'Fe-57');fe58=m.get(ind,'Fe-58')
ni58=m.get(ind,'Ni-58');ni60=m.get(ind,'Ni-60');ni61=m.get(ind,'Ni-61');ni62=m.get(ind,'Ni-62');ni64=m.get(ind,'Ni-64')
zn64=m.get(ind,'Zn-64');zn66=m.get(ind,'Zn-66');zn67=m.get(ind,'Zn-67');zn68=m.get(ind,'Zn-68');zn70=m.get(ind,'Zn-70')
y89=m.get(ind,'Y-89')
zr90=m.get(ind,'Zr-90');zr91=m.get(ind,'Zr-91');zr92=m.get(ind,'Zr-92');zr94=m.get(ind,'Zr-94');zr96=m.get(ind,'Zr-96')
ba130=m.get(ind,'Ba-130');ba132=m.get(ind,'Ba-132');ba134=m.get(ind,'Ba-134');ba135=m.get(ind,'Ba-135');ba136=m.get(ind,'Ba-136');ba137=m.get(ind,'Ba-137');ba138=m.get(ind,'Ba-138')
la138=m.get(ind,'La-138');la139=m.get(ind,'La-139')

fe=fe56/56.0+fe57/57.0+fe58/58.0
ni=ni58/58.0+ni60/60.0+ni61/61.0+ni62/62.0+ni64/64.0
zn=zn64/64.0+zn66/66.0+zn67/67.0+zn68/68.0+zn70/70.0
y=y89/89.0
zr=zr90/90.0+zr91/91.0+zr92/92.0+zr94/94.0+zr96/96.0
ba=ba130/130.0+ba132/132.0+ba134/134.0+ba135/135.0+ba136/136.0+ba137/137.0+ba138/138.0
la=la138/138.0+la139/139.0

specs3=[fe,ni,zn,y,zr,ba,la]

for a in range (len(specs3)):
    ax1.semilogy(mass,specs3[a],lines[a], label=specs3_l[a],linewidth=2,markersize=5,alpha=0.60,markevery=15)

ax1.legend(loc='center', bbox_to_anchor=(1.14, 0.42),prop={'size':9.0})
ax1.set_ylabel('Number fraction')
#ax1.set_xlabel('$\mathrm{Mr/M_{\odot}}$')
#ylim(1e-6,1e-4)
xlim(0.44,0.90)
#text(1.0, 1.0, '1a', ha='center', va='center', rotation='horizontal',color='k',fontsize=12)
title('e')
'''