from numpy import * import matplotlib from matplotlib.pylab import * import ppn import mppnp as mp import data_plot import utils import mesa as ms matplotlib.rc('xtick',labelsize=10,color='k') matplotlib.rc('ytick',labelsize=10) yticks(rotation = 'horizontal') matplotlib.rc('ytick.major',pad=3.5,size=4) matplotlib.rc('xtick.major',pad=4,size=4) matplotlib.rc('axes',labelsize=10,labelcolor='black') close(70) ifig=figure(70,figsize=(8.25,11.6)) element=['H','Li','C','N','O','F','Na','Mg','Al','Si','S','Ca','Ti','Fe','Ni','Zn','Y','Zr','Ba','La'] z=[1,3,6,7,8,9,11,12,13,14,16,20,22,26,28,30,39,40,56,57] X_Fe1_1=array([ -1.13500000e+01, -1.53200000e+01, 1.29000000e+00,1.54000000e+00, 1.43000000e+00, 2.06000000e+00, 8.40000000e-01, -0.00000000e+00, -4.00000000e-02, -4.00000000e-02, 1.00000000e-02, -6.00000000e-02, 3.00000000e-02, 0.00000000e+00, -7.00000000e-02, 1.50000000e-01, 5.00000000e-02, 1.00000000e-02, 2.20000000e-01, 1.00000000e-02]) X_Fe1_2=array([ -1.25800000e+01, -1.58400000e+01, 1.50000000e+00, 1.41000000e+00, 1.49000000e+00, 2.33000000e+00, 8.60000000e-01, 1.00000000e-02, -5.00000000e-02, -4.00000000e-02, 1.00000000e-02, -6.00000000e-02, 3.00000000e-02, 0.00000000e+00, -7.00000000e-02, 1.60000000e-01, 6.00000000e-02, 1.00000000e-02, 2.40000000e-01, 1.00000000e-02]) X_S1_1=array([ -1.13600000e+01, -1.53300000e+01, 1.28000000e+00, 1.53000000e+00, 1.42000000e+00, 2.05000000e+00, 8.30000000e-01, -1.00000000e-02, -5.00000000e-02, -5.00000000e-02, 0.00000000e+00, -7.00000000e-02, 2.00000000e-02, -1.00000000e-02, -8.00000000e-02, 1.40000000e-01, 4.00000000e-02, -0.00000000e+00, 2.10000000e-01, -0.00000000e+00]) X_S1_2=array([ -1.25900000e+01, -1.58500000e+01, 1.49000000e+00, 1.40000000e+00, 1.48000000e+00, 2.32000000e+00, 8.50000000e-01, -0.00000000e+00, -5.00000000e-02, -5.00000000e-02, 0.00000000e+00, -7.00000000e-02,2.00000000e-02, -1.00000000e-02, -8.00000000e-02, 1.50000000e-01, 5.00000000e-02, -0.00000000e+00, 2.30000000e-01, -0.00000000e+00]) X_Fe1 = array([ -1.16200000e+01, -1.55100000e+01, 1.41000000e+00, 1.48000000e+00, 1.46000000e+00, 2.21000000e+00, 8.50000000e-01, 0.00000000e+00, -4.00000000e-02, -4.00000000e-02, 1.00000000e-02, -6.00000000e-02, 3.00000000e-02, 0.00000000e+00, -7.00000000e-02, 1.60000000e-01, 6.00000000e-02, 1.00000000e-02, 2.30000000e-01, 1.00000000e-02]) X_S1=array([ -1.16300000e+01, -1.55200000e+01, 1.40000000e+00, 1.47000000e+00, 1.45000000e+00, 2.21000000e+00, 8.40000000e-01, -1.00000000e-02, -5.00000000e-02, -5.00000000e-02, 0.00000000e+00, -7.00000000e-02, 2.00000000e-02, -1.00000000e-02, -8.00000000e-02, 1.50000000e-01, 5.00000000e-02, -0.00000000e+00, 2.20000000e-01, -0.00000000e+00]) X_Fe1_rcb=array([ -1.07600000e+01, -1.51200000e+01, 1.17000000e+00, 1.58000000e+00, 1.40000000e+00, 1.93000000e+00, 8.40000000e-01, -0.00000000e+00, -4.00000000e-02, -4.00000000e-02, 1.00000000e-02, -6.00000000e-02, 3.00000000e-02, 0.00000000e+00, -6.00000000e-02, 1.40000000e-01, 5.00000000e-02, 1.00000000e-02, 2.20000000e-01, 1.00000000e-02]) X_S1_rcb=array([ -1.07700000e+01, -1.51300000e+01, 1.16000000e+00, 1.57000000e+00, 1.39000000e+00, 1.92000000e+00, 8.30000000e-01, -1.00000000e-02, -5.00000000e-02, -5.00000000e-02, 0.00000000e+00, -7.00000000e-02, 2.00000000e-02, -1.00000000e-02, -7.00000000e-02, 1.30000000e-01, 4.00000000e-02, -0.00000000e+00, 2.10000000e-01, -0.00000000e+00]) X_Fe1_env=array([-1.83, -2.96, 1.01, 0.92, 1.09, -1.16, 0.71, -0.02, -0.03, -0.05, 0.01, -0.04, 0.02, 0. , -0.05, 0.05, 0.03, 0.01, 0.11, -0. ]) X_S1_env=array([-1.84, -2.97, 1. , 0.91, 1.08, -1.17, 0.7 , -0.03, -0.04, -0.06, 0. , -0.05, 0.01, -0.01, -0.06, 0.04, 0.02, -0. , 0.1 , -0.01]) E_x1_1= array([ 0.55953046, -14.21530996, 9.7946122 , 9.54734186, 9.33177098, 6.57049327, 7.05107535, 7.50556246, 6.36403802, 7.46291883, 7.1162672 , 6.24923368, 4.93801144, 7.40670499, 6.14057921, 4.65421066, 2.15817049, 2.51319563, 2.22952704, 1.1132668 ]) E_x1_2=array([ -0.67525083, -14.73375308, 10.00537697, 9.41899041, 9.39459573, 6.83268464, 7.06107471, 7.51420411, 6.36091674, 7.46347728, 7.11537898, 6.24489176, 4.93953833, 7.40593768, 6.13615142, 4.67033145, 2.16482055, 2.51369606, 2.24718589, 1.11451109]) E_x1=array([ 0.2830845 , -14.40136897, 9.91265667, 9.4878906 , 9.36431839, 6.72108178, 7.05610381, 7.50990478, 6.36248018, 7.46319815, 7.11582332, 6.24706814, 4.93877556, 7.4063215 , 6.13837096, 4.66234585, 2.16150825, 2.51344592, 2.23844622, 1.11388939]) E_x1_env=array([ 10.07857251, -1.84498109, 9.52227945, 8.92739375, 8.99821276, 3.35111321, 6.9190869 , 7.49287247, 6.38041062, 7.46033055, 7.12020813, 6.26986155, 4.93037618, 7.4101133 , 6.1591922 , 4.55874442, 2.13987879, 2.51996763, 2.12339348, 1.11000805]) E_x1_rcb=array([ 1.14687049, -14.00960925, 9.67892898, 9.58384551, 9.30572345, 6.43751408, 7.04894948, 7.50291537, 6.36533426, 7.46270022, 7.11661024, 6.25095455, 4.93739069, 7.40699829, 6.14222723, 4.64708355, 2.15589712, 2.51289629, 2.22200985, 1.11281023]) specs10=[E_x1_1,E_x1_2,E_x1_rcb,E_x1_env] markers0=['^','^','o','s'] colors0=['b','b','r','k'] mfcs0=['b','None','orange','brown'] specs11=[E_x1,E_x1_rcb,E_x1_env] markers1=['^','o','s'] colors1=['b','r','k'] mfcs1=['b','orange','brown'] specs12=[X_Fe1,X_Fe1_rcb,X_Fe1_env] markers2=['^','o','s'] colors2=['b','r','k'] mfcs2=['b','orange','brown'] subplot(4,3,1) ax1=ifig.add_subplot(4,3,1) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs10)): ax1.plot(z[2:7],specs10[a][2:7],marker=markers0[a],color=colors0[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors0[a],mfc=mfcs0[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(3.0,11.1) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) #text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14) for a in range(len(element[2:7])): if a%2.0 == 0: text(z[2:7][a]-0.3,11.3,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a]-0.3,11.5,element[2:7][a],fontsize=9,weight='bold',color='k') subplot(4,3,2) ax1=ifig.add_subplot(4,3,2) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs11)): ax1.plot(z[7:13],specs11[a][7:13],marker=markers1[a],color=colors1[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors1[a],mfc=mfcs1[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(4.5,8.1) xlim(11.5,22.5) #ax1.xaxis.set_major_locator(MultipleLocator(0.5)) for a in range(len(element[7:13])): if a%2.0 == 0: text(z[7:13][a]-0.3,8.3,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a]-0.3,8.5,element[7:13][a],fontsize=9,weight='bold',color='k') subplot(4,3,3) ax1=ifig.add_subplot(4,3,3) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs12)): ax1.plot(z[13:len(z)],specs12[a][13:len(z)],marker=markers2[a],color=colors2[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors2[a],mfc=mfcs2[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.5,5.1) xlim(22.5,60.5) for a in range(len(element[13:len(z)])): if a%2.0 == 0: text(z[13:len(z)][a]-0.3,5.2,element[13:len(z)][a],fontsize=9,weight='bold',color='k') else: text(z[13:len(z)][a]-0.3,5.5,element[13:len(z)][a],fontsize=9,weight='bold',color='k') ax1.xaxis.set_major_locator(MultipleLocator(5)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #------------------------------------------------------------------------------------------- X_Fe2_1=array([-14.51, -15.48, 1.55, 1.55, 1.78, 2.84, 1.03, 0.2 , -0.04, 0.06, 0.11, -0.05, 0.63, 0. , 0.15, 1.19, 3.01, 3.03, 4. , 3.91]) X_Fe2_2=array([-15.67, -35.66, 1.68, 1.45, 1.84, 2.96, 1.05, 0.21, -0.05, 0.06, 0.12, -0.06, 0.65, 0. , 0.16, 1.22, 3.03, 3.04, 4.02, 3.92]) X_Fe2_env=array([-1.83, -2.96, 1.01, 0.92, 1.09, -1.16, 0.71, -0.02, -0.03, -0.05, 0.01, -0.04, 0.02, 0. , -0.05, 0.05, 0.03, 0.01, 0.11, -0. ]) X_Fe2_rcb=array([-11.17, -15.04, 1.14, 1.65, 1.65, 2.45, 1.02, 0.16, -0.04, 0.04, 0.1 , -0.05, 0.57, 0. , 0.12, 1.07, 2.93, 2.96, 3.93, 3.82]) X_Fe2=array([-14.78, -15.77, 1.62, 1.5 , 1.81, 2.9 , 1.04, 0.21, -0.04, 0.06, 0.12, -0.06, 0.64, 0. , 0.16, 1.2 , 3.02, 3.04, 4.01, 3.91]) X_S2_1=array([-14.62, -15.59, 1.43, 1.43, 1.67, 2.72, 0.92, 0.08, -0.16, -0.06, 0. , -0.17, 0.52, -0.11, 0.04, 1.08, 2.89, 2.92, 3.89, 3.79]) X_S2_2=array([-15.79, -35.78, 1.56, 1.33, 1.72, 2.84, 0.93, 0.09, -0.16, -0.06, 0. , -0.18, 0.53, -0.12, 0.05, 1.1 , 2.91, 2.93, 3.9 , 3.8 ]) X_S2_rcb=array([-11.27, -15.14, 1.04, 1.55, 1.55, 2.35, 0.92, 0.06, -0.14, -0.06, 0. , -0.15, 0.47, -0.1 , 0.02, 0.97, 2.83, 2.86, 3.83, 3.73]) X_S2=array([-14.9 , -15.89, 1.5 , 1.39, 1.69, 2.79, 0.92, 0.09, -0.16, -0.06, 0. , -0.17, 0.52, -0.12, 0.04, 1.09, 2.9 , 2.92, 3.9 , 3.8 ]) X_S2_env = array([-1.84, -2.97, 1. , 0.91, 1.08, -1.17, 0.7 , -0.03, -0.04, -0.06, 0. , -0.05, 0.01, -0.01, -0.06, 0.04, 0.02, -0. , 0.1 , -0.01]) E_x2_1= array([ 0.55953046, -14.21530996, 9.7946122 , 9.54734186, 9.33177098, 6.57049327, 7.05107535, 7.50556246, 6.36403802, 7.46291883, 7.1162672 , 6.24923368, 4.93801144, 7.40670499, 6.14057921, 4.65421066, 2.15817049, 2.51319563, 2.22952704, 1.1132668 ]) E_x2_2=array([ -0.67525083, -14.73375308, 10.00537697, 9.41899041, 9.39459573, 6.83268464, 7.06107471, 7.51420411, 6.36091674, 7.46347728, 7.11537898, 6.24489176, 4.93953833, 7.40593768, 6.13615142, 4.67033145, 2.16482055, 2.51369606, 2.24718589, 1.11451109]) E_x2=array([ 0.2830845 , -14.40136897, 9.91265667, 9.4878906 , 9.36431839, 6.72108178, 7.05610381, 7.50990478, 6.36248018, 7.46319815, 7.11582332, 6.24706814, 4.93877556, 7.4063215 , 6.13837096, 4.66234585, 2.16150825, 2.51344592, 2.23844622, 1.11388939]) E_x2_env=array([ 10.07857251, -1.84498109, 9.52227945, 8.92739375, 8.99821276, 3.35111321, 6.9190869 , 7.49287247, 6.38041062, 7.46033055, 7.12020813, 6.26986155, 4.93037618, 7.4101133 , 6.1591922 , 4.55874442, 2.13987879, 2.51996763, 2.12339348, 1.11000805]) E_x2_rcb=array([ 1.14687049, -14.00960925, 9.67892898, 9.58384551, 9.30572345, 6.43751408, 7.04894948, 7.50291537, 6.36533426, 7.46270022, 7.11661024, 6.25095455, 4.93739069, 7.40699829, 6.14222723, 4.64708355, 2.15589712, 2.51289629, 2.22200985, 1.11281023]) specs10=[E_x2_1,E_x2_2,E_x2_rcb,E_x2_env] markers0=['^','^','o','s'] colors0=['b','b','r','k'] mfcs0=['b','None','orange','brown'] specs11=[E_x2,E_x2_rcb,E_x2_env] markers1=['^','o','s'] colors1=['b','r','k'] mfcs1=['b','orange','brown'] specs12=[X_Fe2,X_Fe2_rcb,X_Fe2_env] markers2=['^','o','s'] colors2=['b','r','k'] mfcs2=['b','orange','brown'] subplot(4,3,4) ax1=ifig.add_subplot(4,3,4) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs10)): ax1.plot(z[2:7],specs10[a][2:7],marker=markers0[a],color=colors0[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors0[a],mfc=mfcs0[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(3.0,11.1) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) #text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14) subplot(4,3,5) ax1=ifig.add_subplot(4,3,5) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs11)): ax1.plot(z[7:13],specs11[a][7:13],marker=markers1[a],color=colors1[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors1[a],mfc=mfcs1[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(4.5,8.1) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) subplot(4,3,6) ax1=ifig.add_subplot(4,3,6) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs12)): ax1.plot(z[13:len(z)],specs12[a][13:len(z)],marker=markers2[a],color=colors2[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors2[a],mfc=mfcs2[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ##ylim(1e-8,1.3) #ylim(-0.5,5.1) xlim(22.5,60.5) ax1.xaxis.set_major_locator(MultipleLocator(5)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #-------------------------------------------------------------------------------------------- X_Fe3_1=array([ -8.81000000e+00, -1.33700000e+01, 1.20000000e+00, 1.30000000e+00, 1.44000000e+00, 1.22000000e+00, 1.17000000e+00, 2.30000000e-01, 1.00000000e-02, -4.00000000e-02, 1.00000000e-02, -5.00000000e-02, 4.00000000e-02, 0.00000000e+00, 0.00000000e+00, 3.80000000e-01, 5.20000000e-01, 4.50000000e-01, 7.70000000e-01, 6.20000000e-01]) X_Fe3_2=array([ -9.93000000e+00, -1.39900000e+01, 1.48000000e+00, 1.19000000e+00, 1.53000000e+00, 1.56000000e+00, 1.17000000e+00, 2.30000000e-01, 1.00000000e-02, -4.00000000e-02, 1.00000000e-02, -5.00000000e-02, 4.00000000e-02, 0.00000000e+00, 0.00000000e+00, 4.00000000e-01, 5.20000000e-01, 4.60000000e-01, 7.80000000e-01, 6.30000000e-01]) X_Fe3=array([ -9.07000000e+00, -1.35800000e+01, 1.36000000e+00, 1.25000000e+00, 1.49000000e+00, 1.43000000e+00, 1.17000000e+00, 2.30000000e-01, 1.00000000e-02, -4.00000000e-02, 1.00000000e-02, -5.00000000e-02, 4.00000000e-02, 0.00000000e+00, 0.00000000e+00, 3.90000000e-01, 5.20000000e-01, 4.50000000e-01, 7.70000000e-01, 6.30000000e-01]) X_Fe3_rcb= array([ -7.76000000e+00, -1.26800000e+01, 8.70000000e-01, 1.36000000e+00, 1.39000000e+00, 8.70000000e-01, 1.17000000e+00, 2.20000000e-01, 1.00000000e-02, -4.00000000e-02, 1.00000000e-02, -5.00000000e-02, 4.00000000e-02, 0.00000000e+00, 0.00000000e+00, 3.70000000e-01, 5.10000000e-01, 4.40000000e-01, 7.60000000e-01, 6.10000000e-01]) X_Fe3_env=array([-1.8 , -2.95, 0.58, 0.94, 1.13, 0.24, 0.73, 0.14, -0.01, -0.04, 0.01, -0.04, 0.03, 0. , -0.02, 0.23, 0.36, 0.32, 0.59, 0.44]) X_S3_1=array([ -8.82000000e+00, -1.33800000e+01, 1.19000000e+00, 1.29000000e+00, 1.43000000e+00, 1.21000000e+00, 1.16000000e+00, 2.10000000e-01, -0.00000000e+00, -5.00000000e-02, 0.00000000e+00, -6.00000000e-02, 3.00000000e-02, -1.00000000e-02, -1.00000000e-02, 3.70000000e-01, 5.00000000e-01, 4.40000000e-01, 7.60000000e-01, 6.10000000e-01]) X_S3_2=array([ -9.95000000e+00, -1.40000000e+01, 1.47000000e+00, 1.17000000e+00, 1.52000000e+00, 1.55000000e+00, 1.16000000e+00, 2.20000000e-01, -0.00000000e+00, -5.00000000e-02, 0.00000000e+00, -6.00000000e-02, 3.00000000e-02, -1.00000000e-02, -1.00000000e-02, 3.90000000e-01, 5.10000000e-01, 4.40000000e-01, 7.60000000e-01, 6.20000000e-01]) X_S3=array([ -9.09000000e+00, -1.35900000e+01, 1.35000000e+00, 1.24000000e+00, 1.48000000e+00, 1.41000000e+00, 1.16000000e+00, 2.20000000e-01, -0.00000000e+00, -5.00000000e-02, 0.00000000e+00, -6.00000000e-02, 3.00000000e-02, -1.00000000e-02, -1.00000000e-02, 3.80000000e-01, 5.10000000e-01, 4.40000000e-01, 7.60000000e-01, 6.10000000e-01]) X_S3_rcb=array([ -7.77000000e+00, -1.26900000e+01, 8.50000000e-01, 1.35000000e+00, 1.37000000e+00, 8.50000000e-01, 1.15000000e+00, 2.10000000e-01, -0.00000000e+00, -5.00000000e-02, 0.00000000e+00, -6.00000000e-02, 3.00000000e-02, -1.00000000e-02, -1.00000000e-02, 3.60000000e-01, 4.90000000e-01, 4.30000000e-01, 7.40000000e-01, 6.00000000e-01]) X_S3_env=array([-1.82, -2.96, 0.57, 0.93, 1.12, 0.23, 0.72, 0.13, -0.02, -0.06, 0. , -0.05, 0.02, -0.01, -0.03, 0.22, 0.35, 0.31, 0.58, 0.43]) E_x3_1=array([ 3.09623779, -12.26695355, 9.70576926, 9.30178092, 9.34464316, 5.72487362, 7.37192505, 7.72848913, 6.40909928, 7.46506296, 7.113623 , 6.25326152, 4.94306805, 7.40145074, 6.20475582, 4.88618833, 2.61649033, 2.95282509, 2.76874508, 1.723017 ]) E_x3_2=array([ 1.96717743, -12.88534171, 9.98067923, 9.18748432, 9.43017854, 6.0631606 , 7.37545653, 7.73549744, 6.40809698, 7.4654677 , 7.11299703, 6.25055488, 4.94423546, 7.4006262 , 6.20250917, 4.90162977, 2.62224194, 2.95616883, 2.7777599 , 1.72944579]) E_x3=array([ 2.82632997, -12.47429042, 9.86462287, 9.24838186, 9.38951326, 5.92615454, 7.37369438, 7.73200742, 6.40859842, 7.46526538, 7.11331013, 6.25191031, 4.94365215, 7.40103867, 6.20363395, 4.89397767, 2.61937566, 2.95450018, 2.77327588, 1.72624329]) E_x3_env=array([ 10.10222097, -1.84492435, 9.08886021, 8.94806658, 9.03726742, 4.74987052, 6.9346631 , 7.64328254, 6.39763107, 7.46227937, 7.11732289, 6.26510067, 4.93613047, 7.40694452, 6.18443768, 4.73397116, 2.46868911, 2.82606089, 2.59715096, 1.544053 ]) E_x3_rcb=array([ 4.14480282, -11.57752859, 9.3671298 , 9.36176352, 9.28802595, 5.36726756, 7.36794004, 7.7221506 , 6.40923256, 7.46474395, 7.11409161, 6.25515121, 4.94216635, 7.4021256 , 6.20520165, 4.87175041, 2.60860366, 2.94647984, 2.75801127, 1.71367322]) specs10=[E_x3_1,E_x3_2,E_x3_rcb,E_x3_env] markers0=['^','^','o','s'] colors0=['b','b','r','k'] mfcs0=['b','None','orange','brown'] specs11=[E_x3,E_x3_rcb,E_x3_env] markers1=['^','o','s'] colors1=['b','r','k'] mfcs1=['b','orange','brown'] specs12=[X_Fe3,X_Fe3_rcb,X_Fe3_env] markers2=['^','o','s'] colors2=['b','r','k'] mfcs2=['b','orange','brown'] subplot(4,3,7) ax1=ifig.add_subplot(4,3,7) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs10)): ax1.plot(z[2:7],specs10[a][2:7],marker=markers0[a],color=colors0[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors0[a],mfc=mfcs0[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(3.0,11.1) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) #text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14) subplot(4,3,8) ax1=ifig.add_subplot(4,3,8) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs11)): ax1.plot(z[7:13],specs11[a][7:13],marker=markers1[a],color=colors1[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors1[a],mfc=mfcs1[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(4.5,8.1) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) subplot(4,3,9) ax1=ifig.add_subplot(4,3,9) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs12)): ax1.plot(z[13:len(z)],specs12[a][13:len(z)],marker=markers2[a],color=colors2[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors2[a],mfc=mfcs2[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) #ylim(1e-8,1.3) ylim(-0.5,5.1) xlim(22.5,60.5) ax1.xaxis.set_major_locator(MultipleLocator(5)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #----------------------------------------------------------------------------------------------- X_Fe4_1=array([-12.04, -15.44, 1.32, 1.54, 1.88, 2.4 , 0.97, 0.07, 0.17, 0.06, 0.05, -0.07, 0.25, 0. , 0.31, 1.5 , 2.22, 1.91, 0.73, 0.44]) X_Fe4_2=array([-13.1 , -33.29, 1.49, 1.46, 1.92, 2.59, 0.98, 0.08, 0.17, 0.07, 0.06, -0.07, 0.27, 0. , 0.33, 1.54, 2.29, 1.96, 0.77, 0.47]) X_Fe4=array([-12.31, -15.74, 1.41, 1.5 , 1.9 , 2.51, 0.98, 0.08, 0.17, 0.07, 0.05, -0.07, 0.26, 0. , 0.32, 1.52, 2.26, 1.94, 0.75, 0.45]) X_Fe4_rcb=array([ -9.44, -14.2 , 0.83, 1.62, 1.77, 1.95, 0.96, 0.04, 0.15, 0.03, 0.04, -0.06, 0.18, 0. , 0.23, 1.32, 2. , 1.69, 0.59, 0.32]) X_Fe4_env=array([-1.82, -2.96, 0.79, 0.93, 0.88, -1.32, 0.71, -0.02, -0.03, -0.05, 0.01, -0.04, 0.02, 0. , -0.05, 0.05, 0.03, 0.01, 0.11, -0. ]) X_S4_1=array([ -1.21000000e+01, -1.54900000e+01, 1.27000000e+00, 1.49000000e+00, 1.83000000e+00, 2.35000000e+00, 9.20000000e-01, 2.00000000e-02, 1.10000000e-01, 1.00000000e-02, 0.00000000e+00, -1.20000000e-01, 2.00000000e-01, -5.00000000e-02, 2.60000000e-01, 1.45000000e+00, 2.17000000e+00, 1.85000000e+00, 6.70000000e-01, 3.80000000e-01]) X_S4_2=array([ -1.31500000e+01, -3.33400000e+01, 1.43000000e+00, 1.40000000e+00, 1.86000000e+00, 2.54000000e+00, 9.30000000e-01, 3.00000000e-02, 1.10000000e-01, 2.00000000e-02, 0.00000000e+00, -1.30000000e-01, 2.20000000e-01, -6.00000000e-02, 2.70000000e-01, 1.49000000e+00, 2.23000000e+00, 1.91000000e+00, 7.10000000e-01, 4.10000000e-01]) #X_S4_rcb=array([ -9.47, -14.24, 0.79, 1.58, 1.73, 1.91, 0.92, -0. , 0.12, -0. , 0. , -0.1 , 0.15, -0.04, 0.19, 1.28, 1.96, 1.66, 0.55, 0.29]) X_S4_env=array([-1.83, -2.97, 0.78, 0.92, 0.87, -1.33, 0.7 , -0.03, -0.04, -0.06, 0. , -0.05, 0.01, -0.01, -0.06, 0.04, 0.02, -0. , 0.1 , -0.01]) X_S4=array([ -1.23600000e+01, -1.57900000e+01, 1.36000000e+00, 1.45000000e+00, 1.85000000e+00, 2.45000000e+00, 9.20000000e-01, 2.00000000e-02, 1.10000000e-01, 1.00000000e-02, 0.00000000e+00, -1.30000000e-01, 2.10000000e-01, -5.00000000e-02, 2.60000000e-01, 1.47000000e+00, 2.20000000e+00, 1.88000000e+00, 6.90000000e-01, 4.00000000e-01]) E_x4_1=array([ -0.21372421, -14.40685412, 9.75360106, 9.47220989, 9.71281164, 6.83016836, 7.10193456, 7.50150421, 6.49689038, 7.49452748, 7.08293577, 6.15991921, 5.08434531, 7.33102494, 6.43839587, 5.93103672, 4.25498218, 4.33728016, 2.6568567 , 1.4675211 ]) E_x4_2=array([ -1.27369059, -32.26527483, 9.9086854 , 9.38256306, 9.74435794, 7.01541244, 7.10803076, 7.5075543 , 6.4909181 , 7.49649263, 7.07970975, 6.14870194, 5.0979172 , 7.32363307, 6.45198376, 5.96730584, 4.30938934, 4.38761067, 2.69413655, 1.49395713]) E_x4=array([ -0.47848347, -14.70788412, 9.83802924, 9.42969548, 9.72887116, 6.93259318, 7.10499335, 7.50453979, 6.49391451, 7.49551116, 7.08132576, 6.15434679, 5.09118427, 7.32734473, 6.44524296, 5.94954979, 4.2830372 , 4.36317411, 2.67589651, 1.48094023]) E_x4_env=array([ 10.09511353, -1.84498514, 9.30050098, 8.94075602, 8.78859979, 3.19466329, 6.91886874, 7.4918407 , 6.38047321, 7.46032088, 7.12022376, 6.26995464, 4.93035802, 7.41014534, 6.15944809, 4.55898109, 2.13988969, 2.52002832, 2.12279669, 1.11005199]) E_x4_rcb=array([ 2.41860236, -13.14792107, 9.28122401, 9.57181663, 9.62507899, 6.40114149, 7.11671085, 7.49215526, 6.50837799, 7.48809085, 7.09244163, 6.19057689, 5.03832589, 7.35469464, 6.38344971, 5.77231434, 4.05425067, 4.14883479, 2.54277019, 1.37867934]) specs10=[E_x4_1,E_x4_2,E_x4_rcb,E_x4_env] markers0=['^','^','o','s'] colors0=['b','b','r','k'] mfcs0=['b','None','orange','brown'] specs11=[E_x4,E_x4_rcb,E_x4_env] markers1=['^','o','s'] colors1=['b','r','k'] mfcs1=['b','orange','brown'] specs12=[X_Fe4,X_Fe4_rcb,X_Fe4_env] markers2=['^','o','s'] colors2=['b','r','k'] mfcs2=['b','orange','brown'] subplot(4,3,10) ax1=ifig.add_subplot(4,3,10) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs10)): ax1.plot(z[2:7],specs10[a][2:7],marker=markers0[a],color=colors0[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors0[a],mfc=mfcs0[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(3.0,11.1) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) #text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14) for a in range(len(element[2:7])): if a%2.0 == 0: text(z[2:7][a]-0.3,1.0,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a]-0.3,1.3,element[2:7][a],fontsize=9,weight='bold',color='k') subplot(4,3,11) ax1=ifig.add_subplot(4,3,11) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs11)): ax1.plot(z[7:13],specs11[a][7:13],marker=markers1[a],color=colors1[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors1[a],mfc=mfcs1[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(4.5,8.1) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) for a in range(len(element[7:13])): if a%2.0 == 0: text(z[7:13][a]-0.3,3.9,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a]-0.3,3.6,element[7:13][a],fontsize=9,weight='bold',color='k') subplot(4,3,12) ax1=ifig.add_subplot(4,3,12) box = ax1.get_position() ax1.set_position([box.x0, box.y0, box.width * 0.90, box.height]) for a in range(len(specs12)): ax1.plot(z[13:len(z)],specs12[a][13:len(z)],marker=markers2[a],color=colors2[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors2[a],mfc=mfcs2[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.5,5.1) xlim(22.5,60.5) for a in range(len(element[13:len(z)])): if a%2.0 == 0: text(z[13:len(z)][a]-0.3,-2.0,element[13:len(z)][a],fontsize=9,weight='bold',color='k') else: text(z[13:len(z)][a]-0.3,-2.5,element[13:len(z)][a],fontsize=9,weight='bold',color='k') ax1.xaxis.set_major_locator(MultipleLocator(5)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) savefig('X_Fe_mod.eps');savefig('X_Fe_mod.pdf');savefig('X_Fe_mod.png')