Libraries


In [1]:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import itertools
import seaborn as sns
from mpl_toolkits.mplot3d import Axes3D

In [2]:
%matplotlib notebook

Dataset


In [3]:
new_probability  = np.loadtxt('./Results/fit_results_3d_grid.csv', delimiter=',', dtype=str)
results          = np.loadtxt('./Results/model_prob_3d_grid.csv', delimiter=',', dtype=str)

In [4]:
my_probability = {}
for i in range(len(new_probability[0, :])):                             # Converting numpy array into dictionary
    my_probability[new_probability[0, i]] = np.array(new_probability[0 + 1:, i], dtype=str)

In [5]:
param = my_probability['parameter']
for i in range(param.size):
    print param[i]


beta[0]
beta[1]
beta[2]
beta[3]
eta[0]
eta[1]
eta[2]
eta[3]
eta[4]
eta[5]
eta[6]
eta[7]
eta[8]
eta[9]
eta[10]
eta[11]
eta[12]
eta[13]
eta[14]
eta[15]
eta[16]
eta[17]
eta[18]
eta[19]
eta[20]
eta[21]
eta[22]
eta[23]
eta[24]
eta[25]
eta[26]
eta[27]
eta[28]
eta[29]
eta[30]
eta[31]
eta[32]
eta[33]
eta[34]
eta[35]
eta[36]
eta[37]
eta[38]
eta[39]
eta[40]
eta[41]
eta[42]
eta[43]
eta[44]
eta[45]
eta[46]
eta[47]
eta[48]
eta[49]
eta[50]
eta[51]
eta[52]
eta[53]
eta[54]
eta[55]
eta[56]
eta[57]
eta[58]
eta[59]
eta[60]
eta[61]
eta[62]
eta[63]
eta[64]
eta[65]
eta[66]
eta[67]
eta[68]
eta[69]
eta[70]
eta[71]
eta[72]
eta[73]
eta[74]
eta[75]
eta[76]
eta[77]
eta[78]
eta[79]
eta[80]
eta[81]
eta[82]
eta[83]
eta[84]
eta[85]
eta[86]
eta[87]
eta[88]
eta[89]
eta[90]
eta[91]
eta[92]
eta[93]
eta[94]
eta[95]
eta[96]
eta[97]
eta[98]
eta[99]
eta[100]
eta[101]
eta[102]
eta[103]
eta[104]
eta[105]
eta[106]
eta[107]
eta[108]
eta[109]
eta[110]
eta[111]
eta[112]
eta[113]
eta[114]
eta[115]
eta[116]
eta[117]
eta[118]
eta[119]
eta[120]
eta[121]
eta[122]
eta[123]
eta[124]
eta[125]
eta[126]
eta[127]
eta[128]
eta[129]
eta[130]
eta[131]
eta[132]
eta[133]
eta[134]
eta[135]
eta[136]
eta[137]
eta[138]
eta[139]
eta[140]
eta[141]
eta[142]
eta[143]
eta[144]
eta[145]
eta[146]
eta[147]
eta[148]
eta[149]
eta[150]
eta[151]
eta[152]
eta[153]
eta[154]
eta[155]
eta[156]
eta[157]
eta[158]
eta[159]
eta[160]
eta[161]
eta[162]
eta[163]
eta[164]
eta[165]
eta[166]
eta[167]
eta[168]
eta[169]
eta[170]
eta[171]
eta[172]
eta[173]
eta[174]
eta[175]
eta[176]
eta[177]
eta[178]
eta[179]
eta[180]
eta[181]
eta[182]
eta[183]
eta[184]
eta[185]
eta[186]
eta[187]
eta[188]
eta[189]
eta[190]
eta[191]
eta[192]
eta[193]
eta[194]
eta[195]
eta[196]
eta[197]
eta[198]
eta[199]
eta[200]
eta[201]
eta[202]
eta[203]
eta[204]
eta[205]
eta[206]
eta[207]
eta[208]
eta[209]
eta[210]
eta[211]
eta[212]
eta[213]
eta[214]
eta[215]
eta[216]
eta[217]
eta[218]
eta[219]
eta[220]
eta[221]
eta[222]
eta[223]
eta[224]
eta[225]
eta[226]
eta[227]
eta[228]
eta[229]
eta[230]
eta[231]
eta[232]
eta[233]
eta[234]
eta[235]
eta[236]
eta[237]
eta[238]
eta[239]
eta[240]
eta[241]
eta[242]
eta[243]
eta[244]
eta[245]
eta[246]
eta[247]
eta[248]
eta[249]
eta[250]
eta[251]
eta[252]
eta[253]
eta[254]
eta[255]
eta[256]
eta[257]
eta[258]
eta[259]
eta[260]
eta[261]
eta[262]
eta[263]
eta[264]
eta[265]
eta[266]
eta[267]
eta[268]
eta[269]
eta[270]
eta[271]
eta[272]
eta[273]
eta[274]
eta[275]
eta[276]
eta[277]
eta[278]
eta[279]
eta[280]
eta[281]
eta[282]
eta[283]
eta[284]
eta[285]
eta[286]
eta[287]
eta[288]
eta[289]
eta[290]
eta[291]
eta[292]
eta[293]
eta[294]
eta[295]
eta[296]
eta[297]
eta[298]
eta[299]
eta[300]
eta[301]
eta[302]
eta[303]
eta[304]
eta[305]
eta[306]
eta[307]
eta[308]
eta[309]
eta[310]
eta[311]
eta[312]
eta[313]
eta[314]
eta[315]
eta[316]
eta[317]
eta[318]
eta[319]
eta[320]
eta[321]
eta[322]
eta[323]
eta[324]
eta[325]
eta[326]
eta[327]
eta[328]
eta[329]
eta[330]
eta[331]
eta[332]
eta[333]
eta[334]
eta[335]
eta[336]
eta[337]
eta[338]
eta[339]
eta[340]
eta[341]
eta[342]
eta[343]
eta[344]
eta[345]
eta[346]
eta[347]
eta[348]
eta[349]
eta[350]
eta[351]
eta[352]
eta[353]
eta[354]
eta[355]
eta[356]
eta[357]
eta[358]
eta[359]
eta[360]
eta[361]
eta[362]
eta[363]
eta[364]
eta[365]
eta[366]
eta[367]
eta[368]
eta[369]
eta[370]
eta[371]
eta[372]
eta[373]
eta[374]
eta[375]
eta[376]
eta[377]
eta[378]
eta[379]
eta[380]
eta[381]
eta[382]
eta[383]
eta[384]
eta[385]
eta[386]
eta[387]
eta[388]
eta[389]
eta[390]
eta[391]
eta[392]
eta[393]
eta[394]
eta[395]
eta[396]
eta[397]
eta[398]
eta[399]
eta[400]
eta[401]
eta[402]
eta[403]
eta[404]
eta[405]
eta[406]
eta[407]
eta[408]
eta[409]
eta[410]
eta[411]
eta[412]
eta[413]
eta[414]
eta[415]
eta[416]
eta[417]
eta[418]
eta[419]
eta[420]
eta[421]
eta[422]
eta[423]
eta[424]
eta[425]
eta[426]
eta[427]
eta[428]
eta[429]
eta[430]
eta[431]
eta[432]
eta[433]
eta[434]
eta[435]
eta[436]
eta[437]
eta[438]
eta[439]
eta[440]
eta[441]
eta[442]
eta[443]
eta[444]
eta[445]
eta[446]
eta[447]
eta[448]
eta[449]
eta[450]
eta[451]
eta[452]
eta[453]
eta[454]
eta[455]
eta[456]
eta[457]
eta[458]
eta[459]
eta[460]
eta[461]
eta[462]
eta[463]
eta[464]
eta[465]
eta[466]
eta[467]
eta[468]
eta[469]
eta[470]
eta[471]
eta[472]
eta[473]
eta[474]
eta[475]
eta[476]
eta[477]
eta[478]
eta[479]
eta[480]
eta[481]
eta[482]
eta[483]
eta[484]
eta[485]
eta[486]
eta[487]
eta[488]
eta[489]
eta[490]
eta[491]
eta[492]
eta[493]
eta[494]
eta[495]
eta[496]
eta[497]
eta[498]
eta[499]
eta[500]
eta[501]
eta[502]
eta[503]
etanew[0]
etanew[1]
etanew[2]
etanew[3]
etanew[4]
etanew[5]
etanew[6]
etanew[7]
etanew[8]
etanew[9]
etanew[10]
etanew[11]
etanew[12]
etanew[13]
etanew[14]
etanew[15]
etanew[16]
etanew[17]
etanew[18]
etanew[19]
etanew[20]
etanew[21]
etanew[22]
etanew[23]
etanew[24]
etanew[25]
etanew[26]
etanew[27]
etanew[28]
etanew[29]
etanew[30]
etanew[31]
etanew[32]
etanew[33]
etanew[34]
etanew[35]
etanew[36]
etanew[37]
etanew[38]
etanew[39]
etanew[40]
etanew[41]
etanew[42]
etanew[43]
etanew[44]
etanew[45]
etanew[46]
etanew[47]
etanew[48]
etanew[49]
etanew[50]
etanew[51]
etanew[52]
etanew[53]
etanew[54]
etanew[55]
etanew[56]
etanew[57]
etanew[58]
etanew[59]
etanew[60]
etanew[61]
etanew[62]
etanew[63]
etanew[64]
etanew[65]
etanew[66]
etanew[67]
etanew[68]
etanew[69]
etanew[70]
etanew[71]
etanew[72]
etanew[73]
etanew[74]
etanew[75]
etanew[76]
etanew[77]
etanew[78]
etanew[79]
etanew[80]
etanew[81]
etanew[82]
etanew[83]
etanew[84]
etanew[85]
etanew[86]
etanew[87]
etanew[88]
etanew[89]
etanew[90]
etanew[91]
etanew[92]
etanew[93]
etanew[94]
etanew[95]
etanew[96]
etanew[97]
etanew[98]
etanew[99]
etanew[100]
etanew[101]
etanew[102]
etanew[103]
etanew[104]
etanew[105]
etanew[106]
etanew[107]
etanew[108]
etanew[109]
etanew[110]
etanew[111]
etanew[112]
etanew[113]
etanew[114]
etanew[115]
etanew[116]
etanew[117]
etanew[118]
etanew[119]
etanew[120]
etanew[121]
etanew[122]
etanew[123]
etanew[124]
etanew[125]
etanew[126]
etanew[127]
etanew[128]
etanew[129]
etanew[130]
etanew[131]
etanew[132]
etanew[133]
etanew[134]
etanew[135]
etanew[136]
etanew[137]
etanew[138]
etanew[139]
etanew[140]
etanew[141]
etanew[142]
etanew[143]
etanew[144]
etanew[145]
etanew[146]
etanew[147]
etanew[148]
etanew[149]
etanew[150]
etanew[151]
etanew[152]
etanew[153]
etanew[154]
etanew[155]
etanew[156]
etanew[157]
etanew[158]
etanew[159]
etanew[160]
etanew[161]
etanew[162]
etanew[163]
etanew[164]
etanew[165]
etanew[166]
etanew[167]
etanew[168]
etanew[169]
etanew[170]
etanew[171]
etanew[172]
etanew[173]
etanew[174]
etanew[175]
etanew[176]
etanew[177]
etanew[178]
etanew[179]
etanew[180]
etanew[181]
etanew[182]
etanew[183]
etanew[184]
etanew[185]
etanew[186]
etanew[187]
etanew[188]
etanew[189]
etanew[190]
etanew[191]
etanew[192]
etanew[193]
etanew[194]
etanew[195]
etanew[196]
etanew[197]
etanew[198]
etanew[199]
etanew[200]
etanew[201]
etanew[202]
etanew[203]
etanew[204]
etanew[205]
etanew[206]
etanew[207]
etanew[208]
etanew[209]
etanew[210]
etanew[211]
etanew[212]
etanew[213]
etanew[214]
etanew[215]
etanew[216]
etanew[217]
etanew[218]
etanew[219]
etanew[220]
etanew[221]
etanew[222]
etanew[223]
etanew[224]
etanew[225]
etanew[226]
etanew[227]
etanew[228]
etanew[229]
etanew[230]
etanew[231]
etanew[232]
etanew[233]
etanew[234]
etanew[235]
etanew[236]
etanew[237]
etanew[238]
etanew[239]
etanew[240]
etanew[241]
etanew[242]
etanew[243]
etanew[244]
etanew[245]
etanew[246]
etanew[247]
etanew[248]
etanew[249]
etanew[250]
etanew[251]
etanew[252]
etanew[253]
etanew[254]
etanew[255]
etanew[256]
etanew[257]
etanew[258]
etanew[259]
etanew[260]
etanew[261]
etanew[262]
etanew[263]
etanew[264]
etanew[265]
etanew[266]
etanew[267]
etanew[268]
etanew[269]
etanew[270]
etanew[271]
etanew[272]
etanew[273]
etanew[274]
etanew[275]
etanew[276]
etanew[277]
etanew[278]
etanew[279]
etanew[280]
etanew[281]
etanew[282]
etanew[283]
etanew[284]
etanew[285]
etanew[286]
etanew[287]
etanew[288]
etanew[289]
etanew[290]
etanew[291]
etanew[292]
etanew[293]
etanew[294]
etanew[295]
etanew[296]
etanew[297]
etanew[298]
etanew[299]
etanew[300]
etanew[301]
etanew[302]
etanew[303]
etanew[304]
etanew[305]
etanew[306]
etanew[307]
etanew[308]
etanew[309]
etanew[310]
etanew[311]
etanew[312]
etanew[313]
etanew[314]
etanew[315]
etanew[316]
etanew[317]
etanew[318]
etanew[319]
etanew[320]
etanew[321]
etanew[322]
etanew[323]
etanew[324]
etanew[325]
etanew[326]
etanew[327]
etanew[328]
etanew[329]
etanew[330]
etanew[331]
etanew[332]
etanew[333]
etanew[334]
etanew[335]
etanew[336]
etanew[337]
etanew[338]
etanew[339]
etanew[340]
etanew[341]
etanew[342]
etanew[343]
etanew[344]
etanew[345]
etanew[346]
etanew[347]
etanew[348]
etanew[349]
etanew[350]
etanew[351]
etanew[352]
etanew[353]
etanew[354]
etanew[355]
etanew[356]
etanew[357]
etanew[358]
etanew[359]
etanew[360]
etanew[361]
etanew[362]
etanew[363]
etanew[364]
etanew[365]
etanew[366]
etanew[367]
etanew[368]
etanew[369]
etanew[370]
etanew[371]
etanew[372]
etanew[373]
etanew[374]
etanew[375]
etanew[376]
etanew[377]
etanew[378]
etanew[379]
etanew[380]
etanew[381]
etanew[382]
etanew[383]
etanew[384]
etanew[385]
etanew[386]
etanew[387]
etanew[388]
etanew[389]
etanew[390]
etanew[391]
etanew[392]
etanew[393]
etanew[394]
etanew[395]
etanew[396]
etanew[397]
etanew[398]
etanew[399]
etanew[400]
etanew[401]
etanew[402]
etanew[403]
etanew[404]
etanew[405]
etanew[406]
etanew[407]
etanew[408]
etanew[409]
etanew[410]
etanew[411]
etanew[412]
etanew[413]
etanew[414]
etanew[415]
etanew[416]
etanew[417]
etanew[418]
etanew[419]
etanew[420]
etanew[421]
etanew[422]
etanew[423]
etanew[424]
etanew[425]
etanew[426]
etanew[427]
etanew[428]
etanew[429]
etanew[430]
etanew[431]
etanew[432]
etanew[433]
etanew[434]
etanew[435]
etanew[436]
etanew[437]
etanew[438]
etanew[439]
etanew[440]
etanew[441]
etanew[442]
etanew[443]
etanew[444]
etanew[445]
etanew[446]
etanew[447]
etanew[448]
etanew[449]
etanew[450]
etanew[451]
etanew[452]
etanew[453]
etanew[454]
etanew[455]
etanew[456]
etanew[457]
etanew[458]
etanew[459]
etanew[460]
etanew[461]
etanew[462]
etanew[463]
etanew[464]
etanew[465]
etanew[466]
etanew[467]
etanew[468]
etanew[469]
etanew[470]
etanew[471]
etanew[472]
etanew[473]
etanew[474]
etanew[475]
etanew[476]
etanew[477]
etanew[478]
etanew[479]
etanew[480]
etanew[481]
etanew[482]
etanew[483]
etanew[484]
etanew[485]
etanew[486]
etanew[487]
etanew[488]
etanew[489]
etanew[490]
etanew[491]
etanew[492]
etanew[493]
etanew[494]
etanew[495]
etanew[496]
etanew[497]
etanew[498]
etanew[499]
etanew[500]
etanew[501]
etanew[502]
etanew[503]
etanew[504]
etanew[505]
etanew[506]
etanew[507]
etanew[508]
etanew[509]
etanew[510]
etanew[511]
etanew[512]
etanew[513]
etanew[514]
etanew[515]
etanew[516]
etanew[517]
etanew[518]
etanew[519]
etanew[520]
etanew[521]
etanew[522]
etanew[523]
etanew[524]
etanew[525]
etanew[526]
etanew[527]
etanew[528]
etanew[529]
etanew[530]
etanew[531]
etanew[532]
etanew[533]
etanew[534]
etanew[535]
etanew[536]
etanew[537]
etanew[538]
etanew[539]
etanew[540]
etanew[541]
etanew[542]
etanew[543]
etanew[544]
etanew[545]
etanew[546]
etanew[547]
etanew[548]
etanew[549]
etanew[550]
etanew[551]
etanew[552]
etanew[553]
etanew[554]
etanew[555]
etanew[556]
etanew[557]
etanew[558]
etanew[559]
etanew[560]
etanew[561]
etanew[562]
etanew[563]
etanew[564]
etanew[565]
etanew[566]
etanew[567]
etanew[568]
etanew[569]
etanew[570]
etanew[571]
etanew[572]
etanew[573]
etanew[574]
etanew[575]
etanew[576]
etanew[577]
etanew[578]
etanew[579]
etanew[580]
etanew[581]
etanew[582]
etanew[583]
etanew[584]
etanew[585]
etanew[586]
etanew[587]
etanew[588]
etanew[589]
etanew[590]
etanew[591]
etanew[592]
etanew[593]
etanew[594]
etanew[595]
etanew[596]
etanew[597]
etanew[598]
etanew[599]
etanew[600]
etanew[601]
etanew[602]
etanew[603]
etanew[604]
etanew[605]
etanew[606]
etanew[607]
etanew[608]
etanew[609]
etanew[610]
etanew[611]
etanew[612]
etanew[613]
etanew[614]
etanew[615]
etanew[616]
etanew[617]
etanew[618]
etanew[619]
etanew[620]
etanew[621]
etanew[622]
etanew[623]
etanew[624]
etanew[625]
etanew[626]
etanew[627]
etanew[628]
etanew[629]
etanew[630]
etanew[631]
etanew[632]
etanew[633]
etanew[634]
etanew[635]
etanew[636]
etanew[637]
etanew[638]
etanew[639]
etanew[640]
etanew[641]
etanew[642]
etanew[643]
etanew[644]
etanew[645]
etanew[646]
etanew[647]
etanew[648]
etanew[649]
etanew[650]
etanew[651]
etanew[652]
etanew[653]
etanew[654]
etanew[655]
etanew[656]
etanew[657]
etanew[658]
etanew[659]
etanew[660]
etanew[661]
etanew[662]
etanew[663]
etanew[664]
etanew[665]
etanew[666]
etanew[667]
etanew[668]
etanew[669]
etanew[670]
etanew[671]
etanew[672]
etanew[673]
etanew[674]
etanew[675]
etanew[676]
etanew[677]
etanew[678]
etanew[679]
etanew[680]
etanew[681]
etanew[682]
etanew[683]
etanew[684]
etanew[685]
etanew[686]
etanew[687]
etanew[688]
etanew[689]
etanew[690]
etanew[691]
etanew[692]
etanew[693]
etanew[694]
etanew[695]
etanew[696]
etanew[697]
etanew[698]
etanew[699]
etanew[700]
etanew[701]
etanew[702]
etanew[703]
etanew[704]
etanew[705]
etanew[706]
etanew[707]
etanew[708]
etanew[709]
etanew[710]
etanew[711]
etanew[712]
etanew[713]
etanew[714]
etanew[715]
etanew[716]
etanew[717]
etanew[718]
etanew[719]
etanew[720]
etanew[721]
etanew[722]
etanew[723]
etanew[724]
etanew[725]
etanew[726]
etanew[727]
etanew[728]
etanew[729]
etanew[730]
etanew[731]
etanew[732]
etanew[733]
etanew[734]
etanew[735]
etanew[736]
etanew[737]
etanew[738]
etanew[739]
etanew[740]
etanew[741]
etanew[742]
etanew[743]
etanew[744]
etanew[745]
etanew[746]
etanew[747]
etanew[748]
etanew[749]
etanew[750]
etanew[751]
etanew[752]
etanew[753]
etanew[754]
etanew[755]
etanew[756]
etanew[757]
etanew[758]
etanew[759]
etanew[760]
etanew[761]
etanew[762]
etanew[763]
etanew[764]
etanew[765]
etanew[766]
etanew[767]
etanew[768]
etanew[769]
etanew[770]
etanew[771]
etanew[772]
etanew[773]
etanew[774]
etanew[775]
etanew[776]
etanew[777]
etanew[778]
etanew[779]
etanew[780]
etanew[781]
etanew[782]
etanew[783]
etanew[784]
etanew[785]
etanew[786]
etanew[787]
etanew[788]
etanew[789]
etanew[790]
etanew[791]
etanew[792]
etanew[793]
etanew[794]
etanew[795]
etanew[796]
etanew[797]
etanew[798]
etanew[799]
etanew[800]
etanew[801]
etanew[802]
etanew[803]
etanew[804]
etanew[805]
etanew[806]
etanew[807]
etanew[808]
etanew[809]
etanew[810]
etanew[811]
etanew[812]
etanew[813]
etanew[814]
etanew[815]
etanew[816]
etanew[817]
etanew[818]
etanew[819]
etanew[820]
etanew[821]
etanew[822]
etanew[823]
etanew[824]
etanew[825]
etanew[826]
etanew[827]
etanew[828]
etanew[829]
etanew[830]
etanew[831]
etanew[832]
etanew[833]
etanew[834]
etanew[835]
etanew[836]
etanew[837]
etanew[838]
etanew[839]
etanew[840]
etanew[841]
etanew[842]
etanew[843]
etanew[844]
etanew[845]
etanew[846]
etanew[847]
etanew[848]
etanew[849]
etanew[850]
etanew[851]
etanew[852]
etanew[853]
etanew[854]
etanew[855]
etanew[856]
etanew[857]
etanew[858]
etanew[859]
etanew[860]
etanew[861]
etanew[862]
etanew[863]
etanew[864]
etanew[865]
etanew[866]
etanew[867]
etanew[868]
etanew[869]
etanew[870]
etanew[871]
etanew[872]
etanew[873]
etanew[874]
etanew[875]
etanew[876]
etanew[877]
etanew[878]
etanew[879]
etanew[880]
etanew[881]
etanew[882]
etanew[883]
etanew[884]
etanew[885]
etanew[886]
etanew[887]
etanew[888]
etanew[889]
etanew[890]
etanew[891]
etanew[892]
etanew[893]
etanew[894]
etanew[895]
etanew[896]
etanew[897]
etanew[898]
etanew[899]
etanew[900]
etanew[901]
etanew[902]
etanew[903]
etanew[904]
etanew[905]
etanew[906]
etanew[907]
etanew[908]
etanew[909]
etanew[910]
etanew[911]
etanew[912]
etanew[913]
etanew[914]
etanew[915]
etanew[916]
etanew[917]
etanew[918]
etanew[919]
etanew[920]
etanew[921]
etanew[922]
etanew[923]
etanew[924]
etanew[925]
etanew[926]
etanew[927]
etanew[928]
etanew[929]
etanew[930]
etanew[931]
etanew[932]
etanew[933]
etanew[934]
etanew[935]
etanew[936]
etanew[937]
etanew[938]
etanew[939]
etanew[940]
etanew[941]
etanew[942]
etanew[943]
etanew[944]
etanew[945]
etanew[946]
etanew[947]
etanew[948]
etanew[949]
etanew[950]
etanew[951]
etanew[952]
etanew[953]
etanew[954]
etanew[955]
etanew[956]
etanew[957]
etanew[958]
etanew[959]
etanew[960]
etanew[961]
etanew[962]
etanew[963]
etanew[964]
etanew[965]
etanew[966]
etanew[967]
etanew[968]
etanew[969]
etanew[970]
etanew[971]
etanew[972]
etanew[973]
etanew[974]
etanew[975]
etanew[976]
etanew[977]
etanew[978]
etanew[979]
etanew[980]
etanew[981]
etanew[982]
etanew[983]
etanew[984]
etanew[985]
etanew[986]
etanew[987]
etanew[988]
etanew[989]
etanew[990]
etanew[991]
etanew[992]
etanew[993]
etanew[994]
etanew[995]
etanew[996]
etanew[997]
etanew[998]
etanew[999]
etanew[1000]
etanew[1001]
etanew[1002]
etanew[1003]
etanew[1004]
etanew[1005]
etanew[1006]
etanew[1007]
etanew[1008]
etanew[1009]
etanew[1010]
etanew[1011]
etanew[1012]
etanew[1013]
etanew[1014]
etanew[1015]
etanew[1016]
etanew[1017]
etanew[1018]
etanew[1019]
etanew[1020]
etanew[1021]
etanew[1022]
etanew[1023]
etanew[1024]
etanew[1025]
etanew[1026]
etanew[1027]
etanew[1028]
etanew[1029]
etanew[1030]
etanew[1031]
etanew[1032]
etanew[1033]
etanew[1034]
etanew[1035]
etanew[1036]
etanew[1037]
etanew[1038]
etanew[1039]
etanew[1040]
etanew[1041]
etanew[1042]
etanew[1043]
etanew[1044]
etanew[1045]
etanew[1046]
etanew[1047]
etanew[1048]
etanew[1049]
etanew[1050]
etanew[1051]
etanew[1052]
etanew[1053]
etanew[1054]
etanew[1055]
etanew[1056]
etanew[1057]
etanew[1058]
etanew[1059]
etanew[1060]
etanew[1061]
etanew[1062]
etanew[1063]
etanew[1064]
etanew[1065]
etanew[1066]
etanew[1067]
etanew[1068]
etanew[1069]
etanew[1070]
etanew[1071]
etanew[1072]
etanew[1073]
etanew[1074]
etanew[1075]
etanew[1076]
etanew[1077]
etanew[1078]
etanew[1079]
etanew[1080]
etanew[1081]
etanew[1082]
etanew[1083]
etanew[1084]
etanew[1085]
etanew[1086]
etanew[1087]
etanew[1088]
etanew[1089]
etanew[1090]
etanew[1091]
etanew[1092]
etanew[1093]
etanew[1094]
etanew[1095]
etanew[1096]
etanew[1097]
etanew[1098]
etanew[1099]
etanew[1100]
etanew[1101]
etanew[1102]
etanew[1103]
etanew[1104]
etanew[1105]
etanew[1106]
etanew[1107]
etanew[1108]
etanew[1109]
etanew[1110]
etanew[1111]
etanew[1112]
etanew[1113]
etanew[1114]
etanew[1115]
etanew[1116]
etanew[1117]
etanew[1118]
etanew[1119]
etanew[1120]
etanew[1121]
etanew[1122]
etanew[1123]
etanew[1124]
etanew[1125]
etanew[1126]
etanew[1127]
etanew[1128]
etanew[1129]
etanew[1130]
etanew[1131]
etanew[1132]
etanew[1133]
etanew[1134]
etanew[1135]
etanew[1136]
etanew[1137]
etanew[1138]
etanew[1139]
etanew[1140]
etanew[1141]
etanew[1142]
etanew[1143]
etanew[1144]
etanew[1145]
etanew[1146]
etanew[1147]
etanew[1148]
etanew[1149]
etanew[1150]
etanew[1151]
etanew[1152]
etanew[1153]
etanew[1154]
etanew[1155]
etanew[1156]
etanew[1157]
etanew[1158]
etanew[1159]
etanew[1160]
etanew[1161]
etanew[1162]
etanew[1163]
etanew[1164]
etanew[1165]
etanew[1166]
etanew[1167]
etanew[1168]
etanew[1169]
etanew[1170]
etanew[1171]
etanew[1172]
etanew[1173]
etanew[1174]
etanew[1175]
etanew[1176]
etanew[1177]
etanew[1178]
etanew[1179]
etanew[1180]
etanew[1181]
etanew[1182]
etanew[1183]
etanew[1184]
etanew[1185]
etanew[1186]
etanew[1187]
etanew[1188]
etanew[1189]
etanew[1190]
etanew[1191]
etanew[1192]
etanew[1193]
etanew[1194]
etanew[1195]
etanew[1196]
etanew[1197]
etanew[1198]
etanew[1199]
etanew[1200]
etanew[1201]
etanew[1202]
etanew[1203]
etanew[1204]
etanew[1205]
etanew[1206]
etanew[1207]
etanew[1208]
etanew[1209]
etanew[1210]
etanew[1211]
etanew[1212]
etanew[1213]
etanew[1214]
etanew[1215]
etanew[1216]
etanew[1217]
etanew[1218]
etanew[1219]
etanew[1220]
etanew[1221]
etanew[1222]
etanew[1223]
etanew[1224]
etanew[1225]
etanew[1226]
etanew[1227]
etanew[1228]
etanew[1229]
etanew[1230]
etanew[1231]
etanew[1232]
etanew[1233]
etanew[1234]
etanew[1235]
etanew[1236]
etanew[1237]
etanew[1238]
etanew[1239]
etanew[1240]
etanew[1241]
etanew[1242]
etanew[1243]
etanew[1244]
etanew[1245]
etanew[1246]
etanew[1247]
etanew[1248]
etanew[1249]
etanew[1250]
etanew[1251]
etanew[1252]
etanew[1253]
etanew[1254]
etanew[1255]
etanew[1256]
etanew[1257]
etanew[1258]
etanew[1259]
etanew[1260]
etanew[1261]
etanew[1262]
etanew[1263]
etanew[1264]
etanew[1265]
etanew[1266]
etanew[1267]
etanew[1268]
etanew[1269]
etanew[1270]
etanew[1271]
etanew[1272]
etanew[1273]
etanew[1274]
etanew[1275]
etanew[1276]
etanew[1277]
etanew[1278]
etanew[1279]
etanew[1280]
etanew[1281]
etanew[1282]
etanew[1283]
etanew[1284]
etanew[1285]
etanew[1286]
etanew[1287]
etanew[1288]
etanew[1289]
etanew[1290]
etanew[1291]
etanew[1292]
etanew[1293]
etanew[1294]
etanew[1295]
etanew[1296]
etanew[1297]
etanew[1298]
etanew[1299]
etanew[1300]
etanew[1301]
etanew[1302]
etanew[1303]
etanew[1304]
etanew[1305]
etanew[1306]
etanew[1307]
etanew[1308]
etanew[1309]
etanew[1310]
etanew[1311]
etanew[1312]
etanew[1313]
etanew[1314]
etanew[1315]
etanew[1316]
etanew[1317]
etanew[1318]
etanew[1319]
etanew[1320]
etanew[1321]
etanew[1322]
etanew[1323]
etanew[1324]
etanew[1325]
etanew[1326]
etanew[1327]
etanew[1328]
etanew[1329]
etanew[1330]
etanew[1331]
etanew[1332]
etanew[1333]
etanew[1334]
etanew[1335]
etanew[1336]
etanew[1337]
etanew[1338]
etanew[1339]
etanew[1340]
etanew[1341]
etanew[1342]
etanew[1343]
etanew[1344]
etanew[1345]
etanew[1346]
etanew[1347]
etanew[1348]
etanew[1349]
etanew[1350]
etanew[1351]
etanew[1352]
etanew[1353]
etanew[1354]
etanew[1355]
etanew[1356]
etanew[1357]
etanew[1358]
etanew[1359]
etanew[1360]
etanew[1361]
etanew[1362]
etanew[1363]
etanew[1364]
etanew[1365]
etanew[1366]
etanew[1367]
etanew[1368]
etanew[1369]
etanew[1370]
etanew[1371]
etanew[1372]
etanew[1373]
etanew[1374]
etanew[1375]
etanew[1376]
etanew[1377]
etanew[1378]
etanew[1379]
etanew[1380]
etanew[1381]
etanew[1382]
etanew[1383]
etanew[1384]
etanew[1385]
etanew[1386]
etanew[1387]
etanew[1388]
etanew[1389]
etanew[1390]
etanew[1391]
etanew[1392]
etanew[1393]
etanew[1394]
etanew[1395]
etanew[1396]
etanew[1397]
etanew[1398]
etanew[1399]
etanew[1400]
etanew[1401]
etanew[1402]
etanew[1403]
etanew[1404]
etanew[1405]
etanew[1406]
etanew[1407]
etanew[1408]
etanew[1409]
etanew[1410]
etanew[1411]
etanew[1412]
etanew[1413]
etanew[1414]
etanew[1415]
etanew[1416]
etanew[1417]
etanew[1418]
etanew[1419]
etanew[1420]
etanew[1421]
etanew[1422]
etanew[1423]
etanew[1424]
etanew[1425]
etanew[1426]
etanew[1427]
etanew[1428]
etanew[1429]
etanew[1430]
etanew[1431]
etanew[1432]
etanew[1433]
etanew[1434]
etanew[1435]
etanew[1436]
etanew[1437]
etanew[1438]
etanew[1439]
etanew[1440]
etanew[1441]
etanew[1442]
etanew[1443]
etanew[1444]
etanew[1445]
etanew[1446]
etanew[1447]
etanew[1448]
etanew[1449]
etanew[1450]
etanew[1451]
etanew[1452]
etanew[1453]
etanew[1454]
etanew[1455]
etanew[1456]
etanew[1457]
etanew[1458]
etanew[1459]
etanew[1460]
etanew[1461]
etanew[1462]
etanew[1463]
etanew[1464]
etanew[1465]
etanew[1466]
etanew[1467]
etanew[1468]
etanew[1469]
etanew[1470]
etanew[1471]
etanew[1472]
etanew[1473]
etanew[1474]
etanew[1475]
etanew[1476]
etanew[1477]
etanew[1478]
etanew[1479]
etanew[1480]
etanew[1481]
etanew[1482]
etanew[1483]
etanew[1484]
etanew[1485]
etanew[1486]
etanew[1487]
etanew[1488]
etanew[1489]
etanew[1490]
etanew[1491]
etanew[1492]
etanew[1493]
etanew[1494]
etanew[1495]
etanew[1496]
etanew[1497]
etanew[1498]
etanew[1499]
etanew[1500]
etanew[1501]
etanew[1502]
etanew[1503]
etanew[1504]
etanew[1505]
etanew[1506]
etanew[1507]
etanew[1508]
etanew[1509]
etanew[1510]
etanew[1511]
etanew[1512]
etanew[1513]
etanew[1514]
etanew[1515]
etanew[1516]
etanew[1517]
etanew[1518]
etanew[1519]
etanew[1520]
etanew[1521]
etanew[1522]
etanew[1523]
etanew[1524]
etanew[1525]
etanew[1526]
etanew[1527]
etanew[1528]
etanew[1529]
etanew[1530]
etanew[1531]
etanew[1532]
etanew[1533]
etanew[1534]
etanew[1535]
etanew[1536]
etanew[1537]
etanew[1538]
etanew[1539]
etanew[1540]
etanew[1541]
etanew[1542]
etanew[1543]
etanew[1544]
etanew[1545]
etanew[1546]
etanew[1547]
etanew[1548]
etanew[1549]
etanew[1550]
etanew[1551]
etanew[1552]
etanew[1553]
etanew[1554]
etanew[1555]
etanew[1556]
etanew[1557]
etanew[1558]
etanew[1559]
etanew[1560]
etanew[1561]
etanew[1562]
etanew[1563]
etanew[1564]
etanew[1565]
etanew[1566]
etanew[1567]
etanew[1568]
etanew[1569]
etanew[1570]
etanew[1571]
etanew[1572]
etanew[1573]
etanew[1574]
etanew[1575]
etanew[1576]
etanew[1577]
etanew[1578]
etanew[1579]
etanew[1580]
etanew[1581]
etanew[1582]
etanew[1583]
etanew[1584]
etanew[1585]
etanew[1586]
etanew[1587]
etanew[1588]
etanew[1589]
etanew[1590]
etanew[1591]
etanew[1592]
etanew[1593]
etanew[1594]
etanew[1595]
etanew[1596]
etanew[1597]
etanew[1598]
etanew[1599]
etanew[1600]
etanew[1601]
etanew[1602]
etanew[1603]
etanew[1604]
etanew[1605]
etanew[1606]
etanew[1607]
etanew[1608]
etanew[1609]
etanew[1610]
etanew[1611]
etanew[1612]
etanew[1613]
etanew[1614]
etanew[1615]
etanew[1616]
etanew[1617]
etanew[1618]
etanew[1619]
etanew[1620]
etanew[1621]
etanew[1622]
etanew[1623]
etanew[1624]
etanew[1625]
etanew[1626]
etanew[1627]
etanew[1628]
etanew[1629]
etanew[1630]
etanew[1631]
etanew[1632]
etanew[1633]
etanew[1634]
etanew[1635]
etanew[1636]
etanew[1637]
etanew[1638]
etanew[1639]
etanew[1640]
etanew[1641]
etanew[1642]
etanew[1643]
etanew[1644]
etanew[1645]
etanew[1646]
etanew[1647]
etanew[1648]
etanew[1649]
etanew[1650]
etanew[1651]
etanew[1652]
etanew[1653]
etanew[1654]
etanew[1655]
etanew[1656]
etanew[1657]
etanew[1658]
etanew[1659]
etanew[1660]
etanew[1661]
etanew[1662]
etanew[1663]
etanew[1664]
etanew[1665]
etanew[1666]
etanew[1667]
etanew[1668]
etanew[1669]
etanew[1670]
etanew[1671]
etanew[1672]
etanew[1673]
etanew[1674]
etanew[1675]
etanew[1676]
etanew[1677]
etanew[1678]
etanew[1679]
etanew[1680]
etanew[1681]
etanew[1682]
etanew[1683]
etanew[1684]
etanew[1685]
etanew[1686]
etanew[1687]
etanew[1688]
etanew[1689]
etanew[1690]
etanew[1691]
etanew[1692]
etanew[1693]
etanew[1694]
etanew[1695]
etanew[1696]
etanew[1697]
etanew[1698]
etanew[1699]
etanew[1700]
etanew[1701]
etanew[1702]
etanew[1703]
etanew[1704]
etanew[1705]
etanew[1706]
etanew[1707]
etanew[1708]
etanew[1709]
etanew[1710]
etanew[1711]
etanew[1712]
etanew[1713]
etanew[1714]
etanew[1715]
etanew[1716]
etanew[1717]
etanew[1718]
etanew[1719]
etanew[1720]
etanew[1721]
etanew[1722]
etanew[1723]
etanew[1724]
etanew[1725]
etanew[1726]
etanew[1727]
etanew[1728]
etanew[1729]
etanew[1730]
etanew[1731]
etanew[1732]
etanew[1733]
etanew[1734]
etanew[1735]
etanew[1736]
etanew[1737]
etanew[1738]
etanew[1739]
etanew[1740]
etanew[1741]
etanew[1742]
etanew[1743]
etanew[1744]
etanew[1745]
etanew[1746]
etanew[1747]
etanew[1748]
etanew[1749]
etanew[1750]
etanew[1751]
etanew[1752]
etanew[1753]
etanew[1754]
etanew[1755]
etanew[1756]
etanew[1757]
etanew[1758]
etanew[1759]
etanew[1760]
etanew[1761]
etanew[1762]
etanew[1763]
etanew[1764]
etanew[1765]
etanew[1766]
etanew[1767]
etanew[1768]
etanew[1769]
etanew[1770]
etanew[1771]
etanew[1772]
etanew[1773]
etanew[1774]
etanew[1775]
etanew[1776]
etanew[1777]
etanew[1778]
etanew[1779]
etanew[1780]
etanew[1781]
etanew[1782]
etanew[1783]
etanew[1784]
etanew[1785]
etanew[1786]
etanew[1787]
etanew[1788]
etanew[1789]
etanew[1790]
etanew[1791]
etanew[1792]
etanew[1793]
etanew[1794]
etanew[1795]
etanew[1796]
etanew[1797]
etanew[1798]
etanew[1799]
etanew[1800]
etanew[1801]
etanew[1802]
etanew[1803]
etanew[1804]
etanew[1805]
etanew[1806]
etanew[1807]
etanew[1808]
etanew[1809]
etanew[1810]
etanew[1811]
etanew[1812]
etanew[1813]
etanew[1814]
etanew[1815]
etanew[1816]
etanew[1817]
etanew[1818]
etanew[1819]
etanew[1820]
etanew[1821]
etanew[1822]
etanew[1823]
etanew[1824]
etanew[1825]
etanew[1826]
etanew[1827]
etanew[1828]
etanew[1829]
etanew[1830]
etanew[1831]
etanew[1832]
etanew[1833]
etanew[1834]
etanew[1835]
etanew[1836]
etanew[1837]
etanew[1838]
etanew[1839]
etanew[1840]
etanew[1841]
etanew[1842]
etanew[1843]
etanew[1844]
etanew[1845]
etanew[1846]
etanew[1847]
etanew[1848]
etanew[1849]
etanew[1850]
etanew[1851]
etanew[1852]
etanew[1853]
etanew[1854]
etanew[1855]
etanew[1856]
etanew[1857]
etanew[1858]
etanew[1859]
etanew[1860]
etanew[1861]
etanew[1862]
etanew[1863]
etanew[1864]
etanew[1865]
etanew[1866]
etanew[1867]
etanew[1868]
etanew[1869]
etanew[1870]
etanew[1871]
etanew[1872]
etanew[1873]
etanew[1874]
etanew[1875]
etanew[1876]
etanew[1877]
etanew[1878]
etanew[1879]
etanew[1880]
etanew[1881]
etanew[1882]
etanew[1883]
etanew[1884]
etanew[1885]
etanew[1886]
etanew[1887]
etanew[1888]
etanew[1889]
etanew[1890]
etanew[1891]
etanew[1892]
etanew[1893]
etanew[1894]
etanew[1895]
etanew[1896]
etanew[1897]
etanew[1898]
etanew[1899]
etanew[1900]
etanew[1901]
etanew[1902]
etanew[1903]
etanew[1904]
etanew[1905]
etanew[1906]
etanew[1907]
etanew[1908]
etanew[1909]
etanew[1910]
etanew[1911]
etanew[1912]
etanew[1913]
etanew[1914]
etanew[1915]
etanew[1916]
etanew[1917]
etanew[1918]
etanew[1919]
etanew[1920]
etanew[1921]
etanew[1922]
etanew[1923]
etanew[1924]
etanew[1925]
etanew[1926]
etanew[1927]
etanew[1928]
etanew[1929]
etanew[1930]
etanew[1931]
etanew[1932]
etanew[1933]
etanew[1934]
etanew[1935]
etanew[1936]
etanew[1937]
etanew[1938]
etanew[1939]
etanew[1940]
etanew[1941]
etanew[1942]
etanew[1943]
etanew[1944]
etanew[1945]
etanew[1946]
etanew[1947]
etanew[1948]
etanew[1949]
etanew[1950]
etanew[1951]
etanew[1952]
etanew[1953]
etanew[1954]
etanew[1955]
etanew[1956]
etanew[1957]
etanew[1958]
etanew[1959]
etanew[1960]
etanew[1961]
etanew[1962]
etanew[1963]
etanew[1964]
etanew[1965]
etanew[1966]
etanew[1967]
etanew[1968]
etanew[1969]
etanew[1970]
etanew[1971]
etanew[1972]
etanew[1973]
etanew[1974]
etanew[1975]
etanew[1976]
etanew[1977]
etanew[1978]
etanew[1979]
etanew[1980]
etanew[1981]
etanew[1982]
etanew[1983]
etanew[1984]
etanew[1985]
etanew[1986]
etanew[1987]
etanew[1988]
etanew[1989]
etanew[1990]
etanew[1991]
etanew[1992]
etanew[1993]
etanew[1994]
etanew[1995]
etanew[1996]
etanew[1997]
etanew[1998]
etanew[1999]
etanew[2000]
etanew[2001]
etanew[2002]
etanew[2003]
etanew[2004]
etanew[2005]
etanew[2006]
etanew[2007]
etanew[2008]
etanew[2009]
etanew[2010]
etanew[2011]
etanew[2012]
etanew[2013]
etanew[2014]
etanew[2015]
etanew[2016]
etanew[2017]
etanew[2018]
etanew[2019]
etanew[2020]
etanew[2021]
etanew[2022]
etanew[2023]
etanew[2024]
etanew[2025]
etanew[2026]
etanew[2027]
etanew[2028]
etanew[2029]
etanew[2030]
etanew[2031]
etanew[2032]
etanew[2033]
etanew[2034]
etanew[2035]
etanew[2036]
etanew[2037]
etanew[2038]
etanew[2039]
etanew[2040]
etanew[2041]
etanew[2042]
etanew[2043]
etanew[2044]
etanew[2045]
etanew[2046]
etanew[2047]
etanew[2048]
etanew[2049]
etanew[2050]
etanew[2051]
etanew[2052]
etanew[2053]
etanew[2054]
etanew[2055]
etanew[2056]
etanew[2057]
etanew[2058]
etanew[2059]
etanew[2060]
etanew[2061]
etanew[2062]
etanew[2063]
etanew[2064]
etanew[2065]
etanew[2066]
etanew[2067]
etanew[2068]
etanew[2069]
etanew[2070]
etanew[2071]
etanew[2072]
etanew[2073]
etanew[2074]
etanew[2075]
etanew[2076]
etanew[2077]
etanew[2078]
etanew[2079]
etanew[2080]
etanew[2081]
etanew[2082]
etanew[2083]
etanew[2084]
etanew[2085]
etanew[2086]
etanew[2087]
etanew[2088]
etanew[2089]
etanew[2090]
etanew[2091]
etanew[2092]
etanew[2093]
etanew[2094]
etanew[2095]
etanew[2096]
etanew[2097]
etanew[2098]
etanew[2099]
etanew[2100]
etanew[2101]
etanew[2102]
etanew[2103]
etanew[2104]
etanew[2105]
etanew[2106]
etanew[2107]
etanew[2108]
etanew[2109]
etanew[2110]
etanew[2111]
etanew[2112]
etanew[2113]
etanew[2114]
etanew[2115]
etanew[2116]
etanew[2117]
etanew[2118]
etanew[2119]
etanew[2120]
etanew[2121]
etanew[2122]
etanew[2123]
etanew[2124]
etanew[2125]
etanew[2126]
etanew[2127]
etanew[2128]
etanew[2129]
etanew[2130]
etanew[2131]
etanew[2132]
etanew[2133]
etanew[2134]
etanew[2135]
etanew[2136]
etanew[2137]
etanew[2138]
etanew[2139]
etanew[2140]
etanew[2141]
etanew[2142]
etanew[2143]
etanew[2144]
etanew[2145]
etanew[2146]
etanew[2147]
etanew[2148]
etanew[2149]
etanew[2150]
etanew[2151]
etanew[2152]
etanew[2153]
etanew[2154]
etanew[2155]
etanew[2156]
etanew[2157]
etanew[2158]
etanew[2159]
etanew[2160]
etanew[2161]
etanew[2162]
etanew[2163]
etanew[2164]
etanew[2165]
etanew[2166]
etanew[2167]
etanew[2168]
etanew[2169]
etanew[2170]
etanew[2171]
etanew[2172]
etanew[2173]
etanew[2174]
etanew[2175]
etanew[2176]
etanew[2177]
etanew[2178]
etanew[2179]
etanew[2180]
etanew[2181]
etanew[2182]
etanew[2183]
etanew[2184]
etanew[2185]
etanew[2186]
etanew[2187]
etanew[2188]
etanew[2189]
etanew[2190]
etanew[2191]
etanew[2192]
etanew[2193]
etanew[2194]
etanew[2195]
etanew[2196]
etanew[2197]
etanew[2198]
etanew[2199]
etanew[2200]
etanew[2201]
etanew[2202]
etanew[2203]
etanew[2204]
etanew[2205]
etanew[2206]
etanew[2207]
etanew[2208]
etanew[2209]
etanew[2210]
etanew[2211]
etanew[2212]
etanew[2213]
etanew[2214]
etanew[2215]
etanew[2216]
etanew[2217]
etanew[2218]
etanew[2219]
etanew[2220]
etanew[2221]
etanew[2222]
etanew[2223]
etanew[2224]
etanew[2225]
etanew[2226]
etanew[2227]
etanew[2228]
etanew[2229]
etanew[2230]
etanew[2231]
etanew[2232]
etanew[2233]
etanew[2234]
etanew[2235]
etanew[2236]
etanew[2237]
etanew[2238]
etanew[2239]
etanew[2240]
etanew[2241]
etanew[2242]
etanew[2243]
etanew[2244]
etanew[2245]
etanew[2246]
etanew[2247]
etanew[2248]
etanew[2249]
etanew[2250]
etanew[2251]
etanew[2252]
etanew[2253]
etanew[2254]
etanew[2255]
etanew[2256]
etanew[2257]
etanew[2258]
etanew[2259]
etanew[2260]
etanew[2261]
etanew[2262]
etanew[2263]
etanew[2264]
etanew[2265]
etanew[2266]
etanew[2267]
etanew[2268]
etanew[2269]
etanew[2270]
etanew[2271]
etanew[2272]
etanew[2273]
etanew[2274]
etanew[2275]
etanew[2276]
etanew[2277]
etanew[2278]
etanew[2279]
etanew[2280]
etanew[2281]
etanew[2282]
etanew[2283]
etanew[2284]
etanew[2285]
etanew[2286]
etanew[2287]
etanew[2288]
etanew[2289]
etanew[2290]
etanew[2291]
etanew[2292]
etanew[2293]
etanew[2294]
etanew[2295]
etanew[2296]
etanew[2297]
etanew[2298]
etanew[2299]
etanew[2300]
etanew[2301]
etanew[2302]
etanew[2303]
etanew[2304]
etanew[2305]
etanew[2306]
etanew[2307]
etanew[2308]
etanew[2309]
etanew[2310]
etanew[2311]
etanew[2312]
etanew[2313]
etanew[2314]
etanew[2315]
etanew[2316]
etanew[2317]
etanew[2318]
etanew[2319]
etanew[2320]
etanew[2321]
etanew[2322]
etanew[2323]
etanew[2324]
etanew[2325]
etanew[2326]
etanew[2327]
etanew[2328]
etanew[2329]
etanew[2330]
etanew[2331]
etanew[2332]
etanew[2333]
etanew[2334]
etanew[2335]
etanew[2336]
etanew[2337]
etanew[2338]
etanew[2339]
etanew[2340]
etanew[2341]
etanew[2342]
etanew[2343]
etanew[2344]
etanew[2345]
etanew[2346]
etanew[2347]
etanew[2348]
etanew[2349]
etanew[2350]
etanew[2351]
etanew[2352]
etanew[2353]
etanew[2354]
etanew[2355]
etanew[2356]
etanew[2357]
etanew[2358]
etanew[2359]
etanew[2360]
etanew[2361]
etanew[2362]
etanew[2363]
etanew[2364]
etanew[2365]
etanew[2366]
etanew[2367]
etanew[2368]
etanew[2369]
etanew[2370]
etanew[2371]
etanew[2372]
etanew[2373]
etanew[2374]
etanew[2375]
etanew[2376]
etanew[2377]
etanew[2378]
etanew[2379]
etanew[2380]
etanew[2381]
etanew[2382]
etanew[2383]
etanew[2384]
etanew[2385]
etanew[2386]
etanew[2387]
etanew[2388]
etanew[2389]
etanew[2390]
etanew[2391]
etanew[2392]
etanew[2393]
etanew[2394]
etanew[2395]
etanew[2396]
etanew[2397]
etanew[2398]
etanew[2399]
etanew[2400]
etanew[2401]
etanew[2402]
etanew[2403]
etanew[2404]
etanew[2405]
etanew[2406]
etanew[2407]
etanew[2408]
etanew[2409]
etanew[2410]
etanew[2411]
etanew[2412]
etanew[2413]
etanew[2414]
etanew[2415]
etanew[2416]
etanew[2417]
etanew[2418]
etanew[2419]
etanew[2420]
etanew[2421]
etanew[2422]
etanew[2423]
etanew[2424]
etanew[2425]
etanew[2426]
etanew[2427]
etanew[2428]
etanew[2429]
etanew[2430]
etanew[2431]
etanew[2432]
etanew[2433]
etanew[2434]
etanew[2435]
etanew[2436]
etanew[2437]
etanew[2438]
etanew[2439]
etanew[2440]
etanew[2441]
etanew[2442]
etanew[2443]
etanew[2444]
etanew[2445]
etanew[2446]
etanew[2447]
etanew[2448]
etanew[2449]
etanew[2450]
etanew[2451]
etanew[2452]
etanew[2453]
etanew[2454]
etanew[2455]
etanew[2456]
etanew[2457]
etanew[2458]
etanew[2459]
etanew[2460]
etanew[2461]
etanew[2462]
etanew[2463]
etanew[2464]
etanew[2465]
etanew[2466]
etanew[2467]
etanew[2468]
etanew[2469]
etanew[2470]
etanew[2471]
etanew[2472]
etanew[2473]
etanew[2474]
etanew[2475]
etanew[2476]
etanew[2477]
etanew[2478]
etanew[2479]
etanew[2480]
etanew[2481]
etanew[2482]
etanew[2483]
etanew[2484]
etanew[2485]
etanew[2486]
etanew[2487]
etanew[2488]
etanew[2489]
etanew[2490]
etanew[2491]
etanew[2492]
etanew[2493]
etanew[2494]
etanew[2495]
etanew[2496]
etanew[2497]
etanew[2498]
etanew[2499]
pnew[0]
pnew[1]
pnew[2]
pnew[3]
pnew[4]
pnew[5]
pnew[6]
pnew[7]
pnew[8]
pnew[9]
pnew[10]
pnew[11]
pnew[12]
pnew[13]
pnew[14]
pnew[15]
pnew[16]
pnew[17]
pnew[18]
pnew[19]
pnew[20]
pnew[21]
pnew[22]
pnew[23]
pnew[24]
pnew[25]
pnew[26]
pnew[27]
pnew[28]
pnew[29]
pnew[30]
pnew[31]
pnew[32]
pnew[33]
pnew[34]
pnew[35]
pnew[36]
pnew[37]
pnew[38]
pnew[39]
pnew[40]
pnew[41]
pnew[42]
pnew[43]
pnew[44]
pnew[45]
pnew[46]
pnew[47]
pnew[48]
pnew[49]
pnew[50]
pnew[51]
pnew[52]
pnew[53]
pnew[54]
pnew[55]
pnew[56]
pnew[57]
pnew[58]
pnew[59]
pnew[60]
pnew[61]
pnew[62]
pnew[63]
pnew[64]
pnew[65]
pnew[66]
pnew[67]
pnew[68]
pnew[69]
pnew[70]
pnew[71]
pnew[72]
pnew[73]
pnew[74]
pnew[75]
pnew[76]
pnew[77]
pnew[78]
pnew[79]
pnew[80]
pnew[81]
pnew[82]
pnew[83]
pnew[84]
pnew[85]
pnew[86]
pnew[87]
pnew[88]
pnew[89]
pnew[90]
pnew[91]
pnew[92]
pnew[93]
pnew[94]
pnew[95]
pnew[96]
pnew[97]
pnew[98]
pnew[99]
pnew[100]
pnew[101]
pnew[102]
pnew[103]
pnew[104]
pnew[105]
pnew[106]
pnew[107]
pnew[108]
pnew[109]
pnew[110]
pnew[111]
pnew[112]
pnew[113]
pnew[114]
pnew[115]
pnew[116]
pnew[117]
pnew[118]
pnew[119]
pnew[120]
pnew[121]
pnew[122]
pnew[123]
pnew[124]
pnew[125]
pnew[126]
pnew[127]
pnew[128]
pnew[129]
pnew[130]
pnew[131]
pnew[132]
pnew[133]
pnew[134]
pnew[135]
pnew[136]
pnew[137]
pnew[138]
pnew[139]
pnew[140]
pnew[141]
pnew[142]
pnew[143]
pnew[144]
pnew[145]
pnew[146]
pnew[147]
pnew[148]
pnew[149]
pnew[150]
pnew[151]
pnew[152]
pnew[153]
pnew[154]
pnew[155]
pnew[156]
pnew[157]
pnew[158]
pnew[159]
pnew[160]
pnew[161]
pnew[162]
pnew[163]
pnew[164]
pnew[165]
pnew[166]
pnew[167]
pnew[168]
pnew[169]
pnew[170]
pnew[171]
pnew[172]
pnew[173]
pnew[174]
pnew[175]
pnew[176]
pnew[177]
pnew[178]
pnew[179]
pnew[180]
pnew[181]
pnew[182]
pnew[183]
pnew[184]
pnew[185]
pnew[186]
pnew[187]
pnew[188]
pnew[189]
pnew[190]
pnew[191]
pnew[192]
pnew[193]
pnew[194]
pnew[195]
pnew[196]
pnew[197]
pnew[198]
pnew[199]
pnew[200]
pnew[201]
pnew[202]
pnew[203]
pnew[204]
pnew[205]
pnew[206]
pnew[207]
pnew[208]
pnew[209]
pnew[210]
pnew[211]
pnew[212]
pnew[213]
pnew[214]
pnew[215]
pnew[216]
pnew[217]
pnew[218]
pnew[219]
pnew[220]
pnew[221]
pnew[222]
pnew[223]
pnew[224]
pnew[225]
pnew[226]
pnew[227]
pnew[228]
pnew[229]
pnew[230]
pnew[231]
pnew[232]
pnew[233]
pnew[234]
pnew[235]
pnew[236]
pnew[237]
pnew[238]
pnew[239]
pnew[240]
pnew[241]
pnew[242]
pnew[243]
pnew[244]
pnew[245]
pnew[246]
pnew[247]
pnew[248]
pnew[249]
pnew[250]
pnew[251]
pnew[252]
pnew[253]
pnew[254]
pnew[255]
pnew[256]
pnew[257]
pnew[258]
pnew[259]
pnew[260]
pnew[261]
pnew[262]
pnew[263]
pnew[264]
pnew[265]
pnew[266]
pnew[267]
pnew[268]
pnew[269]
pnew[270]
pnew[271]
pnew[272]
pnew[273]
pnew[274]
pnew[275]
pnew[276]
pnew[277]
pnew[278]
pnew[279]
pnew[280]
pnew[281]
pnew[282]
pnew[283]
pnew[284]
pnew[285]
pnew[286]
pnew[287]
pnew[288]
pnew[289]
pnew[290]
pnew[291]
pnew[292]
pnew[293]
pnew[294]
pnew[295]
pnew[296]
pnew[297]
pnew[298]
pnew[299]
pnew[300]
pnew[301]
pnew[302]
pnew[303]
pnew[304]
pnew[305]
pnew[306]
pnew[307]
pnew[308]
pnew[309]
pnew[310]
pnew[311]
pnew[312]
pnew[313]
pnew[314]
pnew[315]
pnew[316]
pnew[317]
pnew[318]
pnew[319]
pnew[320]
pnew[321]
pnew[322]
pnew[323]
pnew[324]
pnew[325]
pnew[326]
pnew[327]
pnew[328]
pnew[329]
pnew[330]
pnew[331]
pnew[332]
pnew[333]
pnew[334]
pnew[335]
pnew[336]
pnew[337]
pnew[338]
pnew[339]
pnew[340]
pnew[341]
pnew[342]
pnew[343]
pnew[344]
pnew[345]
pnew[346]
pnew[347]
pnew[348]
pnew[349]
pnew[350]
pnew[351]
pnew[352]
pnew[353]
pnew[354]
pnew[355]
pnew[356]
pnew[357]
pnew[358]
pnew[359]
pnew[360]
pnew[361]
pnew[362]
pnew[363]
pnew[364]
pnew[365]
pnew[366]
pnew[367]
pnew[368]
pnew[369]
pnew[370]
pnew[371]
pnew[372]
pnew[373]
pnew[374]
pnew[375]
pnew[376]
pnew[377]
pnew[378]
pnew[379]
pnew[380]
pnew[381]
pnew[382]
pnew[383]
pnew[384]
pnew[385]
pnew[386]
pnew[387]
pnew[388]
pnew[389]
pnew[390]
pnew[391]
pnew[392]
pnew[393]
pnew[394]
pnew[395]
pnew[396]
pnew[397]
pnew[398]
pnew[399]
pnew[400]
pnew[401]
pnew[402]
pnew[403]
pnew[404]
pnew[405]
pnew[406]
pnew[407]
pnew[408]
pnew[409]
pnew[410]
pnew[411]
pnew[412]
pnew[413]
pnew[414]
pnew[415]
pnew[416]
pnew[417]
pnew[418]
pnew[419]
pnew[420]
pnew[421]
pnew[422]
pnew[423]
pnew[424]
pnew[425]
pnew[426]
pnew[427]
pnew[428]
pnew[429]
pnew[430]
pnew[431]
pnew[432]
pnew[433]
pnew[434]
pnew[435]
pnew[436]
pnew[437]
pnew[438]
pnew[439]
pnew[440]
pnew[441]
pnew[442]
pnew[443]
pnew[444]
pnew[445]
pnew[446]
pnew[447]
pnew[448]
pnew[449]
pnew[450]
pnew[451]
pnew[452]
pnew[453]
pnew[454]
pnew[455]
pnew[456]
pnew[457]
pnew[458]
pnew[459]
pnew[460]
pnew[461]
pnew[462]
pnew[463]
pnew[464]
pnew[465]
pnew[466]
pnew[467]
pnew[468]
pnew[469]
pnew[470]
pnew[471]
pnew[472]
pnew[473]
pnew[474]
pnew[475]
pnew[476]
pnew[477]
pnew[478]
pnew[479]
pnew[480]
pnew[481]
pnew[482]
pnew[483]
pnew[484]
pnew[485]
pnew[486]
pnew[487]
pnew[488]
pnew[489]
pnew[490]
pnew[491]
pnew[492]
pnew[493]
pnew[494]
pnew[495]
pnew[496]
pnew[497]
pnew[498]
pnew[499]
pnew[500]
pnew[501]
pnew[502]
pnew[503]
pnew[504]
pnew[505]
pnew[506]
pnew[507]
pnew[508]
pnew[509]
pnew[510]
pnew[511]
pnew[512]
pnew[513]
pnew[514]
pnew[515]
pnew[516]
pnew[517]
pnew[518]
pnew[519]
pnew[520]
pnew[521]
pnew[522]
pnew[523]
pnew[524]
pnew[525]
pnew[526]
pnew[527]
pnew[528]
pnew[529]
pnew[530]
pnew[531]
pnew[532]
pnew[533]
pnew[534]
pnew[535]
pnew[536]
pnew[537]
pnew[538]
pnew[539]
pnew[540]
pnew[541]
pnew[542]
pnew[543]
pnew[544]
pnew[545]
pnew[546]
pnew[547]
pnew[548]
pnew[549]
pnew[550]
pnew[551]
pnew[552]
pnew[553]
pnew[554]
pnew[555]
pnew[556]
pnew[557]
pnew[558]
pnew[559]
pnew[560]
pnew[561]
pnew[562]
pnew[563]
pnew[564]
pnew[565]
pnew[566]
pnew[567]
pnew[568]
pnew[569]
pnew[570]
pnew[571]
pnew[572]
pnew[573]
pnew[574]
pnew[575]
pnew[576]
pnew[577]
pnew[578]
pnew[579]
pnew[580]
pnew[581]
pnew[582]
pnew[583]
pnew[584]
pnew[585]
pnew[586]
pnew[587]
pnew[588]
pnew[589]
pnew[590]
pnew[591]
pnew[592]
pnew[593]
pnew[594]
pnew[595]
pnew[596]
pnew[597]
pnew[598]
pnew[599]
pnew[600]
pnew[601]
pnew[602]
pnew[603]
pnew[604]
pnew[605]
pnew[606]
pnew[607]
pnew[608]
pnew[609]
pnew[610]
pnew[611]
pnew[612]
pnew[613]
pnew[614]
pnew[615]
pnew[616]
pnew[617]
pnew[618]
pnew[619]
pnew[620]
pnew[621]
pnew[622]
pnew[623]
pnew[624]
pnew[625]
pnew[626]
pnew[627]
pnew[628]
pnew[629]
pnew[630]
pnew[631]
pnew[632]
pnew[633]
pnew[634]
pnew[635]
pnew[636]
pnew[637]
pnew[638]
pnew[639]
pnew[640]
pnew[641]
pnew[642]
pnew[643]
pnew[644]
pnew[645]
pnew[646]
pnew[647]
pnew[648]
pnew[649]
pnew[650]
pnew[651]
pnew[652]
pnew[653]
pnew[654]
pnew[655]
pnew[656]
pnew[657]
pnew[658]
pnew[659]
pnew[660]
pnew[661]
pnew[662]
pnew[663]
pnew[664]
pnew[665]
pnew[666]
pnew[667]
pnew[668]
pnew[669]
pnew[670]
pnew[671]
pnew[672]
pnew[673]
pnew[674]
pnew[675]
pnew[676]
pnew[677]
pnew[678]
pnew[679]
pnew[680]
pnew[681]
pnew[682]
pnew[683]
pnew[684]
pnew[685]
pnew[686]
pnew[687]
pnew[688]
pnew[689]
pnew[690]
pnew[691]
pnew[692]
pnew[693]
pnew[694]
pnew[695]
pnew[696]
pnew[697]
pnew[698]
pnew[699]
pnew[700]
pnew[701]
pnew[702]
pnew[703]
pnew[704]
pnew[705]
pnew[706]
pnew[707]
pnew[708]
pnew[709]
pnew[710]
pnew[711]
pnew[712]
pnew[713]
pnew[714]
pnew[715]
pnew[716]
pnew[717]
pnew[718]
pnew[719]
pnew[720]
pnew[721]
pnew[722]
pnew[723]
pnew[724]
pnew[725]
pnew[726]
pnew[727]
pnew[728]
pnew[729]
pnew[730]
pnew[731]
pnew[732]
pnew[733]
pnew[734]
pnew[735]
pnew[736]
pnew[737]
pnew[738]
pnew[739]
pnew[740]
pnew[741]
pnew[742]
pnew[743]
pnew[744]
pnew[745]
pnew[746]
pnew[747]
pnew[748]
pnew[749]
pnew[750]
pnew[751]
pnew[752]
pnew[753]
pnew[754]
pnew[755]
pnew[756]
pnew[757]
pnew[758]
pnew[759]
pnew[760]
pnew[761]
pnew[762]
pnew[763]
pnew[764]
pnew[765]
pnew[766]
pnew[767]
pnew[768]
pnew[769]
pnew[770]
pnew[771]
pnew[772]
pnew[773]
pnew[774]
pnew[775]
pnew[776]
pnew[777]
pnew[778]
pnew[779]
pnew[780]
pnew[781]
pnew[782]
pnew[783]
pnew[784]
pnew[785]
pnew[786]
pnew[787]
pnew[788]
pnew[789]
pnew[790]
pnew[791]
pnew[792]
pnew[793]
pnew[794]
pnew[795]
pnew[796]
pnew[797]
pnew[798]
pnew[799]
pnew[800]
pnew[801]
pnew[802]
pnew[803]
pnew[804]
pnew[805]
pnew[806]
pnew[807]
pnew[808]
pnew[809]
pnew[810]
pnew[811]
pnew[812]
pnew[813]
pnew[814]
pnew[815]
pnew[816]
pnew[817]
pnew[818]
pnew[819]
pnew[820]
pnew[821]
pnew[822]
pnew[823]
pnew[824]
pnew[825]
pnew[826]
pnew[827]
pnew[828]
pnew[829]
pnew[830]
pnew[831]
pnew[832]
pnew[833]
pnew[834]
pnew[835]
pnew[836]
pnew[837]
pnew[838]
pnew[839]
pnew[840]
pnew[841]
pnew[842]
pnew[843]
pnew[844]
pnew[845]
pnew[846]
pnew[847]
pnew[848]
pnew[849]
pnew[850]
pnew[851]
pnew[852]
pnew[853]
pnew[854]
pnew[855]
pnew[856]
pnew[857]
pnew[858]
pnew[859]
pnew[860]
pnew[861]
pnew[862]
pnew[863]
pnew[864]
pnew[865]
pnew[866]
pnew[867]
pnew[868]
pnew[869]
pnew[870]
pnew[871]
pnew[872]
pnew[873]
pnew[874]
pnew[875]
pnew[876]
pnew[877]
pnew[878]
pnew[879]
pnew[880]
pnew[881]
pnew[882]
pnew[883]
pnew[884]
pnew[885]
pnew[886]
pnew[887]
pnew[888]
pnew[889]
pnew[890]
pnew[891]
pnew[892]
pnew[893]
pnew[894]
pnew[895]
pnew[896]
pnew[897]
pnew[898]
pnew[899]
pnew[900]
pnew[901]
pnew[902]
pnew[903]
pnew[904]
pnew[905]
pnew[906]
pnew[907]
pnew[908]
pnew[909]
pnew[910]
pnew[911]
pnew[912]
pnew[913]
pnew[914]
pnew[915]
pnew[916]
pnew[917]
pnew[918]
pnew[919]
pnew[920]
pnew[921]
pnew[922]
pnew[923]
pnew[924]
pnew[925]
pnew[926]
pnew[927]
pnew[928]
pnew[929]
pnew[930]
pnew[931]
pnew[932]
pnew[933]
pnew[934]
pnew[935]
pnew[936]
pnew[937]
pnew[938]
pnew[939]
pnew[940]
pnew[941]
pnew[942]
pnew[943]
pnew[944]
pnew[945]
pnew[946]
pnew[947]
pnew[948]
pnew[949]
pnew[950]
pnew[951]
pnew[952]
pnew[953]
pnew[954]
pnew[955]
pnew[956]
pnew[957]
pnew[958]
pnew[959]
pnew[960]
pnew[961]
pnew[962]
pnew[963]
pnew[964]
pnew[965]
pnew[966]
pnew[967]
pnew[968]
pnew[969]
pnew[970]
pnew[971]
pnew[972]
pnew[973]
pnew[974]
pnew[975]
pnew[976]
pnew[977]
pnew[978]
pnew[979]
pnew[980]
pnew[981]
pnew[982]
pnew[983]
pnew[984]
pnew[985]
pnew[986]
pnew[987]
pnew[988]
pnew[989]
pnew[990]
pnew[991]
pnew[992]
pnew[993]
pnew[994]
pnew[995]
pnew[996]
pnew[997]
pnew[998]
pnew[999]
pnew[1000]
pnew[1001]
pnew[1002]
pnew[1003]
pnew[1004]
pnew[1005]
pnew[1006]
pnew[1007]
pnew[1008]
pnew[1009]
pnew[1010]
pnew[1011]
pnew[1012]
pnew[1013]
pnew[1014]
pnew[1015]
pnew[1016]
pnew[1017]
pnew[1018]
pnew[1019]
pnew[1020]
pnew[1021]
pnew[1022]
pnew[1023]
pnew[1024]
pnew[1025]
pnew[1026]
pnew[1027]
pnew[1028]
pnew[1029]
pnew[1030]
pnew[1031]
pnew[1032]
pnew[1033]
pnew[1034]
pnew[1035]
pnew[1036]
pnew[1037]
pnew[1038]
pnew[1039]
pnew[1040]
pnew[1041]
pnew[1042]
pnew[1043]
pnew[1044]
pnew[1045]
pnew[1046]
pnew[1047]
pnew[1048]
pnew[1049]
pnew[1050]
pnew[1051]
pnew[1052]
pnew[1053]
pnew[1054]
pnew[1055]
pnew[1056]
pnew[1057]
pnew[1058]
pnew[1059]
pnew[1060]
pnew[1061]
pnew[1062]
pnew[1063]
pnew[1064]
pnew[1065]
pnew[1066]
pnew[1067]
pnew[1068]
pnew[1069]
pnew[1070]
pnew[1071]
pnew[1072]
pnew[1073]
pnew[1074]
pnew[1075]
pnew[1076]
pnew[1077]
pnew[1078]
pnew[1079]
pnew[1080]
pnew[1081]
pnew[1082]
pnew[1083]
pnew[1084]
pnew[1085]
pnew[1086]
pnew[1087]
pnew[1088]
pnew[1089]
pnew[1090]
pnew[1091]
pnew[1092]
pnew[1093]
pnew[1094]
pnew[1095]
pnew[1096]
pnew[1097]
pnew[1098]
pnew[1099]
pnew[1100]
pnew[1101]
pnew[1102]
pnew[1103]
pnew[1104]
pnew[1105]
pnew[1106]
pnew[1107]
pnew[1108]
pnew[1109]
pnew[1110]
pnew[1111]
pnew[1112]
pnew[1113]
pnew[1114]
pnew[1115]
pnew[1116]
pnew[1117]
pnew[1118]
pnew[1119]
pnew[1120]
pnew[1121]
pnew[1122]
pnew[1123]
pnew[1124]
pnew[1125]
pnew[1126]
pnew[1127]
pnew[1128]
pnew[1129]
pnew[1130]
pnew[1131]
pnew[1132]
pnew[1133]
pnew[1134]
pnew[1135]
pnew[1136]
pnew[1137]
pnew[1138]
pnew[1139]
pnew[1140]
pnew[1141]
pnew[1142]
pnew[1143]
pnew[1144]
pnew[1145]
pnew[1146]
pnew[1147]
pnew[1148]
pnew[1149]
pnew[1150]
pnew[1151]
pnew[1152]
pnew[1153]
pnew[1154]
pnew[1155]
pnew[1156]
pnew[1157]
pnew[1158]
pnew[1159]
pnew[1160]
pnew[1161]
pnew[1162]
pnew[1163]
pnew[1164]
pnew[1165]
pnew[1166]
pnew[1167]
pnew[1168]
pnew[1169]
pnew[1170]
pnew[1171]
pnew[1172]
pnew[1173]
pnew[1174]
pnew[1175]
pnew[1176]
pnew[1177]
pnew[1178]
pnew[1179]
pnew[1180]
pnew[1181]
pnew[1182]
pnew[1183]
pnew[1184]
pnew[1185]
pnew[1186]
pnew[1187]
pnew[1188]
pnew[1189]
pnew[1190]
pnew[1191]
pnew[1192]
pnew[1193]
pnew[1194]
pnew[1195]
pnew[1196]
pnew[1197]
pnew[1198]
pnew[1199]
pnew[1200]
pnew[1201]
pnew[1202]
pnew[1203]
pnew[1204]
pnew[1205]
pnew[1206]
pnew[1207]
pnew[1208]
pnew[1209]
pnew[1210]
pnew[1211]
pnew[1212]
pnew[1213]
pnew[1214]
pnew[1215]
pnew[1216]
pnew[1217]
pnew[1218]
pnew[1219]
pnew[1220]
pnew[1221]
pnew[1222]
pnew[1223]
pnew[1224]
pnew[1225]
pnew[1226]
pnew[1227]
pnew[1228]
pnew[1229]
pnew[1230]
pnew[1231]
pnew[1232]
pnew[1233]
pnew[1234]
pnew[1235]
pnew[1236]
pnew[1237]
pnew[1238]
pnew[1239]
pnew[1240]
pnew[1241]
pnew[1242]
pnew[1243]
pnew[1244]
pnew[1245]
pnew[1246]
pnew[1247]
pnew[1248]
pnew[1249]
pnew[1250]
pnew[1251]
pnew[1252]
pnew[1253]
pnew[1254]
pnew[1255]
pnew[1256]
pnew[1257]
pnew[1258]
pnew[1259]
pnew[1260]
pnew[1261]
pnew[1262]
pnew[1263]
pnew[1264]
pnew[1265]
pnew[1266]
pnew[1267]
pnew[1268]
pnew[1269]
pnew[1270]
pnew[1271]
pnew[1272]
pnew[1273]
pnew[1274]
pnew[1275]
pnew[1276]
pnew[1277]
pnew[1278]
pnew[1279]
pnew[1280]
pnew[1281]
pnew[1282]
pnew[1283]
pnew[1284]
pnew[1285]
pnew[1286]
pnew[1287]
pnew[1288]
pnew[1289]
pnew[1290]
pnew[1291]
pnew[1292]
pnew[1293]
pnew[1294]
pnew[1295]
pnew[1296]
pnew[1297]
pnew[1298]
pnew[1299]
pnew[1300]
pnew[1301]
pnew[1302]
pnew[1303]
pnew[1304]
pnew[1305]
pnew[1306]
pnew[1307]
pnew[1308]
pnew[1309]
pnew[1310]
pnew[1311]
pnew[1312]
pnew[1313]
pnew[1314]
pnew[1315]
pnew[1316]
pnew[1317]
pnew[1318]
pnew[1319]
pnew[1320]
pnew[1321]
pnew[1322]
pnew[1323]
pnew[1324]
pnew[1325]
pnew[1326]
pnew[1327]
pnew[1328]
pnew[1329]
pnew[1330]
pnew[1331]
pnew[1332]
pnew[1333]
pnew[1334]
pnew[1335]
pnew[1336]
pnew[1337]
pnew[1338]
pnew[1339]
pnew[1340]
pnew[1341]
pnew[1342]
pnew[1343]
pnew[1344]
pnew[1345]
pnew[1346]
pnew[1347]
pnew[1348]
pnew[1349]
pnew[1350]
pnew[1351]
pnew[1352]
pnew[1353]
pnew[1354]
pnew[1355]
pnew[1356]
pnew[1357]
pnew[1358]
pnew[1359]
pnew[1360]
pnew[1361]
pnew[1362]
pnew[1363]
pnew[1364]
pnew[1365]
pnew[1366]
pnew[1367]
pnew[1368]
pnew[1369]
pnew[1370]
pnew[1371]
pnew[1372]
pnew[1373]
pnew[1374]
pnew[1375]
pnew[1376]
pnew[1377]
pnew[1378]
pnew[1379]
pnew[1380]
pnew[1381]
pnew[1382]
pnew[1383]
pnew[1384]
pnew[1385]
pnew[1386]
pnew[1387]
pnew[1388]
pnew[1389]
pnew[1390]
pnew[1391]
pnew[1392]
pnew[1393]
pnew[1394]
pnew[1395]
pnew[1396]
pnew[1397]
pnew[1398]
pnew[1399]
pnew[1400]
pnew[1401]
pnew[1402]
pnew[1403]
pnew[1404]
pnew[1405]
pnew[1406]
pnew[1407]
pnew[1408]
pnew[1409]
pnew[1410]
pnew[1411]
pnew[1412]
pnew[1413]
pnew[1414]
pnew[1415]
pnew[1416]
pnew[1417]
pnew[1418]
pnew[1419]
pnew[1420]
pnew[1421]
pnew[1422]
pnew[1423]
pnew[1424]
pnew[1425]
pnew[1426]
pnew[1427]
pnew[1428]
pnew[1429]
pnew[1430]
pnew[1431]
pnew[1432]
pnew[1433]
pnew[1434]
pnew[1435]
pnew[1436]
pnew[1437]
pnew[1438]
pnew[1439]
pnew[1440]
pnew[1441]
pnew[1442]
pnew[1443]
pnew[1444]
pnew[1445]
pnew[1446]
pnew[1447]
pnew[1448]
pnew[1449]
pnew[1450]
pnew[1451]
pnew[1452]
pnew[1453]
pnew[1454]
pnew[1455]
pnew[1456]
pnew[1457]
pnew[1458]
pnew[1459]
pnew[1460]
pnew[1461]
pnew[1462]
pnew[1463]
pnew[1464]
pnew[1465]
pnew[1466]
pnew[1467]
pnew[1468]
pnew[1469]
pnew[1470]
pnew[1471]
pnew[1472]
pnew[1473]
pnew[1474]
pnew[1475]
pnew[1476]
pnew[1477]
pnew[1478]
pnew[1479]
pnew[1480]
pnew[1481]
pnew[1482]
pnew[1483]
pnew[1484]
pnew[1485]
pnew[1486]
pnew[1487]
pnew[1488]
pnew[1489]
pnew[1490]
pnew[1491]
pnew[1492]
pnew[1493]
pnew[1494]
pnew[1495]
pnew[1496]
pnew[1497]
pnew[1498]
pnew[1499]
pnew[1500]
pnew[1501]
pnew[1502]
pnew[1503]
pnew[1504]
pnew[1505]
pnew[1506]
pnew[1507]
pnew[1508]
pnew[1509]
pnew[1510]
pnew[1511]
pnew[1512]
pnew[1513]
pnew[1514]
pnew[1515]
pnew[1516]
pnew[1517]
pnew[1518]
pnew[1519]
pnew[1520]
pnew[1521]
pnew[1522]
pnew[1523]
pnew[1524]
pnew[1525]
pnew[1526]
pnew[1527]
pnew[1528]
pnew[1529]
pnew[1530]
pnew[1531]
pnew[1532]
pnew[1533]
pnew[1534]
pnew[1535]
pnew[1536]
pnew[1537]
pnew[1538]
pnew[1539]
pnew[1540]
pnew[1541]
pnew[1542]
pnew[1543]
pnew[1544]
pnew[1545]
pnew[1546]
pnew[1547]
pnew[1548]
pnew[1549]
pnew[1550]
pnew[1551]
pnew[1552]
pnew[1553]
pnew[1554]
pnew[1555]
pnew[1556]
pnew[1557]
pnew[1558]
pnew[1559]
pnew[1560]
pnew[1561]
pnew[1562]
pnew[1563]
pnew[1564]
pnew[1565]
pnew[1566]
pnew[1567]
pnew[1568]
pnew[1569]
pnew[1570]
pnew[1571]
pnew[1572]
pnew[1573]
pnew[1574]
pnew[1575]
pnew[1576]
pnew[1577]
pnew[1578]
pnew[1579]
pnew[1580]
pnew[1581]
pnew[1582]
pnew[1583]
pnew[1584]
pnew[1585]
pnew[1586]
pnew[1587]
pnew[1588]
pnew[1589]
pnew[1590]
pnew[1591]
pnew[1592]
pnew[1593]
pnew[1594]
pnew[1595]
pnew[1596]
pnew[1597]
pnew[1598]
pnew[1599]
pnew[1600]
pnew[1601]
pnew[1602]
pnew[1603]
pnew[1604]
pnew[1605]
pnew[1606]
pnew[1607]
pnew[1608]
pnew[1609]
pnew[1610]
pnew[1611]
pnew[1612]
pnew[1613]
pnew[1614]
pnew[1615]
pnew[1616]
pnew[1617]
pnew[1618]
pnew[1619]
pnew[1620]
pnew[1621]
pnew[1622]
pnew[1623]
pnew[1624]
pnew[1625]
pnew[1626]
pnew[1627]
pnew[1628]
pnew[1629]
pnew[1630]
pnew[1631]
pnew[1632]
pnew[1633]
pnew[1634]
pnew[1635]
pnew[1636]
pnew[1637]
pnew[1638]
pnew[1639]
pnew[1640]
pnew[1641]
pnew[1642]
pnew[1643]
pnew[1644]
pnew[1645]
pnew[1646]
pnew[1647]
pnew[1648]
pnew[1649]
pnew[1650]
pnew[1651]
pnew[1652]
pnew[1653]
pnew[1654]
pnew[1655]
pnew[1656]
pnew[1657]
pnew[1658]
pnew[1659]
pnew[1660]
pnew[1661]
pnew[1662]
pnew[1663]
pnew[1664]
pnew[1665]
pnew[1666]
pnew[1667]
pnew[1668]
pnew[1669]
pnew[1670]
pnew[1671]
pnew[1672]
pnew[1673]
pnew[1674]
pnew[1675]
pnew[1676]
pnew[1677]
pnew[1678]
pnew[1679]
pnew[1680]
pnew[1681]
pnew[1682]
pnew[1683]
pnew[1684]
pnew[1685]
pnew[1686]
pnew[1687]
pnew[1688]
pnew[1689]
pnew[1690]
pnew[1691]
pnew[1692]
pnew[1693]
pnew[1694]
pnew[1695]
pnew[1696]
pnew[1697]
pnew[1698]
pnew[1699]
pnew[1700]
pnew[1701]
pnew[1702]
pnew[1703]
pnew[1704]
pnew[1705]
pnew[1706]
pnew[1707]
pnew[1708]
pnew[1709]
pnew[1710]
pnew[1711]
pnew[1712]
pnew[1713]
pnew[1714]
pnew[1715]
pnew[1716]
pnew[1717]
pnew[1718]
pnew[1719]
pnew[1720]
pnew[1721]
pnew[1722]
pnew[1723]
pnew[1724]
pnew[1725]
pnew[1726]
pnew[1727]
pnew[1728]
pnew[1729]
pnew[1730]
pnew[1731]
pnew[1732]
pnew[1733]
pnew[1734]
pnew[1735]
pnew[1736]
pnew[1737]
pnew[1738]
pnew[1739]
pnew[1740]
pnew[1741]
pnew[1742]
pnew[1743]
pnew[1744]
pnew[1745]
pnew[1746]
pnew[1747]
pnew[1748]
pnew[1749]
pnew[1750]
pnew[1751]
pnew[1752]
pnew[1753]
pnew[1754]
pnew[1755]
pnew[1756]
pnew[1757]
pnew[1758]
pnew[1759]
pnew[1760]
pnew[1761]
pnew[1762]
pnew[1763]
pnew[1764]
pnew[1765]
pnew[1766]
pnew[1767]
pnew[1768]
pnew[1769]
pnew[1770]
pnew[1771]
pnew[1772]
pnew[1773]
pnew[1774]
pnew[1775]
pnew[1776]
pnew[1777]
pnew[1778]
pnew[1779]
pnew[1780]
pnew[1781]
pnew[1782]
pnew[1783]
pnew[1784]
pnew[1785]
pnew[1786]
pnew[1787]
pnew[1788]
pnew[1789]
pnew[1790]
pnew[1791]
pnew[1792]
pnew[1793]
pnew[1794]
pnew[1795]
pnew[1796]
pnew[1797]
pnew[1798]
pnew[1799]
pnew[1800]
pnew[1801]
pnew[1802]
pnew[1803]
pnew[1804]
pnew[1805]
pnew[1806]
pnew[1807]
pnew[1808]
pnew[1809]
pnew[1810]
pnew[1811]
pnew[1812]
pnew[1813]
pnew[1814]
pnew[1815]
pnew[1816]
pnew[1817]
pnew[1818]
pnew[1819]
pnew[1820]
pnew[1821]
pnew[1822]
pnew[1823]
pnew[1824]
pnew[1825]
pnew[1826]
pnew[1827]
pnew[1828]
pnew[1829]
pnew[1830]
pnew[1831]
pnew[1832]
pnew[1833]
pnew[1834]
pnew[1835]
pnew[1836]
pnew[1837]
pnew[1838]
pnew[1839]
pnew[1840]
pnew[1841]
pnew[1842]
pnew[1843]
pnew[1844]
pnew[1845]
pnew[1846]
pnew[1847]
pnew[1848]
pnew[1849]
pnew[1850]
pnew[1851]
pnew[1852]
pnew[1853]
pnew[1854]
pnew[1855]
pnew[1856]
pnew[1857]
pnew[1858]
pnew[1859]
pnew[1860]
pnew[1861]
pnew[1862]
pnew[1863]
pnew[1864]
pnew[1865]
pnew[1866]
pnew[1867]
pnew[1868]
pnew[1869]
pnew[1870]
pnew[1871]
pnew[1872]
pnew[1873]
pnew[1874]
pnew[1875]
pnew[1876]
pnew[1877]
pnew[1878]
pnew[1879]
pnew[1880]
pnew[1881]
pnew[1882]
pnew[1883]
pnew[1884]
pnew[1885]
pnew[1886]
pnew[1887]
pnew[1888]
pnew[1889]
pnew[1890]
pnew[1891]
pnew[1892]
pnew[1893]
pnew[1894]
pnew[1895]
pnew[1896]
pnew[1897]
pnew[1898]
pnew[1899]
pnew[1900]
pnew[1901]
pnew[1902]
pnew[1903]
pnew[1904]
pnew[1905]
pnew[1906]
pnew[1907]
pnew[1908]
pnew[1909]
pnew[1910]
pnew[1911]
pnew[1912]
pnew[1913]
pnew[1914]
pnew[1915]
pnew[1916]
pnew[1917]
pnew[1918]
pnew[1919]
pnew[1920]
pnew[1921]
pnew[1922]
pnew[1923]
pnew[1924]
pnew[1925]
pnew[1926]
pnew[1927]
pnew[1928]
pnew[1929]
pnew[1930]
pnew[1931]
pnew[1932]
pnew[1933]
pnew[1934]
pnew[1935]
pnew[1936]
pnew[1937]
pnew[1938]
pnew[1939]
pnew[1940]
pnew[1941]
pnew[1942]
pnew[1943]
pnew[1944]
pnew[1945]
pnew[1946]
pnew[1947]
pnew[1948]
pnew[1949]
pnew[1950]
pnew[1951]
pnew[1952]
pnew[1953]
pnew[1954]
pnew[1955]
pnew[1956]
pnew[1957]
pnew[1958]
pnew[1959]
pnew[1960]
pnew[1961]
pnew[1962]
pnew[1963]
pnew[1964]
pnew[1965]
pnew[1966]
pnew[1967]
pnew[1968]
pnew[1969]
pnew[1970]
pnew[1971]
pnew[1972]
pnew[1973]
pnew[1974]
pnew[1975]
pnew[1976]
pnew[1977]
pnew[1978]
pnew[1979]
pnew[1980]
pnew[1981]
pnew[1982]
pnew[1983]
pnew[1984]
pnew[1985]
pnew[1986]
pnew[1987]
pnew[1988]
pnew[1989]
pnew[1990]
pnew[1991]
pnew[1992]
pnew[1993]
pnew[1994]
pnew[1995]
pnew[1996]
pnew[1997]
pnew[1998]
pnew[1999]
pnew[2000]
pnew[2001]
pnew[2002]
pnew[2003]
pnew[2004]
pnew[2005]
pnew[2006]
pnew[2007]
pnew[2008]
pnew[2009]
pnew[2010]
pnew[2011]
pnew[2012]
pnew[2013]
pnew[2014]
pnew[2015]
pnew[2016]
pnew[2017]
pnew[2018]
pnew[2019]
pnew[2020]
pnew[2021]
pnew[2022]
pnew[2023]
pnew[2024]
pnew[2025]
pnew[2026]
pnew[2027]
pnew[2028]
pnew[2029]
pnew[2030]
pnew[2031]
pnew[2032]
pnew[2033]
pnew[2034]
pnew[2035]
pnew[2036]
pnew[2037]
pnew[2038]
pnew[2039]
pnew[2040]
pnew[2041]
pnew[2042]
pnew[2043]
pnew[2044]
pnew[2045]
pnew[2046]
pnew[2047]
pnew[2048]
pnew[2049]
pnew[2050]
pnew[2051]
pnew[2052]
pnew[2053]
pnew[2054]
pnew[2055]
pnew[2056]
pnew[2057]
pnew[2058]
pnew[2059]
pnew[2060]
pnew[2061]
pnew[2062]
pnew[2063]
pnew[2064]
pnew[2065]
pnew[2066]
pnew[2067]
pnew[2068]
pnew[2069]
pnew[2070]
pnew[2071]
pnew[2072]
pnew[2073]
pnew[2074]
pnew[2075]
pnew[2076]
pnew[2077]
pnew[2078]
pnew[2079]
pnew[2080]
pnew[2081]
pnew[2082]
pnew[2083]
pnew[2084]
pnew[2085]
pnew[2086]
pnew[2087]
pnew[2088]
pnew[2089]
pnew[2090]
pnew[2091]
pnew[2092]
pnew[2093]
pnew[2094]
pnew[2095]
pnew[2096]
pnew[2097]
pnew[2098]
pnew[2099]
pnew[2100]
pnew[2101]
pnew[2102]
pnew[2103]
pnew[2104]
pnew[2105]
pnew[2106]
pnew[2107]
pnew[2108]
pnew[2109]
pnew[2110]
pnew[2111]
pnew[2112]
pnew[2113]
pnew[2114]
pnew[2115]
pnew[2116]
pnew[2117]
pnew[2118]
pnew[2119]
pnew[2120]
pnew[2121]
pnew[2122]
pnew[2123]
pnew[2124]
pnew[2125]
pnew[2126]
pnew[2127]
pnew[2128]
pnew[2129]
pnew[2130]
pnew[2131]
pnew[2132]
pnew[2133]
pnew[2134]
pnew[2135]
pnew[2136]
pnew[2137]
pnew[2138]
pnew[2139]
pnew[2140]
pnew[2141]
pnew[2142]
pnew[2143]
pnew[2144]
pnew[2145]
pnew[2146]
pnew[2147]
pnew[2148]
pnew[2149]
pnew[2150]
pnew[2151]
pnew[2152]
pnew[2153]
pnew[2154]
pnew[2155]
pnew[2156]
pnew[2157]
pnew[2158]
pnew[2159]
pnew[2160]
pnew[2161]
pnew[2162]
pnew[2163]
pnew[2164]
pnew[2165]
pnew[2166]
pnew[2167]
pnew[2168]
pnew[2169]
pnew[2170]
pnew[2171]
pnew[2172]
pnew[2173]
pnew[2174]
pnew[2175]
pnew[2176]
pnew[2177]
pnew[2178]
pnew[2179]
pnew[2180]
pnew[2181]
pnew[2182]
pnew[2183]
pnew[2184]
pnew[2185]
pnew[2186]
pnew[2187]
pnew[2188]
pnew[2189]
pnew[2190]
pnew[2191]
pnew[2192]
pnew[2193]
pnew[2194]
pnew[2195]
pnew[2196]
pnew[2197]
pnew[2198]
pnew[2199]
pnew[2200]
pnew[2201]
pnew[2202]
pnew[2203]
pnew[2204]
pnew[2205]
pnew[2206]
pnew[2207]
pnew[2208]
pnew[2209]
pnew[2210]
pnew[2211]
pnew[2212]
pnew[2213]
pnew[2214]
pnew[2215]
pnew[2216]
pnew[2217]
pnew[2218]
pnew[2219]
pnew[2220]
pnew[2221]
pnew[2222]
pnew[2223]
pnew[2224]
pnew[2225]
pnew[2226]
pnew[2227]
pnew[2228]
pnew[2229]
pnew[2230]
pnew[2231]
pnew[2232]
pnew[2233]
pnew[2234]
pnew[2235]
pnew[2236]
pnew[2237]
pnew[2238]
pnew[2239]
pnew[2240]
pnew[2241]
pnew[2242]
pnew[2243]
pnew[2244]
pnew[2245]
pnew[2246]
pnew[2247]
pnew[2248]
pnew[2249]
pnew[2250]
pnew[2251]
pnew[2252]
pnew[2253]
pnew[2254]
pnew[2255]
pnew[2256]
pnew[2257]
pnew[2258]
pnew[2259]
pnew[2260]
pnew[2261]
pnew[2262]
pnew[2263]
pnew[2264]
pnew[2265]
pnew[2266]
pnew[2267]
pnew[2268]
pnew[2269]
pnew[2270]
pnew[2271]
pnew[2272]
pnew[2273]
pnew[2274]
pnew[2275]
pnew[2276]
pnew[2277]
pnew[2278]
pnew[2279]
pnew[2280]
pnew[2281]
pnew[2282]
pnew[2283]
pnew[2284]
pnew[2285]
pnew[2286]
pnew[2287]
pnew[2288]
pnew[2289]
pnew[2290]
pnew[2291]
pnew[2292]
pnew[2293]
pnew[2294]
pnew[2295]
pnew[2296]
pnew[2297]
pnew[2298]
pnew[2299]
pnew[2300]
pnew[2301]
pnew[2302]
pnew[2303]
pnew[2304]
pnew[2305]
pnew[2306]
pnew[2307]
pnew[2308]
pnew[2309]
pnew[2310]
pnew[2311]
pnew[2312]
pnew[2313]
pnew[2314]
pnew[2315]
pnew[2316]
pnew[2317]
pnew[2318]
pnew[2319]
pnew[2320]
pnew[2321]
pnew[2322]
pnew[2323]
pnew[2324]
pnew[2325]
pnew[2326]
pnew[2327]
pnew[2328]
pnew[2329]
pnew[2330]
pnew[2331]
pnew[2332]
pnew[2333]
pnew[2334]
pnew[2335]
pnew[2336]
pnew[2337]
pnew[2338]
pnew[2339]
pnew[2340]
pnew[2341]
pnew[2342]
pnew[2343]
pnew[2344]
pnew[2345]
pnew[2346]
pnew[2347]
pnew[2348]
pnew[2349]
pnew[2350]
pnew[2351]
pnew[2352]
pnew[2353]
pnew[2354]
pnew[2355]
pnew[2356]
pnew[2357]
pnew[2358]
pnew[2359]
pnew[2360]
pnew[2361]
pnew[2362]
pnew[2363]
pnew[2364]
pnew[2365]
pnew[2366]
pnew[2367]
pnew[2368]
pnew[2369]
pnew[2370]
pnew[2371]
pnew[2372]
pnew[2373]
pnew[2374]
pnew[2375]
pnew[2376]
pnew[2377]
pnew[2378]
pnew[2379]
pnew[2380]
pnew[2381]
pnew[2382]
pnew[2383]
pnew[2384]
pnew[2385]
pnew[2386]
pnew[2387]
pnew[2388]
pnew[2389]
pnew[2390]
pnew[2391]
pnew[2392]
pnew[2393]
pnew[2394]
pnew[2395]
pnew[2396]
pnew[2397]
pnew[2398]
pnew[2399]
pnew[2400]
pnew[2401]
pnew[2402]
pnew[2403]
pnew[2404]
pnew[2405]
pnew[2406]
pnew[2407]
pnew[2408]
pnew[2409]
pnew[2410]
pnew[2411]
pnew[2412]
pnew[2413]
pnew[2414]
pnew[2415]
pnew[2416]
pnew[2417]
pnew[2418]
pnew[2419]
pnew[2420]
pnew[2421]
pnew[2422]
pnew[2423]
pnew[2424]
pnew[2425]
pnew[2426]
pnew[2427]
pnew[2428]
pnew[2429]
pnew[2430]
pnew[2431]
pnew[2432]
pnew[2433]
pnew[2434]
pnew[2435]
pnew[2436]
pnew[2437]
pnew[2438]
pnew[2439]
pnew[2440]
pnew[2441]
pnew[2442]
pnew[2443]
pnew[2444]
pnew[2445]
pnew[2446]
pnew[2447]
pnew[2448]
pnew[2449]
pnew[2450]
pnew[2451]
pnew[2452]
pnew[2453]
pnew[2454]
pnew[2455]
pnew[2456]
pnew[2457]
pnew[2458]
pnew[2459]
pnew[2460]
pnew[2461]
pnew[2462]
pnew[2463]
pnew[2464]
pnew[2465]
pnew[2466]
pnew[2467]
pnew[2468]
pnew[2469]
pnew[2470]
pnew[2471]
pnew[2472]
pnew[2473]
pnew[2474]
pnew[2475]
pnew[2476]
pnew[2477]
pnew[2478]
pnew[2479]
pnew[2480]
pnew[2481]
pnew[2482]
pnew[2483]
pnew[2484]
pnew[2485]
pnew[2486]
pnew[2487]
pnew[2488]
pnew[2489]
pnew[2490]
pnew[2491]
pnew[2492]
pnew[2493]
pnew[2494]
pnew[2495]
pnew[2496]
pnew[2497]
pnew[2498]
pnew[2499]

