from numpy import * import matplotlib from matplotlib.pylab import * import ppne 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=(10,12)) 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]) X1 = [-11.72,-15.6 ,1.31 ,1.39 ,1.36 ,2.12 ,0.76 ,-0.09 ,-0.14 ,-0.14 ,-0.08 ,-0.15 ,-0.06 ,-0.09 ,-0.16 ,0.06 ,-0.04 ,-0.09 ,0.14 ,-0.09] X1_1 = [-11.44,-15.42 ,1.19 ,1.45 ,1.33 ,1.97 ,0.75 ,-0.09 ,-0.14 ,-0.14 ,-0.08 ,-0.15 ,-0.06 ,-0.09 ,-0.16 ,0.05 ,-0.04 ,-0.09 ,0.13 ,-0.09] X1_2 = [-12.68,-15.93 ,1.41 ,1.32 ,1.39 ,2.23 ,0.76 ,-0.09 ,-0.14 ,-0.14 ,-0.08 ,-0.16 ,-0.06 ,-0.09 ,-0.16 ,0.07 ,-0.04 ,-0.09 ,0.15 ,-0.09] X1_env = [-1.92,-3.04 ,0.92 ,0.83 ,1. ,-1.25 ,0.62,-0.11,-0.12,-0.14,-0.08,-0.13,-0.07,-0.09,-0.14,-0.04,-0.06,-0.08 ,0.02,-0.09] X1_rcb = array([-10.52, -15.02, 1.01, 1.5 , 1.29, 1.76, 0.75, -0.1 , -0.13, -0.14, -0.08, -0.15, -0.06, -0.09, -0.16, 0.04, -0.05, -0.09, 0.12, -0.09]) specs1=[X1,X1_rcb,X1_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs11=[X1_1,X1_2,X1_rcb,X1_env] markers1=['^','^','o','s'] colors1=['b','b','r','k'] mfcs1=['b','None','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(specs11)): ax1.plot(z[2:7],specs11[a][2:7],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(-1.6,2.6) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #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],2.8,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a],3.0,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(specs1)): ax1.plot(z[7:13],specs1[a][7:13],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.31,0.2) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) ax1.yaxis.set_major_locator(MultipleLocator(0.1)) for a in range(len(element[7:13])): if a%2.0 == 0: text(z[7:13][a],0.22,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a],0.25,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(specs1)): ax1.plot(z[13:len(z)],specs1[a][13:len(z)],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ##ylim(1e-8,1.3) ylim(-0.5,2.5) 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],2.6,element[13:len(z)][a],fontsize=9,weight='bold',color='k') else: text(z[13:len(z)][a],2.9,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. ]) XFe2_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]) X2 = array([-12.52, -15.82, 1.38, 1.33, 1.49, 2.25, 0.76, -0.07, -0.21, -0.13, -0.11, -0.23, 0.09, -0.17, 0.21, 1.26, 1.97, 1.69, 1.02, 0.63]) X2_1 =array([-12.24, -15.51, 1.29, 1.38, 1.46, 2.15, 0.76, -0.07, -0.21, -0.13, -0.11, -0.22, 0.09, -0.17, 0.21, 1.24, 1.95, 1.68, 1.01, 0.62]) X2_2 = array([-13.54, -21.77, 1.46, 1.26, 1.53, 2.33, 0.77, -0.06, -0.22, -0.13, -0.11, -0.23, 0.1 , -0.17, 0.21, 1.27, 1.98, 1.7 , 1.04, 0.64]) X2_env = array([-1.92, -3.04, 0.92, 0.83, 1. , -1.25, 0.62, -0.11, -0.12, -0.14, -0.08, -0.13, -0.07, -0.09, -0.14, -0.04, -0.06, -0.08, 0.02, -0.09]) X2_rcb=array([-10.09, -14.67, 0.96, 1.47, 1.35, 1.84, 0.76, -0.08, -0.2 , -0.13, -0.11, -0.21, 0.07, -0.16, 0.19, 1.19, 1.89, 1.63, 0.95, 0.57]) specs2=[X2,X2_rcb,X2_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs21=[X2_1,X2_2,X2_rcb,X2_env] markers1=['^','^','o','s'] colors1=['b','b','r','k'] mfcs1=['b','None','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(specs21)): ax1.plot(z[2:7],specs21[a][2:7],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(-1.6,2.6) ##ylim(1e-8,1.3) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #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(specs2)): ax1.plot(z[7:13],specs2[a][7:13],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.31,0.2) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) ax1.yaxis.set_major_locator(MultipleLocator(0.1)) 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(specs2)): ax1.plot(z[13:len(z)],specs2[a][13:len(z)],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ##ylim(1e-8,1.3) ylim(-0.5,2.5) 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]) X3 = array([ -9.17, -13.67, 1.26, 1.15, 1.39, 1.33, 1.07, 0.13, -0.09, -0.13, -0.09, -0.15, -0.06, -0.1 , -0.1 , 0.29, 0.42, 0.35, 0.67, 0.53]) X3_1 = array([ -8.9 , -13.47, 1.11, 1.2 , 1.34, 1.12, 1.07, 0.13, -0.09, -0.13, -0.09, -0.15, -0.06, -0.1 , -0.1 , 0.29, 0.42, 0.35, 0.67, 0.52]) X3_2 =array([-10.03, -14.09, 1.38, 1.09, 1.43, 1.46, 1.08, 0.14, -0.09, -0.13, -0.09, -0.15, -0.06, -0.1 , -0.1 , 0.3 , 0.42, 0.36, 0.68, 0.53]) X3_env =array([-1.9 , -3.04, 0.49, 0.85, 1.04, 0.15, 0.63, 0.04, -0.1 , -0.14, -0.08, -0.13, -0.06, -0.09, -0.12, 0.13, 0.27, 0.23, 0.5 , 0.34]) X3_rcb =array([ -7.86, -12.78, 0.77, 1.26, 1.29, 0.77, 1.07, 0.12, -0.09, -0.14, -0.09, -0.14, -0.06, -0.1 , -0.09, 0.27, 0.41, 0.35, 0.66, 0.51]) specs3=[X3,X3_rcb,X3_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs31=[X3_1,X3_2,X3_rcb,X3_env] markers1=['^','^','o','s'] colors1=['b','b','r','k'] mfcs1=['b','None','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(specs31)): ax1.plot(z[2:7],specs31[a][2:7],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(-1.6,2.6) ##ylim(1e-8,1.3) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #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(specs3)): ax1.plot(z[7:13],specs3[a][7:13],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.31,0.2) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) ax1.yaxis.set_major_locator(MultipleLocator(0.1)) 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(specs3)): ax1.plot(z[13:len(z)],specs3[a][13:len(z)],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.5,2.5) 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]) X4 = array([ -1.24800000e+01, -1.59100000e+01, 1.24000000e+00, 1.33000000e+00, 1.73000000e+00, 2.33000000e+00, 8.00000000e-01, -1.00000000e-01, -1.00000000e-02, -1.00000000e-01, -1.20000000e-01, -2.50000000e-01, 9.00000000e-02, -1.70000000e-01, 1.50000000e-01, 1.35000000e+00, 2.08000000e+00, 1.76000000e+00, 5.80000000e-01, 2.80000000e-01]) X4_1 = array([-12.21, -15.61, 1.15, 1.37, 1.71, 2.23, 0.8 , -0.1 , -0. , -0.11, -0.12, -0.24, 0.08, -0.17, 0.14, 1.33, 2.05, 1.74, 0.56, 0.27]) X4_2 =array([ -1.32700000e+01, -3.34700000e+01, 1.31000000e+00, 1.28000000e+00, 1.74000000e+00, 2.42000000e+00, 8.10000000e-01, -9.00000000e-02, -1.00000000e-02, -1.00000000e-01, -1.20000000e-01, -2.50000000e-01, 1.00000000e-01, -1.80000000e-01, 1.50000000e-01, 1.37000000e+00, 2.11000000e+00, 1.79000000e+00, 5.90000000e-01, 2.90000000e-01]) X4_env= array([-1.9 , -3.04, 0.7 , 0.84, 0.79, -1.41, 0.62, -0.11, -0.12, -0.14, -0.08, -0.13, -0.07, -0.09, -0.14, -0.04, -0.06, -0.08, 0.02, -0.09]) X4_rcb= array([ -9.58000000e+00, -1.43500000e+01, 6.80000000e-01, 1.47000000e+00, 1.63000000e+00, 1.80000000e+00, 8.20000000e-01, -1.10000000e-01, 1.00000000e-02, -1.10000000e-01, -1.10000000e-01, -2.10000000e-01, 4.00000000e-02, -1.50000000e-01, 8.00000000e-02, 1.17000000e+00, 1.85000000e+00, 1.55000000e+00, 4.40000000e-01, 1.80000000e-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]) specs4=[X4,X4_env] markers=['^','s'] colors=['b','k'] mfcs=['b','brown'] specs41=[X4_1,X4_2,X4_env] markers1=['^','^','s'] colors1=['b','b','k'] mfcs1=['b','None','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(specs41)): ax1.plot(z[2:7],specs41[a][2:7],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(-1.6,2.6) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) for a in range(len(element[2:7])): if a%2.0 == 0: text(z[2:7][a],-2.1,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a],-2.4,element[2:7][a],fontsize=9,weight='bold',color='k') #text(0.35, 1.0, 'a', ha='center', va='center', rotation='horizontal',color='k',fontsize=14) 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(specs4)): ax1.plot(z[7:13],specs4[a][7:13],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.31,0.2) xlim(11.5,22.5) ax1.xaxis.set_major_locator(MultipleLocator(2)) ax1.yaxis.set_major_locator(MultipleLocator(0.1)) for a in range(len(element[7:13])): if a%2.0 == 0: text(z[7:13][a],-0.38,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a],-0.41,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(specs4)): ax1.plot(z[13:len(z)],specs4[a][13:len(z)],marker=markers[a],color=colors[a],linestyle='',markersize=5.5,alpha=0.7,mec=colors[a],mfc=mfcs[a],mew=1.4) setp(ax1.get_xticklabels(), visible=False) ylim(-0.5,2.5) xlim(22.5,60.5) ax1.xaxis.set_major_locator(MultipleLocator(5)) ax1.yaxis.set_major_locator(MultipleLocator(0.5)) #for a in range(len(element[2:])): #text(z[2:][a],-0.8,element[2:][a],fontsize=11,weight='bold',color='k') for a in range(len(element[13:len(z)])): if a%2.0 == 0: text(z[13:len(z)][a],-0.85,element[13:len(z)][a],fontsize=9,weight='bold',color='k') else: text(z[13:len(z)][a],-1.15,element[13:len(z)][a],fontsize=9,weight='bold',color='k') savefig('X_mod.eps');savefig('X_mod.pdf');savefig('X_mod.png') X_master=[X1_rcb,X1_2,X2_rcb,X2_2,X3_rcb,X3_2,X4_rcb,X4_2] X_master_label=['X1_rcb','X1_2','X2_rcb','X2_2','X3_rcb','X3_2','X4_rcb','X4_2'] X_min=[] X_max=[] ele_temp=[] X_master_min=[] X_master_max=[] fe_mod_min=[] fe_mod_max=[] for i in range(len(X1)): for j in range(len(X_master)): ele_temp.append(X_master[j][i]) X_min.append(min(ele_temp)) X_max.append(max(ele_temp)) X_master_min.append(X_master_label[ele_temp.index(min(ele_temp))]) X_master_max.append(X_master_label[ele_temp.index(max(ele_temp))]) fe_mod_min.append(X_master[ele_temp.index(min(ele_temp))][13]) fe_mod_max.append(X_master[ele_temp.index(max(ele_temp))][13]) ele_temp=[]