File:Zagoskin lake alaska selected pollen paleotemperature.svg

1. 1. temperature

indata_1 = np.genfromtxt(infiili_1, delimiter=';', dtype=None) muoto_1=np.shape(indata_1) rowis_1=muoto_1[0] colis_1=muoto_1[1] data_1=np.zeros((rowis_1, colis_1))

for i in xrange(1,rowis_1):

     for j in xrange(0,colis_1):
           data_1[i-1][j]=np.float(indata_1[i][j])

rowis_1=rowis_1-1

pdata_1=data_1.transpose() pdata_1[0][rowis_1]=pdata_1[0][rowis_1-1] pdata_1[1][rowis_1]=pdata_1[1][rowis_1-1]

1. 1. pollen data

indata_2 = np.genfromtxt(infiili_2, delimiter=';', dtype=None) muoto_2=np.shape(indata_2) rowis_2=muoto_2[0] colis_2=muoto_2[1] data_2=np.zeros((rowis_2, colis_2))

for i in xrange(1,rowis_2):

1. for i in xrange(1,50):
2. print ">",i

     for j in xrange(0,colis_2):
           data_2[i-1][j]=np.float(indata_2[i][j])

1. print indata_2[i][j]

rowis_2=rowis_2-1

pdata_2=data_2.transpose() pdata_2[0][rowis_2]=pdata_2[0][rowis_2-1] pdata_2[1][rowis_2]=pdata_2[1][rowis_2-1]

1. 1. 1. plottingi

fig = plt.figure()

ax2=plt.subplot()

x = pdata_2[0] y = np.row_stack(( pdata_2[1],pdata_2[2], pdata_2[3], pdata_2[4],pdata_2[5], pdata_2[6], pdata_2[7], pdata_2[8]))

percent = y / y.sum(axis=0).astype(float) * 100

1. Cal_yr_BP;Alnus;Betula;Salix;Ericales;Cyperaceae;Poaceae;Artemisia;Sphagnum

1. ax2.stackplot(x,y, colors=['#377EB8','#55BA87','#008000','#FFFF00','#FFFF80','#800000'])

ax2.stackplot(x,percent, colors=['#377EB8','#55BA87','#008000','#00F000', '#FFFF00','#FFFF80','#800000', '#8000F0'],

edgecolors=['#377EB8'], linewidth=0

)

ax2.set_ylabel('Pollen percent \n from selected datasets', color='#004000', fontsize=15) for tl in ax2.get_yticklabels():

   tl.set_color('#004000')

ax2.set_ylim([0, 100]) ax2.set_xlim([0, 25000])

1. Shrink current axis's height by 10% on the bottom

box = ax2.get_position()

ax2.set_position([box.x0, box.y0, box.width * 0.8, box.height])

ax2.legend( [mpatches.Patch(color='#377EB8'),

            mpatches.Patch(color='#55BA87'), 
            mpatches.Patch(color='#008000'),
            mpatches.Patch(color='#00F000'), 
            mpatches.Patch(color='#FFFF00'),
            mpatches.Patch(color='#FFFF80'),
            mpatches.Patch(color='#800000'),
            mpatches.Patch(color='#8000F0')],  ['Alnus','Betula','Salix','Ericaceae', 'Cyperaceae','Poaceae', 'Artemisia', 'Sphagnum'],

loc='center left', bbox_to_anchor=(1.1, 0.75) )

1. 1. Cyperaceae #FFFF00
   2. Poaceae #FFFF80
   3. Artemisia #800000
   4. Brassicaceae ##008000
   5. Salix #008000

for label in (ax2.get_xticklabels() + ax2.get_yticklabels()):

   label.set_fontname('Arial')
   label.set_fontsize(15)

ax1 = ax2.twinx()

ax1.set_xlim([0, 25000])

degree_sign= u'\N{DEGREE SIGN}'

tempilabel1='TJuly ' tempilabel2=tempilabel1.encode('utf8')+degree_sign+'C';

ax1.plot(pdata_1[0], pdata_1[1], label=tempilabel2,color='#200000', linewidth=5)

box = ax1.get_position()

ax1.set_position([box.x0, box.y0, box.width * 0.8, box.height])

1. Set the tick labels font

for label in (ax1.get_xticklabels() + ax1.get_yticklabels()):

   label.set_fontname('Arial')
   label.set_fontsize(15)

ax1.set_xlabel('Time (cal BP)', fontsize=15)

1. ax1.get_xticklabels().set_fontsize(20)

ax1.set_ylabel(tempilabel2, color='#a00000', fontsize=15) for tl in ax1.get_yticklabels():

   tl.set_color('#400000')

kaptioni="Zagoskin Lake, Alaska \n paleotemperature and pollen diagram"

plt.title(kaptioni, fontsize=18)

ax1.legend(loc='upper left')

start, end = ax1.get_xlim() ax1.xaxis.set_ticks(np.arange(start, end, 2500))

ax1.xaxis.set_tick_params(width=4, length=10) ax1.yaxis.set_tick_params(width=4) ax2.xaxis.set_tick_params(width=4) ax2.yaxis.set_tick_params(width=4)

plt.gca().invert_xaxis()

plt.show()

fig.set_size_inches(16, 6) fig.savefig('out.svg', dpi=80)

1. fig.savefig('test.svg')