In [6]:
print results.shape


(2501, 3)

In [7]:
print results[0,:]


['pnew' 'redshift' 'stellar_mass']

In [8]:
my_results = {}
for i in range(len(results[0, :])):                             # Converting numpy array into dictionary
    my_results[results[0, i]] = np.array(results[0 + 1:, i], dtype=str)

Analyzing the dataset


In [9]:
print new_probability.shape
print new_probability[0,:]


(5509, 11)
['parameter' 'mean' 'se_mean' 'sd' '2.5%' '25%' '50%' '75%' '97.5%' 'n_eff'
 'Rhat']

In [10]:
parameter    = my_probability['parameter'].astype(str)
quartile_1st = my_probability['25%'].astype(float)
quartile_3rd = my_probability['75%'].astype(float)
quantile_25  = my_probability['2.5%'].astype(float)
quantile_97  = my_probability['97.5%'].astype(float)
mean         = my_probability['mean'].astype(float)

In [11]:
print np.unique(parameter).size


5508

In [12]:
probability   = my_results['pnew'].astype(float)
redshift      = my_results['redshift'].astype(float)
stellar_mass  = my_results['stellar_mass'].astype(float)

In [13]:
print mean.size
print probability.size
print redshift.size
print stellar_mass.size


5508
2500
2500
2500

