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(60) ifig=figure(60,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]) specs1=[X_S1,X_S1_rcb,X_S1_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs11=[X_S1_1,X_S1_2,X_S1_rcb,X_S1_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.5,3.0) 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,3.1,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a]-0.3,3.3,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.2,0.61) 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,0.62,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a]-0.3,0.68,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,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.6,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]) specs2=[X_S2,X_S2_rcb,X_S2_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs21=[X_S2_1,X_S2_2,X_S2_rcb,X_S2_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.5,3.0) #ylim(1e-8,1.3) 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(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.2,0.61) 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(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,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]) specs3=[X_S3,X_S3_rcb,X_S3_env] markers=['^','o','s'] colors=['b','r','k'] mfcs=['b','orange','brown'] specs31=[X_S3_1,X_S3_2,X_S3_rcb,X_S3_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.5,3.0) #ylim(1e-8,1.3) 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(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.2,0.60) 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(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,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]) specs4=[X_S4,X_S4_env] markers=['^','s'] colors=['b','k'] mfcs=['b','brown'] specs41=[X_S4_1,X_S4_2,X_S4_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.5,3.0) xlim(5.5,11.5) ax1.xaxis.set_major_locator(MultipleLocator(1)) for a in range(len(element[2:7])): if a%2.0 == 0: text(z[2:7][a],-1.9,element[2:7][a],fontsize=9,weight='bold',color='k') else: text(z[2:7][a],-2.2,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.2,0.61) 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.28,element[7:13][a],fontsize=9,weight='bold',color='k') else: text(z[7:13][a],-0.32,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,5.1) 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],-1.0,element[13:len(z)][a],fontsize=9,weight='bold',color='k') else: text(z[13:len(z)][a],-1.5,element[13:len(z)][a],fontsize=9,weight='bold',color='k') savefig('X_S_mod.eps');savefig('X_S_mod.pdf');savefig('X_S_mod.png')