In [14]:
idx_beg = []
idx_end = []
for i in range(parameter.size):
    if (parameter[i]=='pnew[0]'):
        idx_beg = int(i)
        print idx_beg
        print 'ok'
    elif (parameter[i]=='pnew[2499]'):
        idx_end = int(i)
        print idx_end
        print 'ok'
    else:
        continue


3008
ok
5507
ok

In [15]:
new_q1 = quartile_1st[idx_beg:idx_end+1]
new_q3 = quartile_3rd[idx_beg:idx_end+1]
new_25 = quantile_25[idx_beg:idx_end+1]
new_97 = quantile_97[idx_beg:idx_end+1]
new_param = parameter[idx_beg:idx_end+1]
new_mean = mean[idx_beg:idx_end+1]

In [16]:
print new_mean.size


2500

In [17]:
print redshift.size


2500

In [18]:
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
# plt.scatter(redshift, sersic_gal, mean, c = '#feb24c', alpha=0.7)
ax.plot_trisurf(redshift, stellar_mass, new_mean, color='#a6611a', alpha=0.8, linewidth=0, antialiased=False)
ax.set_xlabel("\n z", fontsize=14)
ax.set_ylabel("\n Log M$_*$ (M$_\odot$)", fontsize=14)
ax.set_zlabel("\n Probability of \n UV upturn", fontsize=14)
for t in ax.yaxis.get_major_ticks(): t.label.set_fontsize(10)
ax.set_yticks([9.75, 10.25, 10.75, 11.25, 11.75])
ax.view_init(elev=22., azim=315)
# ax.zaxis.set_rotate_label(False)
plt.savefig('./../../Figs/logit3D_grid.pdf', dpi=100)

# for ii in xrange(180,360,10):
#     ax.view_init(elev=10., azim=ii)
#     plt.savefig("./Results/movie_short/movie_0%d_0.png" % ii)
    
# plt.savefig("./movie_0%d.png" % 359)

plt.show()



In [19]:
plt.scatter(redshift, new_mean, alpha=0.2)
plt.xlabel("Redshift")
plt.ylabel("Probability")
plt.savefig("./Results/dependency01.png")
plt.show()



In [20]:
plt.scatter(stellar_mass, new_mean, alpha=0.2)
plt.xlabel("Stellar Mass")
plt.ylabel("Probability")
plt.savefig("./Results/dependency02.png")
plt.show()



In [21]:
print np.unique(redshift).size
print redshift.size
print np.unique(redshift)


50
2500
[ 0.06794     0.07366163  0.07938327  0.0851049   0.09082653  0.09654816
  0.1022698   0.10799143  0.11371306  0.11943469  0.12515633  0.13087796
  0.13659959  0.14232122  0.14804286  0.15376449  0.15948612  0.16520776
  0.17092939  0.17665102  0.18237265  0.18809429  0.19381592  0.19953755
  0.20525918  0.21098082  0.21670245  0.22242408  0.22814571  0.23386735
  0.23958898  0.24531061  0.25103224  0.25675388  0.26247551  0.26819714
  0.27391878  0.27964041  0.28536204  0.29108367  0.29680531  0.30252694
  0.30824857  0.3139702   0.31969184  0.32541347  0.3311351   0.33685673
  0.34257837  0.3483    ]

Let's find in which indices the redshifts are the same, so we can estimate medians (of prob., 1st quartile, etc.)


In [22]:
unique_redshifts = np.unique(redshift)

In [23]:
idxs_unique = []
matrix = []
for i in range(unique_redshifts.size):
    idxs_unique.append("z_%d" % i)
    idxs_unique.append(unique_redshifts[i])
    for j in range(redshift.size):
        if unique_redshifts[i] == redshift[j]:
            idxs_unique.append(j)
        else:
            continue
new_dataset =  np.array(idxs_unique).reshape(50,52).T

In [24]:
print new_dataset


[['z_0' 'z_1' 'z_2' ..., 'z_47' 'z_48' 'z_49']
 ['0.06794' '0.0736616326531' '0.0793832653061' ..., '0.336856734694'
  '0.342578367347' '0.3483']
 ['0' '1' '2' ..., '47' '48' '49']
 ..., 
 ['2350' '2351' '2352' ..., '2397' '2398' '2399']
 ['2400' '2401' '2402' ..., '2447' '2448' '2449']
 ['2450' '2451' '2452' ..., '2497' '2498' '2499']]

In [25]:
unique_mass = np.unique(stellar_mass)
print unique_mass


[  9.729482     9.76925931   9.80903661   9.84881392   9.88859122
   9.92836853   9.96814584  10.00792314  10.04770045  10.08747776
  10.12725506  10.16703237  10.20680967  10.24658698  10.28636429
  10.32614159  10.3659189   10.4056962   10.44547351  10.48525082
  10.52502812  10.56480543  10.60458273  10.64436004  10.68413735
  10.72391465  10.76369196  10.80346927  10.84324657  10.88302388
  10.92280118  10.96257849  11.0023558   11.0421331   11.08191041
  11.12168771  11.16146502  11.20124233  11.24101963  11.28079694
  11.32057424  11.36035155  11.40012886  11.43990616  11.47968347
  11.51946078  11.55923808  11.59901539  11.63879269  11.67857   ]

In [26]:
print new_dataset[:,0].size


52

In [27]:
mean_prob = []
prob_i = []
q1_i   = []
q3_i   = []
q97_i  = []
q25_i  = []
matrix = []

for k in range(new_dataset[2:,0].size):
    column   = new_dataset[:,k]
    idxs_col = column[2:]
    mean_prob = list(column[1:2])
#     print mean_prob
    for m in idxs_col:
        prob_i.append(new_mean[int(m)])
        q1_i.append(new_q1[int(m)])
        q3_i.append(new_q3[int(m)])
        q97_i.append(new_97[int(m)])
        q25_i.append(new_25[int(m)])
    prob = np.mean(np.array(prob_i))
    q1   = np.mean(np.array(q1_i))
    q3   = np.mean(np.array(q3_i))
    q97  = np.mean(np.array(q97_i))
    q25  = np.mean(np.array(q25_i))
    mean_prob.append(prob)
    mean_prob.append(q1)
    mean_prob.append(q3)
    mean_prob.append(q25)
    mean_prob.append(q97)
    matrix.append(mean_prob)

In [28]:
matrix = np.array(matrix)

In [29]:
print matrix


[['0.06794' '0.211242' '0.152192' '0.259632' '0.087704' '0.390316']
 ['0.0736616326531' '0.217716' '0.15948' '0.265868' '0.094138' '0.393526']
 ['0.0793832653061' '0.224304' '0.166928' '0.272176' '0.100824' '0.396978']
 ['0.0851048979592' '0.2309925' '0.174504' '0.2785485' '0.107767'
  '0.400537']
 ['0.0908265306122' '0.237766' '0.1821728' '0.2849876' '0.1149044'
  '0.4042716']
 ['0.0965481632653' '0.244606' '0.189893333333' '0.291455333333' '0.122215'
  '0.408158666667']
 ['0.102269795918' '0.251494' '0.197649714286' '0.29796' '0.129661142857'
  '0.412159142857']
 ['0.107991428571' '0.25841075' '0.205428' '0.3045025' '0.137192'
  '0.41625775']
 ['0.113713061224' '0.265338444444' '0.213185333333' '0.311051333333'
  '0.144749333333' '0.420517555556']
 ['0.119434693878' '0.272257' '0.2208938' '0.3175894' '0.1522854' '0.42498']
 ['0.125156326531' '0.279149636364' '0.228527090909' '0.324105636364'
  '0.159761818182' '0.429627636364']
 ['0.130877959184' '0.2859955' '0.236041666667' '0.330598666667'
  '0.167121833333' '0.434415833333']
 ['0.136599591837' '0.292778923077' '0.243415384615' '0.337069846154'
  '0.174322' '0.439352615385']
 ['0.14232122449' '0.299482714286' '0.250633857143' '0.343506285714'
  '0.181342142857' '0.444416285714']
 ['0.148042857143' '0.306092133333' '0.257682133333' '0.349894133333'
  '0.1881696' '0.449614933333']
 ['0.153764489796' '0.31259225' '0.264554625' '0.3562245' '0.194752375'
  '0.454917']
 ['0.159486122449' '0.318969882353' '0.271247058824' '0.362486'
  '0.201098235294' '0.460284']
 ['0.165207755102' '0.325213111111' '0.277746888889' '0.368663333333'
  '0.207213666667' '0.465676444444']
 ['0.170929387755' '0.331311263158' '0.284051263158' '0.374740947368'
  '0.213099894737' '0.471074105263']
 ['0.176651020408' '0.3372546' '0.2901486' '0.3807002' '0.2187523'
  '0.4764495']
 ['0.182372653061' '0.343034285714' '0.296037142857' '0.38652752381'
  '0.224176380952' '0.481781333333']
 ['0.188094285714' '0.348642272727' '0.301716545455' '0.392213'
  '0.229379818182' '0.487040909091']
 ['0.193815918367' '0.354072608696' '0.307185217391' '0.397744695652'
  '0.234367826087' '0.492206434783']
 ['0.19953755102' '0.359319833333' '0.312442666667' '0.403113666667'
  '0.23914325' '0.497268']
 ['0.205259183673' '0.364378' '0.31748928' '0.40831456' '0.24370728'
  '0.50220088']
 ['0.210980816327' '0.369243153846' '0.322322230769' '0.413337923077'
  '0.248057153846' '0.506989769231']
 ['0.21670244898' '0.373911703704' '0.326939777778' '0.418176296296'
  '0.252188296296' '0.511627777778']
 ['0.222424081633' '0.378379928571' '0.3313415' '0.422824071429'
  '0.256106285714' '0.516107428571']
 ['0.228145714286' '0.382645862069' '0.335526896552' '0.427280206897'
  '0.259813931034' '0.520425517241']
 ['0.233867346939' '0.386706466667' '0.339493466667' '0.4315402'
  '0.263295133333' '0.524582666667']
 ['0.239588979592' '0.390559741935' '0.343234580645' '0.435598064516'
  '0.266545419355' '0.52858683871']
 ['0.245310612245' '0.3942038125' '0.346745875' '0.43945825' '0.2695650625'
  '0.5324446875']
 ['0.251032244898' '0.397637575758' '0.350022181818' '0.443126666667'
  '0.272339454545' '0.536165818182']
 ['0.256753877551' '0.400860117647' '0.353059235294' '0.446604529412'
  '0.274855352941' '0.539745411765']
 ['0.262475510204' '0.403869714286' '0.355854914286' '0.449891371429'
  '0.277104514286' '0.543191885714']
 ['0.268197142857' '0.406665666667' '0.358403888889' '0.452987166667'
  '0.279083611111' '0.546514388889']
 ['0.27391877551' '0.409247243243' '0.360702756757' '0.455897189189'
  '0.280775027027' '0.549723891892']
 ['0.279640408163' '0.411614052632' '0.362745263158' '0.458623526316'
  '0.282160736842' '0.552828263158']
 ['0.285362040816' '0.413765794872' '0.364527333333' '0.461170717949'
  '0.283236' '0.555836923077']
 ['0.291083673469' '0.4157024' '0.3660485' '0.46354575' '0.283998'
  '0.55876205']
 ['0.296805306122' '0.417423902439' '0.367304585366' '0.465752878049'
  '0.284436878049' '0.561611219512']
 ['0.302526938776' '0.418931142857' '0.368295333333' '0.467793238095'
  '0.284556' '0.564389952381']
 ['0.308248571429' '0.420224697674' '0.369014046512' '0.46966827907'
  '0.284362604651' '0.567106837209']
 ['0.313970204082' '0.421305818182' '0.369462954545' '0.471379318182'
  '0.283862272727' '0.569772']
 ['0.319691836735' '0.422176177778' '0.369639422222' '0.472928977778'
  '0.283055955556' '0.5723912']
 ['0.325413469388' '0.422837913043' '0.369545869565' '0.474323782609'
  '0.281948913043' '0.574966130435']
 ['0.331135102041' '0.423293659574' '0.369184765957' '0.475561191489'
  '0.280552425532' '0.577503574468']
 ['0.336856734694' '0.423546375' '0.368558125' '0.476641166667'
  '0.278878416667' '0.580007583333']
 ['0.342578367347' '0.423600122449' '0.36766877551' '0.477567428571'
  '0.276945306122' '0.582483673469']
 ['0.3483' '0.423459' '0.36652168' '0.47833868' '0.27477096' '0.58492932']]

Let's do a similar thing for the mass

Beeing completely over cautious


In [30]:
idxs_unique_m = []
matrix_m = []
for i in range(unique_mass.size):
    idxs_unique_m.append("m_%d" % i)
    idxs_unique_m.append(unique_mass[i])
    for j in range(stellar_mass.size):
        if unique_mass[i] == stellar_mass[j]:
            idxs_unique_m.append(j)
        else:
            continue
new_dataset_m = np.array(idxs_unique_m).reshape(50,52).T

In [31]:
print new_dataset_m


[['m_0' 'm_1' 'm_2' ..., 'm_47' 'm_48' 'm_49']
 ['9.729482' '9.76925930612' '9.80903661224' ..., '11.5990153878'
  '11.6387926939' '11.67857']
 ['0' '50' '100' ..., '2350' '2400' '2450']
 ..., 
 ['47' '97' '147' ..., '2397' '2447' '2497']
 ['48' '98' '148' ..., '2398' '2448' '2498']
 ['49' '99' '149' ..., '2399' '2449' '2499']]

In [32]:
mean_prob_m = []
prob_i_m = []
q1_i_m   = []
q3_i_m   = []
q97_i_m  = []
q25_i_m  = []
matrix_m = []

for k in range(new_dataset_m[2:,0].size):
    column      = new_dataset_m[:,k]
    idxs_col    = column[2:]
    mean_prob_m = list(column[1:2])
#     print mean_prob
    for m in idxs_col:
        prob_i_m.append(new_mean[int(m)])
        q1_i_m.append(new_q1[int(m)])
        q3_i_m.append(new_q3[int(m)])
        q97_i_m.append(new_97[int(m)])
        q25_i_m.append(new_25[int(m)])
    prob = np.mean(np.array(prob_i_m))
    q1   = np.mean(np.array(q1_i_m))
    q3   = np.mean(np.array(q3_i_m))
    q97  = np.mean(np.array(q97_i_m))
    q25  = np.mean(np.array(q25_i_m))
    mean_prob_m.append(prob)
    mean_prob_m.append(q1)
    mean_prob_m.append(q3)
    mean_prob_m.append(q25)
    mean_prob_m.append(q97)
    matrix_m.append(mean_prob_m)

In [33]:
matrix_m = np.array(matrix_m)

In [34]:
print matrix_m


[['9.729482' '0.383428' '0.293846' '0.468616' '0.165238' '0.63433']
 ['9.76925930612' '0.38409' '0.295721' '0.4682' '0.16807' '0.631701']
 ['9.80903661224' '0.384756666667' '0.297605333333' '0.467822666667'
  '0.170953333333' '0.629067333333']
 ['9.84881391837' '0.385432' '0.29948' '0.467464' '0.1738775' '0.626472']
 ['9.88859122449' '0.3861144' '0.3013604' '0.4670956' '0.176832' '0.623874']
 ['9.92836853061' '0.386804' '0.303246' '0.466733333333' '0.179830666667'
  '0.621308666667']
 ['9.96814583673' '0.387502857143' '0.305136857143' '0.466373714286'
  '0.182862571429' '0.618751714286']
 ['10.0079231429' '0.38820875' '0.30702925' '0.4660245' '0.18591775'
  '0.6162025']
 ['10.047700449' '0.388922666667' '0.308919333333' '0.465679111111'
  '0.189001111111' '0.613675777778']
 ['10.0874777551' '0.3896436' '0.3108108' '0.4653496' '0.1921144'
  '0.6111764']
 ['10.1272550612' '0.390371818182' '0.312697636364' '0.465028727273'
  '0.195241090909' '0.608701636364']
 ['10.1670323673' '0.3911075' '0.314579' '0.464721166667' '0.198387166667'
  '0.606240833333']
 ['10.2068096735' '0.391850923077' '0.316454307692' '0.464426'
  '0.201547384615' '0.603806307692']
 ['10.2465869796' '0.392601571429' '0.318324142857' '0.464147'
  '0.204710428571' '0.601403714286']
 ['10.2863642857' '0.3933596' '0.3201844' '0.463884266667' '0.207882'
  '0.599039333333']
 ['10.3261415918' '0.394125125' '0.322036875' '0.463641625' '0.211052'
  '0.596725625']
 ['10.365918898' '0.394897411765' '0.323878470588' '0.463418235294'
  '0.214218352941' '0.59446']
 ['10.4056962041' '0.395677222222' '0.325710333333' '0.463214666667'
  '0.217375666667' '0.592261555556']
 ['10.4454735102' '0.396462947368' '0.327529263158' '0.463031789474'
  '0.220518947368' '0.590129263158']
 ['10.4852508163' '0.3972561' '0.3293307' '0.4628754' '0.2236378'
  '0.5880689']
 ['10.5250281224' '0.398056285714' '0.331114952381' '0.46274352381'
  '0.226722190476' '0.58608552381']
 ['10.5648054286' '0.398862545455' '0.332877909091' '0.462639545455'
  '0.229768181818' '0.584181']
 ['10.6045827347' '0.399675391304' '0.334616695652' '0.462568347826'
  '0.232770608696' '0.582370782609']
 ['10.6443600408' '0.400494666667' '0.336331916667' '0.462532833333'
  '0.235717333333' '0.580660166667']
 ['10.6841373469' '0.40132008' '0.33801968' '0.46253384' '0.2385984'
  '0.57905296']
 ['10.7239146531' '0.402151692308' '0.339676461538' '0.462574846154'
  '0.241413923077' '0.577563769231']
 ['10.7636919592' '0.402989111111' '0.341298222222' '0.462658296296'
  '0.244144888889' '0.57620437037']
 ['10.8034692653' '0.403832428571' '0.342884214286' '0.462785'
  '0.246785214286' '0.574977928571']
 ['10.8432465714' '0.404681241379' '0.344432896552' '0.462959034483'
  '0.249326206897' '0.573884896552']
 ['10.8830238776' '0.4055354' '0.3459392' '0.463183' '0.251760066667'
  '0.572932733333']
 ['10.9228011837' '0.406394967742' '0.347403483871' '0.463457935484'
  '0.254081483871' '0.572120774194']
 ['10.9625784898' '0.407259375' '0.348824375' '0.463784' '0.2562823125'
  '0.5714595625']
 ['11.0023557959' '0.408128787879' '0.350201636364' '0.464162'
  '0.258362242424' '0.570947757576']
 ['11.042133102' '0.409003' '0.351532941176' '0.464592' '0.260318294118'
  '0.570590058824']
 ['11.0819104082' '0.4098816' '0.352817942857' '0.465076857143'
  '0.262146228571' '0.570394285714']
 ['11.1216877143' '0.410764555556' '0.354054722222' '0.465614388889'
  '0.263843277778' '0.570357777778']
 ['11.1614650204' '0.411651621622' '0.355243297297' '0.466204648649'
  '0.265412432432' '0.570477135135']
 ['11.2012423265' '0.412542421053' '0.356381842105' '0.466848368421'
  '0.266852578947' '0.570754368421']
 ['11.2410196327' '0.413437076923' '0.357469692308' '0.467545897436'
  '0.268159128205' '0.571185948718']
 ['11.2807969388' '0.41433515' '0.3585086' '0.4682957' '0.2693367'
  '0.57176865']
 ['11.3205742449' '0.415236439024' '0.359500243902' '0.469094536585'
  '0.270386682927' '0.572502195122']
 ['11.360351551' '0.416141047619' '0.360445285714' '0.469942238095'
  '0.271313571429' '0.573381333333']
 ['11.4001288571' '0.417048186047' '0.361346837209' '0.470837674419'
  '0.27212055814' '0.574400325581']
 ['11.4399061633' '0.417957909091' '0.362204681818' '0.471780454545'
  '0.272815045455' '0.575550727273']
 ['11.4796834694' '0.418869955556' '0.363020977778' '0.472770044444'
  '0.273398888889' '0.576829866667']
 ['11.5194607755' '0.419784347826' '0.363798217391' '0.473803826087'
  '0.273874652174' '0.578232695652']
 ['11.5592380816' '0.420700638298' '0.364536042553' '0.474878'
  '0.274245361702' '0.579750638298']
 ['11.5990153878' '0.421618625' '0.365235333333' '0.475991833333'
  '0.274515' '0.5813755']
 ['11.6387926939' '0.422538122449' '0.365896653061' '0.477146285714'
  '0.274687632653' '0.583102693878']
 ['11.67857' '0.423459' '0.36652168' '0.47833868' '0.27477096' '0.58492932']]

In [35]:
print new_dataset_m[1]
print unique_mass


['9.729482' '9.76925930612' '9.80903661224' '9.84881391837' '9.88859122449'
 '9.92836853061' '9.96814583673' '10.0079231429' '10.047700449'
 '10.0874777551' '10.1272550612' '10.1670323673' '10.2068096735'
 '10.2465869796' '10.2863642857' '10.3261415918' '10.365918898'
 '10.4056962041' '10.4454735102' '10.4852508163' '10.5250281224'
 '10.5648054286' '10.6045827347' '10.6443600408' '10.6841373469'
 '10.7239146531' '10.7636919592' '10.8034692653' '10.8432465714'
 '10.8830238776' '10.9228011837' '10.9625784898' '11.0023557959'
 '11.042133102' '11.0819104082' '11.1216877143' '11.1614650204'
 '11.2012423265' '11.2410196327' '11.2807969388' '11.3205742449'
 '11.360351551' '11.4001288571' '11.4399061633' '11.4796834694'
 '11.5194607755' '11.5592380816' '11.5990153878' '11.6387926939' '11.67857']
[  9.729482     9.76925931   9.80903661   9.84881392   9.88859122
   9.92836853   9.96814584  10.00792314  10.04770045  10.08747776
  10.12725506  10.16703237  10.20680967  10.24658698  10.28636429
  10.32614159  10.3659189   10.4056962   10.44547351  10.48525082
  10.52502812  10.56480543  10.60458273  10.64436004  10.68413735
  10.72391465  10.76369196  10.80346927  10.84324657  10.88302388
  10.92280118  10.96257849  11.0023558   11.0421331   11.08191041
  11.12168771  11.16146502  11.20124233  11.24101963  11.28079694
  11.32057424  11.36035155  11.40012886  11.43990616  11.47968347
  11.51946078  11.55923808  11.59901539  11.63879269  11.67857   ]

In [36]:
remodled_data_z = {}
remodled_data_z['redshift']         = np.array(matrix).T[0]
remodled_data_z['new_probability']  = np.array(matrix).T[1]
remodled_data_z['new_1st_quartile'] = np.array(matrix).T[2]
remodled_data_z['new_3rd_quartile'] = np.array(matrix).T[3]
remodled_data_z['new_2.5_quantile'] = np.array(matrix).T[4]
remodled_data_z['new_97_quantile']  = np.array(matrix).T[5]

In [37]:
remodled_data_m = {}
remodled_data_m['mass']             = np.array(matrix_m).T[0]
remodled_data_m['new_probability']  = np.array(matrix_m).T[1]
remodled_data_m['new_1st_quartile'] = np.array(matrix_m).T[2]
remodled_data_m['new_3rd_quartile'] = np.array(matrix_m).T[3]
remodled_data_m['new_2.5_quantile'] = np.array(matrix_m).T[4]
remodled_data_m['new_97_quantile']  = np.array(matrix_m).T[5]

In [38]:
x   = remodled_data_z['redshift'].astype('float')
y0  = remodled_data_z['new_1st_quartile'].astype('float')
y1  = remodled_data_z['new_probability'].astype('float')
y2  = remodled_data_z['new_3rd_quartile'].astype('float')
y25 = remodled_data_z['new_2.5_quantile'].astype('float')
y97 = remodled_data_z['new_97_quantile'] .astype('float')

In [39]:
palette = ['#1b9e77','#d95f02','#7570b3']

In [40]:
sns.set_style("white")
plt.rcParams["axes.edgecolor"] = "0.15"
plt.rcParams["axes.linewidth"] = 1.

plt.subplots(1,1, figsize=(8,5.1))
plot25, = plt.plot(x, y25, '-.', color=palette[1], label="2.5%")
plotq1, = plt.plot(x, y0,  ':',  color=palette[1], label="25%")
plotq2, = plt.plot(x, y1,  '-',  color=palette[2], label='Median')
plotq3, = plt.plot(x, y2,  ':',  color=palette[1], label="75%")
plot97, = plt.plot(x, y97, '-.', color=palette[1], label="97.5%")

# l1 = plt.legend([plot25, plotq1], ["2.5%", "1$^{\mathrm{st}}$ quartile (25%)"], 
#                 numpoints=1, loc='lower center', fontsize=12, frameon=True, framealpha=0.85)
# l2 = plt.legend([plotq2, plotq3, plot97], ["Median", "3$^{\mathrm{rd}}$ quartile (75%)", "97.5%"], numpoints=1, 
#                 loc='lower right', fontsize=12, frameon=True, framealpha=0.85)
# l1.get_frame().set_edgecolor('black')
# l2.get_frame().set_edgecolor('black')
# plt.gca().add_artist(l1)
plt.fill_between(x, y0, y2,   facecolor=palette[0], alpha=0.1) # credible intervals
plt.fill_between(x, y25, y97, facecolor=palette[0], alpha=0.1) # credible intervals
plt.ylim([0,0.7])
plt.xlabel("Redshift", fontsize=15)
plt.ylabel("Probability of UV upturn", fontsize=15)
plt.legend(numpoints=1, loc='lower right', fontsize=12, frameon=True, framealpha=1.)
plt.tick_params('both', labelsize='14')
plt.savefig('./Results/prob_z.png')
plt.savefig('./Results/prob_z.pdf')
plt.show()



In [41]:
z    = remodled_data_m['mass'].astype('float')
y0m  = remodled_data_m['new_1st_quartile'].astype('float')
y1m  = remodled_data_m['new_probability'].astype('float')
y2m  = remodled_data_m['new_3rd_quartile'].astype('float')
y25m = remodled_data_m['new_2.5_quantile'].astype('float')
y97m = remodled_data_m['new_97_quantile'] .astype('float')

In [42]:
sns.set_style("white")
plt.rcParams["axes.edgecolor"] = "0.15"
plt.rcParams["axes.linewidth"] = 1.

plt.subplots(1,1, figsize=(8,5.1))
plot25, = plt.plot(z, y25m, '-.', color=palette[1], label="2.5%")
plotq1, = plt.plot(z, y0m,  ':',  color=palette[1], label="25%")
plotq2, = plt.plot(z, y1m,  '-',  color=palette[2], label='Median')
plotq3, = plt.plot(z, y2m,  ':',  color=palette[1], label="75%")
plot97, = plt.plot(z, y97m, '-.', color=palette[1], label="97.5%")
plt.fill_between(z, y0m, y2m,   facecolor=palette[0], alpha=0.1) # credible intervals
plt.fill_between(z, y25m, y97m, facecolor=palette[0], alpha=0.1) # credible intervals
plt.xlabel("$\log$ M$_*$ (M$_\odot$)", fontsize=15)
plt.ylabel("Probability of UV upturn", fontsize=15)
plt.legend(numpoints=1, loc='lower right', fontsize=12, frameon=True, framealpha=1)
plt.tick_params('both', labelsize='14')
plt.ylim([0., 0.7])
plt.savefig('./Results/prob_mass.png')
plt.savefig('./Results/prob_mass.pdf')
plt.show()



In [ ]:
x1 = unique_mass
x2 = unique_redshifts
plt.plot(x, z, 'o', alpha=0.8)
plt.xlabel("Redshift")
plt.ylabel("Mass")
plt.savefig('./Results/mass_z.png')
plt.show()

In [ ]: