MALAPP Age Depth Modeling and Pollen Diagrams
Benjamin Davies, David K. Wright, Sarah J. Ivory, Jago J. Birk, Jeong-Heon Choi, Sabine Fiedler, Jacob Davis, Potiphar Kaliba, Jessica C. Thompson
This code is used to generate age-depth models and pollen diagrams for a forthcoming publication:
David K. Wright, Sarah J. Ivory, Jago J. Birk, Jeong-Heon Choi, Benjamin Davies, Sabine Fiedler, Jacob Davis, Potiphar Kaliba, Jessica C. Thompson.”Palaeoenvironmental data track Late Quaternary anthropogenic impacts on vegetation and landscapes in Mzimba, northern Malawi” Submitted to Frontiers in Environmental Archaeology, May 2023.
This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.
Load packages
rm(list = ls())
library(rbacon)Loading required package: IntCallibrary(rioja)This is rioja 1.0-5AU3: Generate age depth model and assign pollen depths
Bacon(core="AU3",d.min=0,d.max=67,accept.suggestions = TRUE,verbose=FALSE)The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 49.0+-964000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69826: 49.0+-29.0 d=14 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69825: 267.0+-29.0 d=19 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69824: 964.0+-29.0 d=28 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: NA: 964.0+-964000000.0 d=67 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=15, H=0, dim=17, Seed=161655, Dc=5.000000, c(0)=0.000000, c(K)=75.000000
twalk: 3740000 iterations to run, Sat Jun 3 07:54:09 2023
twalk thinning: 1 out of every 85 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.out
twalk: Finished, 2.2% of moved pars per iteration (ratio 83855.705883/3740000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU3/AU3_15.out,
Sat Jun 3 07:54:13 2023
bacon: 184049 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -96.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 17000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
|= | 1%
|
|== | 3%
|
|=== | 4%
|
|==== | 6%
|
|===== | 7%
|
|====== | 9%
|
|======= | 10%
|
|======== | 12%
|
|========= | 13%
|
|========== | 15%
|
|=========== | 16%
|
|============ | 18%
|
|============= | 19%
|
|============== | 21%
|
|=============== | 22%
|
|================ | 24%
|
|================= | 25%
|
|================== | 27%
|
|=================== | 28%
|
|==================== | 30%
|
|===================== | 31%
|
|====================== | 33%
|
|======================= | 34%
|
|======================== | 36%
|
|========================= | 37%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 42%
|
|============================= | 43%
|
|============================== | 45%
|
|=============================== | 46%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 51%
|
|=================================== | 52%
|
|==================================== | 54%
|
|===================================== | 55%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|================================================ | 72%
|
|================================================= | 73%
|
|================================================== | 75%
|
|=================================================== | 76%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 84%
|
|========================================================= | 85%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 90%
|
|============================================================= | 91%
|
|============================================================== | 93%
|
|=============================================================== | 94%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2 
core<-rep("AU3",9)
pollenNo<-c("POL11","POL10","POL9","POL8","POL12","POL13","POL16","POL14","POL15")
depth<-c(3,14,19,28,34,40,46,55,67)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG3data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG3data)write.csv(AUG3data,"AUG3data.csv",row.names = FALSE)AU7: Generate age depth model and assign pollen depths
Bacon(core="AU7",d.min=0,d.max=77,accept.suggestions = TRUE,verbose=FALSE)The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 149.0+-2388000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69829: 149.0+-29.0 d=15 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69526: 792.0+-27.0 d=32 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69527: 1233.0+-28.0 d=51 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69528: 927.0+-28.0 d=54 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69529: 2388.0+-37.0 d=77 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=17, H=0, dim=19, Seed=814200, Dc=5.000000, c(0)=0.000000, c(K)=85.000000
twalk: 4180000 iterations to run, Sat Jun 3 07:54:25 2023
twalk thinning: 1 out of every 95 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.out
twalk: Finished, 1.2% of moved pars per iteration (ratio 49705.526316/4180000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU7/AU7_17.out,
Sat Jun 3 07:54:31 2023
bacon: 185239 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -73.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 19000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
|= | 1%
|
|== | 3%
|
|=== | 4%
|
|=== | 5%
|
|==== | 6%
|
|===== | 8%
|
|====== | 9%
|
|======= | 10%
|
|======== | 12%
|
|========= | 13%
|
|========== | 14%
|
|========== | 16%
|
|=========== | 17%
|
|============ | 18%
|
|============= | 19%
|
|============== | 21%
|
|=============== | 22%
|
|================ | 23%
|
|================= | 25%
|
|================= | 26%
|
|================== | 27%
|
|=================== | 29%
|
|==================== | 30%
|
|===================== | 31%
|
|====================== | 32%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 36%
|
|========================= | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 42%
|
|============================= | 43%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 51%
|
|=================================== | 52%
|
|==================================== | 53%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================== | 62%
|
|=========================================== | 64%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 68%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|================================================ | 71%
|
|================================================= | 73%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 83%
|
|========================================================= | 84%
|
|========================================================= | 86%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 90%
|
|============================================================= | 91%
|
|============================================================== | 92%
|
|=============================================================== | 94%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2 
core<-rep("AU7",9)
pollenNo<-c("POL36","POL37","POL38","POL39","POL40","POL41","POL42","POL43","POL44")
depth<-c(6,15,26,32,42,51,54,66,77)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG7data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG7data)write.csv(AUG7data,"AUG7data.csv",row.names = FALSE)AU8: Generate age depth model and assign pollen depths
#Run age-depth model between 0 and max depth
Bacon(core="AU8",d.min=0,d.max=146,depth.unit = "cm",accept.suggestions = TRUE,verbose=FALSE)The run's files will be put in this folder: C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8
Reading C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.bacon
Constant calibration curve.
IntCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_intcal20.14C
Marine20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_marine20.14C
SHCal20: Reading from file: C:\Users\bdav_\AppData\Local\R\win-library\4.2\IntCal\extdata/3Col_shcal20.14C
Added det: NA: 1794.0+-15742000.0 d=0 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69524: 1794.0+-32.0 d=11 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69523: 7759.0+-124.0 d=59 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69828: 9690.0+-51.0 d=74 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: Ua-69522: 15742.0+-59.0 d=86 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
Added det: NA: 15742.0+-15742000000.0 d=146 ResCorr= 0.0+-0.0 a=3 b=4 cc=SHCal20
BaconFixed: Bacon jumps model with fixed c's.
K=31, H=0, dim=33, Seed=970277, Dc=5.000000, c(0)=0.000000, c(K)=155.000000
twalk: 7260000 iterations to run, Sat Jun 3 07:54:59 2023
twalk thinning: 1 out of every 165 accepted iterations will be saved in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.out
twalk: Finished, 0.5% of moved pars per iteration (ratio 38029.303031/7260000). Output in file C:\Users\bdav_\Dropbox\Yale Paleo\AgeDepth\MALAPP_AgeDepthPollen\MALAPP_AgeDepthPollen_2JUN23\Cores/AU8/AU8_31.out,
Sat Jun 3 07:55:13 2023
bacon: 173536 WarnBeyondLimits warnings:
bacon: WARNING: calibration attempted beyond MinYr= -73.000000 or MaxYr= 1000000.000000
bacon: burn in (initial iterations which will be removed): 33000
Eso es to...eso es to...eso es to...eso es toooodo amigos!
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 3%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 12%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 18%
|
|============= | 19%
|
|============= | 20%
|
|============== | 21%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================= | 25%
|
|================= | 26%
|
|================== | 27%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|====================== | 34%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 73%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 79%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 82%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 88%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 97%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|== | 4%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|==== | 7%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|====== | 10%
|
|======= | 10%
|
|======= | 11%
|
|======== | 11%
|
|======== | 12%
|
|======== | 13%
|
|========= | 13%
|
|========= | 14%
|
|========== | 14%
|
|========== | 15%
|
|========== | 16%
|
|=========== | 16%
|
|=========== | 17%
|
|============ | 17%
|
|============ | 18%
|
|============ | 19%
|
|============= | 19%
|
|============= | 20%
|
|============== | 20%
|
|============== | 21%
|
|============== | 22%
|
|=============== | 22%
|
|=============== | 23%
|
|================ | 23%
|
|================ | 24%
|
|================ | 25%
|
|================= | 25%
|
|================= | 26%
|
|================== | 26%
|
|================== | 27%
|
|================== | 28%
|
|=================== | 28%
|
|=================== | 29%
|
|==================== | 29%
|
|==================== | 30%
|
|==================== | 31%
|
|===================== | 31%
|
|===================== | 32%
|
|====================== | 32%
|
|====================== | 33%
|
|======================= | 34%
|
|======================= | 35%
|
|======================== | 35%
|
|======================== | 36%
|
|========================= | 37%
|
|========================= | 38%
|
|========================== | 38%
|
|========================== | 39%
|
|=========================== | 40%
|
|=========================== | 41%
|
|============================ | 41%
|
|============================ | 42%
|
|============================= | 43%
|
|============================= | 44%
|
|============================== | 44%
|
|============================== | 45%
|
|=============================== | 46%
|
|=============================== | 47%
|
|================================ | 47%
|
|================================ | 48%
|
|================================= | 49%
|
|================================= | 50%
|
|================================== | 50%
|
|================================== | 51%
|
|=================================== | 52%
|
|=================================== | 53%
|
|==================================== | 53%
|
|==================================== | 54%
|
|===================================== | 55%
|
|===================================== | 56%
|
|====================================== | 56%
|
|====================================== | 57%
|
|======================================= | 58%
|
|======================================= | 59%
|
|======================================== | 59%
|
|======================================== | 60%
|
|========================================= | 61%
|
|========================================= | 62%
|
|========================================== | 62%
|
|========================================== | 63%
|
|=========================================== | 64%
|
|=========================================== | 65%
|
|============================================ | 65%
|
|============================================ | 66%
|
|============================================= | 67%
|
|============================================= | 68%
|
|============================================== | 68%
|
|============================================== | 69%
|
|=============================================== | 69%
|
|=============================================== | 70%
|
|=============================================== | 71%
|
|================================================ | 71%
|
|================================================ | 72%
|
|================================================= | 72%
|
|================================================= | 73%
|
|================================================= | 74%
|
|================================================== | 74%
|
|================================================== | 75%
|
|=================================================== | 75%
|
|=================================================== | 76%
|
|=================================================== | 77%
|
|==================================================== | 77%
|
|==================================================== | 78%
|
|===================================================== | 78%
|
|===================================================== | 79%
|
|===================================================== | 80%
|
|====================================================== | 80%
|
|====================================================== | 81%
|
|======================================================= | 81%
|
|======================================================= | 82%
|
|======================================================= | 83%
|
|======================================================== | 83%
|
|======================================================== | 84%
|
|========================================================= | 84%
|
|========================================================= | 85%
|
|========================================================= | 86%
|
|========================================================== | 86%
|
|========================================================== | 87%
|
|=========================================================== | 87%
|
|=========================================================== | 88%
|
|=========================================================== | 89%
|
|============================================================ | 89%
|
|============================================================ | 90%
|
|============================================================= | 90%
|
|============================================================= | 91%
|
|============================================================= | 92%
|
|============================================================== | 92%
|
|============================================================== | 93%
|
|=============================================================== | 93%
|
|=============================================================== | 94%
|
|=============================================================== | 95%
|
|================================================================ | 95%
|
|================================================================ | 96%
|
|================================================================= | 96%
|
|================================================================= | 97%
|
|================================================================= | 98%
|
|================================================================== | 98%
|
|================================================================== | 99%
|
|===================================================================| 99%
|
|===================================================================| 100%png
2 
#Create age-depth table
core<-rep("AU8",7)
pollenNo<-c("POL35","POL34","POL33","POL32","POL31","POL30","POL29")
depth<-c(11,45,59,74,86,110,146)
ageBP<-c()
ageMin<-c()
ageMax<-c()
for (i in depth){
ageBP<-append(ageBP,round(mean(Bacon.Age.d(i))))
ageMin<-append(ageMin,round(min(Bacon.Age.d(i))))
ageMax<-append(ageMax,round(max(Bacon.Age.d(i))))
}
AUG8data<-data.frame(core,pollenNo,depth,ageBP,ageMin,ageMax)
print(AUG8data)write.csv(AUG8data,"AUG8data.csv",row.names = FALSE)Generate pollen diagrams
Notes: 1) Terrestrial pollen data are expressed as percentages of all terrestrial species (excl. aquatics). Aquatics (Cyperaceae and Ascolepis) are plotted as percentage of all species (terrestrial + aquatics). 2) In order to view microcharcoal concentrations in the same plot as percentages, these values were transformed using log10, and then multiplied by 10 to exaggerate (e.g., conc of 10000 -> transformed log10(10000) = 4 -> exaggerated 4 * 10 = 40). The x-axis tick mark labels were then manually replaced with the appropriate scientific notation (e.g., 40 = 10000 = 10^4). 3) Age-depth modeling was not used for composite terrace sequence; ages are estimated using adjacent OSL samples (included as “TERRACEdata.csv”). 4) The convention used for the sequences is “ka”, and several of the dates are less than 1000 years and are expressed as decimals (e.g., 0.4 ka for 400 years). The rioja package does not seem to plot ages as non-integer values, so these were plotted here in terms of integer years (e.g., 400 for 400 years), but y-axis tick mark labels were edited manually to accommodate the convention. Y-axis titles (“Years ka”) were left off as well and were added manually. 5) A red line indicating disturbance in AUG7 was added manually in the published figure. A function is included here to mark the location of the disturbance, but this will look slightly different.
LS0tDQp0aXRsZTogIk1BTEFQUCBBZ2UgRGVwdGggTW9kZWxpbmcgYW5kIFBvbGxlbiBEaWFncmFtcyINCmF1dGhvcjogIkJlbmphbWluIERhdmllcywgRGF2aWQgSy4gV3JpZ2h0LCBTYXJhaCBKLiBJdm9yeSwgSmFnbyBKLiBCaXJrLCBKZW9uZy1IZW9uIENob2ksIFNhYmluZSBGaWVkbGVyLCBKYWNvYiBEYXZpcywgUG90aXBoYXIgS2FsaWJhLCBKZXNzaWNhIEMuIFRob21wc29uIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KVGhpcyBjb2RlIGlzIHVzZWQgdG8gZ2VuZXJhdGUgYWdlLWRlcHRoIG1vZGVscyBhbmQgcG9sbGVuIGRpYWdyYW1zIGZvciBhIGZvcnRoY29taW5nIHB1YmxpY2F0aW9uOiANCg0KRGF2aWQgSy4gV3JpZ2h0LCBTYXJhaCBKLiBJdm9yeSwgSmFnbyBKLiBCaXJrLCBKZW9uZy1IZW9uIENob2ksIEJlbmphbWluIERhdmllcywgU2FiaW5lIEZpZWRsZXIsIEphY29iIERhdmlzLCBQb3RpcGhhciBLYWxpYmEsIEplc3NpY2EgQy4gVGhvbXBzb24uIlBhbGFlb2Vudmlyb25tZW50YWwgZGF0YSB0cmFjayBMYXRlIFF1YXRlcm5hcnkgYW50aHJvcG9nZW5pYyBpbXBhY3RzIG9uIHZlZ2V0YXRpb24gYW5kIGxhbmRzY2FwZXMgaW4gTXppbWJhLCBub3J0aGVybiBNYWxhd2kiIFN1Ym1pdHRlZCB0byBGcm9udGllcnMgaW4gRW52aXJvbm1lbnRhbCBBcmNoYWVvbG9neSwgTWF5IDIwMjMuDQoNClRoaXMgaXMgYW4gW1IgTWFya2Rvd25dKGh0dHA6Ly9ybWFya2Rvd24ucnN0dWRpby5jb20pIE5vdGVib29rLiBXaGVuIHlvdSBleGVjdXRlIGNvZGUgd2l0aGluIHRoZSBub3RlYm9vaywgdGhlIHJlc3VsdHMgYXBwZWFyIGJlbmVhdGggdGhlIGNvZGUuIA0KDQojIyBMb2FkIHBhY2thZ2VzDQpgYGB7ciB3YXJuaW5nPUZBTFNFfQ0Kcm0obGlzdCA9IGxzKCkpDQpsaWJyYXJ5KHJiYWNvbikNCmxpYnJhcnkocmlvamEpDQpgYGANCg0KIyMjIEFVMzogR2VuZXJhdGUgYWdlIGRlcHRoIG1vZGVsIGFuZCBhc3NpZ24gcG9sbGVuIGRlcHRocw0KYGBge3IgbWVzc2FnZT1GQUxTRX0NCkJhY29uKGNvcmU9IkFVMyIsZC5taW49MCxkLm1heD02NyxhY2NlcHQuc3VnZ2VzdGlvbnMgPSBUUlVFLHZlcmJvc2U9RkFMU0UpDQpjb3JlPC1yZXAoIkFVMyIsOSkNCnBvbGxlbk5vPC1jKCJQT0wxMSIsIlBPTDEwIiwiUE9MOSIsIlBPTDgiLCJQT0wxMiIsIlBPTDEzIiwiUE9MMTYiLCJQT0wxNCIsIlBPTDE1IikNCmRlcHRoPC1jKDMsMTQsMTksMjgsMzQsNDAsNDYsNTUsNjcpDQphZ2VCUDwtYygpDQphZ2VNaW48LWMoKQ0KYWdlTWF4PC1jKCkNCmZvciAoaSBpbiBkZXB0aCl7DQogIGFnZUJQPC1hcHBlbmQoYWdlQlAscm91bmQobWVhbihCYWNvbi5BZ2UuZChpKSkpKQ0KICBhZ2VNaW48LWFwcGVuZChhZ2VNaW4scm91bmQobWluKEJhY29uLkFnZS5kKGkpKSkpDQogIGFnZU1heDwtYXBwZW5kKGFnZU1heCxyb3VuZChtYXgoQmFjb24uQWdlLmQoaSkpKSkNCn0NCkFVRzNkYXRhPC1kYXRhLmZyYW1lKGNvcmUscG9sbGVuTm8sZGVwdGgsYWdlQlAsYWdlTWluLGFnZU1heCkNCnByaW50KEFVRzNkYXRhKQ0Kd3JpdGUuY3N2KEFVRzNkYXRhLCJBVUczZGF0YS5jc3YiLHJvdy5uYW1lcyA9IEZBTFNFKQ0KYGBgDQoNCiMjIyBBVTc6IEdlbmVyYXRlIGFnZSBkZXB0aCBtb2RlbCBhbmQgYXNzaWduIHBvbGxlbiBkZXB0aHMNCmBgYHtyIG1lc3NhZ2U9RkFMU0V9DQpCYWNvbihjb3JlPSJBVTciLGQubWluPTAsZC5tYXg9NzcsYWNjZXB0LnN1Z2dlc3Rpb25zID0gVFJVRSx2ZXJib3NlPUZBTFNFKQ0KDQpjb3JlPC1yZXAoIkFVNyIsOSkNCnBvbGxlbk5vPC1jKCJQT0wzNiIsIlBPTDM3IiwiUE9MMzgiLCJQT0wzOSIsIlBPTDQwIiwiUE9MNDEiLCJQT0w0MiIsIlBPTDQzIiwiUE9MNDQiKQ0KZGVwdGg8LWMoNiwxNSwyNiwzMiw0Miw1MSw1NCw2Niw3NykNCmFnZUJQPC1jKCkNCmFnZU1pbjwtYygpDQphZ2VNYXg8LWMoKQ0KZm9yIChpIGluIGRlcHRoKXsNCiAgYWdlQlA8LWFwcGVuZChhZ2VCUCxyb3VuZChtZWFuKEJhY29uLkFnZS5kKGkpKSkpDQogIGFnZU1pbjwtYXBwZW5kKGFnZU1pbixyb3VuZChtaW4oQmFjb24uQWdlLmQoaSkpKSkNCiAgYWdlTWF4PC1hcHBlbmQoYWdlTWF4LHJvdW5kKG1heChCYWNvbi5BZ2UuZChpKSkpKQ0KfQ0KQVVHN2RhdGE8LWRhdGEuZnJhbWUoY29yZSxwb2xsZW5ObyxkZXB0aCxhZ2VCUCxhZ2VNaW4sYWdlTWF4KQ0KcHJpbnQoQVVHN2RhdGEpDQp3cml0ZS5jc3YoQVVHN2RhdGEsIkFVRzdkYXRhLmNzdiIscm93Lm5hbWVzID0gRkFMU0UpDQpgYGANCg0KIyMgQVU4OiBHZW5lcmF0ZSBhZ2UgZGVwdGggbW9kZWwgYW5kIGFzc2lnbiBwb2xsZW4gZGVwdGhzDQpgYGB7ciBtZXNzYWdlPUZBTFNFfQ0KI1J1biBhZ2UtZGVwdGggbW9kZWwgYmV0d2VlbiAwIGFuZCBtYXggZGVwdGgNCkJhY29uKGNvcmU9IkFVOCIsZC5taW49MCxkLm1heD0xNDYsZGVwdGgudW5pdCA9ICJjbSIsYWNjZXB0LnN1Z2dlc3Rpb25zID0gVFJVRSx2ZXJib3NlPUZBTFNFKQ0KI0NyZWF0ZSBhZ2UtZGVwdGggdGFibGUgDQpjb3JlPC1yZXAoIkFVOCIsNykNCnBvbGxlbk5vPC1jKCJQT0wzNSIsIlBPTDM0IiwiUE9MMzMiLCJQT0wzMiIsIlBPTDMxIiwiUE9MMzAiLCJQT0wyOSIpDQpkZXB0aDwtYygxMSw0NSw1OSw3NCw4NiwxMTAsMTQ2KQ0KYWdlQlA8LWMoKQ0KYWdlTWluPC1jKCkNCmFnZU1heDwtYygpDQpmb3IgKGkgaW4gZGVwdGgpew0KICBhZ2VCUDwtYXBwZW5kKGFnZUJQLHJvdW5kKG1lYW4oQmFjb24uQWdlLmQoaSkpKSkNCiAgYWdlTWluPC1hcHBlbmQoYWdlTWluLHJvdW5kKG1pbihCYWNvbi5BZ2UuZChpKSkpKQ0KICBhZ2VNYXg8LWFwcGVuZChhZ2VNYXgscm91bmQobWF4KEJhY29uLkFnZS5kKGkpKSkpDQp9DQpBVUc4ZGF0YTwtZGF0YS5mcmFtZShjb3JlLHBvbGxlbk5vLGRlcHRoLGFnZUJQLGFnZU1pbixhZ2VNYXgpDQpwcmludChBVUc4ZGF0YSkNCndyaXRlLmNzdihBVUc4ZGF0YSwiQVVHOGRhdGEuY3N2Iixyb3cubmFtZXMgPSBGQUxTRSkNCmBgYA0KDQoNCiMjIyBHZW5lcmF0ZSBwb2xsZW4gZGlhZ3JhbXMNCk5vdGVzOiANCjEpIFRlcnJlc3RyaWFsIHBvbGxlbiBkYXRhIGFyZSBleHByZXNzZWQgYXMgcGVyY2VudGFnZXMgb2YgYWxsIHRlcnJlc3RyaWFsIHNwZWNpZXMgKGV4Y2wuIGFxdWF0aWNzKS4gQXF1YXRpY3MgKEN5cGVyYWNlYWUgYW5kIEFzY29sZXBpcykgYXJlIHBsb3R0ZWQgYXMgcGVyY2VudGFnZSBvZiBhbGwgc3BlY2llcyAodGVycmVzdHJpYWwgKyBhcXVhdGljcykuDQoyKSBJbiBvcmRlciB0byB2aWV3IG1pY3JvY2hhcmNvYWwgY29uY2VudHJhdGlvbnMgaW4gdGhlIHNhbWUgcGxvdCBhcyBwZXJjZW50YWdlcywgdGhlc2UgdmFsdWVzIHdlcmUgdHJhbnNmb3JtZWQgdXNpbmcgbG9nMTAsIGFuZCB0aGVuIG11bHRpcGxpZWQgYnkgMTAgdG8gZXhhZ2dlcmF0ZSAoZS5nLiwgY29uYyBvZiAxMDAwMCAtPiB0cmFuc2Zvcm1lZCBsb2cxMCgxMDAwMCkgPSA0IC0+IGV4YWdnZXJhdGVkIDQgKiAxMCA9IDQwKS4gVGhlIHgtYXhpcyB0aWNrIG1hcmsgbGFiZWxzIHdlcmUgdGhlbiBtYW51YWxseSByZXBsYWNlZCB3aXRoIHRoZSBhcHByb3ByaWF0ZSBzY2llbnRpZmljIG5vdGF0aW9uIChlLmcuLCA0MCA9IDEwMDAwID0gMTBeNCkuIA0KMykgQWdlLWRlcHRoIG1vZGVsaW5nIHdhcyBub3QgdXNlZCBmb3IgY29tcG9zaXRlIHRlcnJhY2Ugc2VxdWVuY2U7IGFnZXMgYXJlIGVzdGltYXRlZCB1c2luZyBhZGphY2VudCBPU0wgc2FtcGxlcyAoaW5jbHVkZWQgYXMgIlRFUlJBQ0VkYXRhLmNzdiIpLg0KNCkgVGhlIGNvbnZlbnRpb24gdXNlZCBmb3IgdGhlIHNlcXVlbmNlcyBpcyAia2EiLCBhbmQgc2V2ZXJhbCBvZiB0aGUgZGF0ZXMgYXJlIGxlc3MgdGhhbiAxMDAwIHllYXJzIGFuZCBhcmUgZXhwcmVzc2VkIGFzIGRlY2ltYWxzIChlLmcuLCAwLjQga2EgZm9yIDQwMCB5ZWFycykuIFRoZSByaW9qYSBwYWNrYWdlIGRvZXMgbm90IHNlZW0gdG8gcGxvdCBhZ2VzIGFzIG5vbi1pbnRlZ2VyIHZhbHVlcywgc28gdGhlc2Ugd2VyZSBwbG90dGVkIGhlcmUgaW4gdGVybXMgb2YgaW50ZWdlciB5ZWFycyAoZS5nLiwgNDAwIGZvciA0MDAgeWVhcnMpLCBidXQgeS1heGlzIHRpY2sgbWFyayBsYWJlbHMgd2VyZSBlZGl0ZWQgbWFudWFsbHkgdG8gYWNjb21tb2RhdGUgdGhlIGNvbnZlbnRpb24uIFktYXhpcyB0aXRsZXMgKCJZZWFycyBrYSIpIHdlcmUgbGVmdCBvZmYgYXMgd2VsbCBhbmQgd2VyZSBhZGRlZCBtYW51YWxseS4NCjUpIEEgcmVkIGxpbmUgaW5kaWNhdGluZyBkaXN0dXJiYW5jZSBpbiBBVUc3IHdhcyBhZGRlZCBtYW51YWxseSBpbiB0aGUgcHVibGlzaGVkIGZpZ3VyZS4gQSBmdW5jdGlvbiBpcyBpbmNsdWRlZCBoZXJlIHRvIG1hcmsgdGhlIGxvY2F0aW9uIG9mIHRoZSBkaXN0dXJiYW5jZSwgYnV0IHRoaXMgd2lsbCBsb29rIHNsaWdodGx5IGRpZmZlcmVudC4gDQpgYGB7cn0NCg0KIyNBVUczDQpuYW1lczwtYygnUG9hY2VhZScsJ0FzdGVyYWNlYWUnLGV4cHJlc3Npb24oaXRhbGljKCdQb2RvY2FycHVzJykpLGV4cHJlc3Npb24oaXRhbGljKCdQaW51cycpKSxleHByZXNzaW9uKGl0YWxpYygnQnJhY2h5c3RlZ2lhJykpLCAnQ3lwZXJhY2VhZScsZXhwcmVzc2lvbihpdGFsaWMoJ0FzY29sZXBpcycpKSwiRmVybnMiLCJUcmVlcyIsIkhlcmJzIixleHByZXNzaW9uKHBhc3RlKCJNaWNyb2NoYXJjb2FsICgjL2NtIl4zLCIpIikpKQ0KDQpjb2xvcnM8LWMoImdvbGQiLCJnb2xkIiwiZm9yZXN0Z3JlZW4iLCJmb3Jlc3RncmVlbiIsImZvcmVzdGdyZWVuIiwiZGFya2N5YW4iLCJkYXJrY3lhbiIsImRhcmtncmV5IiwiZm9yZXN0Z3JlZW4iLCJnb2xkIiwiYmxhY2siKQ0KDQoNCkFVRzNEQVRFUzwtcmVhZC5jc3YoIkFVRzNkYXRhLmNzdiIsc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQ0KYWdlQlA8LUFVRzNEQVRFUyRhZ2VCUA0KQVVHM3BvbGxlbjwtcmVhZC5jc3YoIkF1ZzNQb2xsZW5fMS5jc3YiLHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkNCmNoYXI8LUFVRzNwb2xsZW4kTWljcm9jaGFyY29hbA0KY2hhcjwtKGxvZzEwKGNoYXIpKSoxMA0KY2hhclshaXMuZmluaXRlKGNoYXIpXSA8LSAwDQpBVUczcG9sbGVuJE1pY3JvY2hhcmNvYWw8LWNoYXINCnBvbGxlbkRhdGE8LUFVRzNwb2xsZW5bLDY6bmNvbChBVUczcG9sbGVuKV0NCmNzPC1jb2xTdW1zKHBvbGxlbkRhdGEpDQpwb2xsZW5EYXRhW2lzLm5hKHBvbGxlbkRhdGEpXSA8LSAwDQpwb2xsZW5EYXRhMiA8LSBwb2xsZW5EYXRhDQpleDwtYyhGQUxTRSxUUlVFLEZBTFNFLEZBTFNFLFRSVUUsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UpDQptdWx0PC01DQoNCg0KI3BuZygiQVVHM3BvbGxlbl9GSU5BTC5wbmciLHdpZHRoPTYsaGVpZ2h0PTMsdW5pdHM9ImluIixyZXM9NjAwKQ0KcGFyKG9tYT1jKDQsNCwzLDQpKQ0Kc3RyYXQucGxvdChwb2xsZW5EYXRhMix5dmFyPXJvdW5kKGFnZUJQKSwgeS50a3M9YygwLDEwMDAsMjAwMCwzMDAwKSwgeS5yZXY9VFJVRSwgcGxvdC5saW5lPUZBTFNFLCBwbG90LnBvbHk9VFJVRSwgcGxvdC5iYXI9RkFMU0UsIGNvbC5wb2x5PWNvbG9ycywgY29sLnBvbHkubGluZT0id2hpdGUiLCBsd2QucG9seT0wLjAxLCBzY2FsZS5wZXJjZW50PVRSVUUsIHhTcGFjZT0wLjAyLHgucGMuaW5jPWMoMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTApLCB4LnBjLm9taXQwPVRSVUUsY2V4LnhsYWJlbD0wLjUsY2V4LmF4aXM9MC41LGNleC55YXhpcz0wLjUseWxhYmVsID0gIiIsY2V4LnlsYWJlbD0wLjUsc3J0LnhsYWJlbD00NSxsYXM9MyxtZ3A9YygzLDAuNiwwKSxleGFnPWV4LGV4YWcubXVsdCA9IDUsY29sLmV4YWc9ImF1dG8iLGV4YWcuYWxwaGEgPSAwLjUseC5uYW1lcz1uYW1lcykNCiNkZXYub2ZmKCkNCg0KDQojQVVHNw0KI0Z1bmN0aW9uIHRvIGdlbmVyYXRlIGxpbmUgbWFya2luZyBkaXN0dXJiYW5jZSBhdCBQT0w0MiAofjEuMTI5IGthKQ0KZGlzdHVyYmVkTWFyayA8LSBmdW5jdGlvbih4LCB5LCBpLCBubSkgew0KICAjdG1wIDwtIGRhdGEuZnJhbWUoeD15LCB5PXgpDQogICN0bXAgPC0gbmEub21pdCh0bXApDQogIGxpbmVzKGMoMCwxMDApLCBjKDExMjksMTEyOSksIGNvbD0icmVkIiwgbHdkPTEpDQogIGlmIChpID09IDEyKXsNCiAgdGV4dCgyNSwxMjAwLGxhYmVscz0iRGlzdHVyYmVkIixjb2w9InJlZCIsY2V4PTAuNSkNCiAgfQ0KfQ0KDQpuYW1lczwtYygnUG9hY2VhZScsJ0FzdGVyYWNlYWUnLGV4cHJlc3Npb24oaXRhbGljKCdDb21tZWxpbmEnKSksZXhwcmVzc2lvbihpdGFsaWMoJ1BvZG9jYXJwdXMnKSksZXhwcmVzc2lvbihpdGFsaWMoJ1BpbnVzJykpLGV4cHJlc3Npb24oaXRhbGljKCdCcmFjaHlzdGVnaWEnKSksICdDeXBlcmFjZWFlJyxleHByZXNzaW9uKGl0YWxpYygnQXNjb2xlcGlzJykpLCJGZXJucyIsIlRyZWVzIiwiSGVyYnMiLGV4cHJlc3Npb24ocGFzdGUoIk1pY3JvY2hhcmNvYWwgKCMvY20iXjMsIikiKSkpDQoNCmNvbG9yczwtYygiZ29sZCIsImdvbGQiLCJnb2xkIiwiZm9yZXN0Z3JlZW4iLCJmb3Jlc3RncmVlbiIsImZvcmVzdGdyZWVuIiwiZGFya2N5YW4iLCJkYXJrY3lhbiIsImRhcmtncmV5IiwiZm9yZXN0Z3JlZW4iLCJnb2xkIiwiYmxhY2siKQ0KDQpBVUc3REFURVM8LXJlYWQuY3N2KCJBVUc3ZGF0YS5jc3YiLHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkNCmFnZUJQPC1BVUc3REFURVMkYWdlQlANCkFVRzdwb2xsZW48LXJlYWQuY3N2KCJBdWc3UG9sbGVuXzEuY3N2IixzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpDQpjaGFyPC1BVUc3cG9sbGVuJE1pY3JvY2hhcmNvYWwNCmNoYXI8LShsb2cxMChjaGFyKSkqMTANCmNoYXJbIWlzLmZpbml0ZShjaGFyKV0gPC0gMA0KQVVHN3BvbGxlbiRNaWNyb2NoYXJjb2FsPC1jaGFyDQpwb2xsZW5EYXRhPC1BVUc3cG9sbGVuWyw2Om5jb2woQVVHN3BvbGxlbildDQpjczwtY29sU3Vtcyhwb2xsZW5EYXRhKQ0KcG9sbGVuRGF0YVtpcy5uYShwb2xsZW5EYXRhKV0gPC0gMA0KcG9sbGVuRGF0YTIgPC0gcG9sbGVuRGF0YQ0KZXg8LWMoRkFMU0UsVFJVRSxGQUxTRSxGQUxTRSxUUlVFLFRSVUUsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UsRkFMU0UpDQptdWx0PC01DQoNCg0KI3BuZygiQVU3cG9sbGVuX0ZJTkFMLnBuZyIsd2lkdGg9NixoZWlnaHQ9Myx1bml0cz0iaW4iLHJlcz02MDApDQpwYXIob21hPWMoNCw0LDMsNCkpDQpzdHJhdC5wbG90KHBvbGxlbkRhdGEyLHl2YXI9cm91bmQoYWdlQlApLCB5LnRrcz1jKDAsNTAwLDEwMDAsMTUwMCwyMDAwLDI1MDApLCB5LnJldj1UUlVFLCBwbG90LmxpbmU9RkFMU0UsIHBsb3QucG9seT1UUlVFLCBwbG90LmJhcj1GQUxTRSwgY29sLnBvbHk9Y29sb3JzLCBjb2wucG9seS5saW5lPSJ3aGl0ZSIsIGx3ZC5wb2x5PTAuMDEsIHNjYWxlLnBlcmNlbnQ9VFJVRSwgeFNwYWNlPTAuMDIseC5wYy5pbmM9YygxMCwxMCwxMCwxMCwxMCwxMCwxMCwxMCwxMCwxMCwxMCwxMCksIHgucGMub21pdDA9VFJVRSxjZXgueGxhYmVsPTAuNSxjZXguYXhpcz0wLjUsY2V4LnlheGlzPTAuNSx5bGFiZWwgPSAiIixjZXgueWxhYmVsPTAuNSxzcnQueGxhYmVsPTQ1LGxhcz0zLG1ncD1jKDMsMC42LDApLGV4YWc9ZXgsZXhhZy5tdWx0ID0gNSxjb2wuZXhhZz0iYXV0byIsZXhhZy5hbHBoYSA9IDAuNSx4Lm5hbWVzPW5hbWVzLGZ1bjI9ZGlzdHJ1cmJlZE1hcmspDQoNCiNkZXYub2ZmKCkNCg0KDQojQVVHIDgNCg0KbmFtZXM8LWMoJ1BvYWNlYWUnLCdBc3RlcmFjZWFlJyxleHByZXNzaW9uKGl0YWxpYygnUG9kb2NhcnB1cycpKSxleHByZXNzaW9uKGl0YWxpYygnUGludXMnKSksJ0N5cGVyYWNlYWUnLGV4cHJlc3Npb24oaXRhbGljKCdBc2NvbGVwaXMnKSksIkZlcm5zIiwiVHJlZXMiLCJIZXJicyIsZXhwcmVzc2lvbihwYXN0ZSgiTWljcm9jaGFyY29hbCAoIy9jbSJeMywiKSIpKSkNCg0KY29sb3JzPC1jKCJnb2xkIiwiZ29sZCIsImZvcmVzdGdyZWVuIiwiZm9yZXN0Z3JlZW4iLCJkYXJrY3lhbiIsImRhcmtjeWFuIiwiZGFya2dyZXkiLCJmb3Jlc3RncmVlbiIsImdvbGQiLCJibGFjayIpDQoNCkFVRzhEQVRFUzwtcmVhZC5jc3YoIkFVRzhkYXRhLmNzdiIsc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQ0KYWdlQlA8LUFVRzhEQVRFUyRhZ2VCUA0KQVVHOHBvbGxlbjwtcmVhZC5jc3YoIkF1ZzhQb2xsZW5fMS5jc3YiLHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkNCmNoYXI8LUFVRzhwb2xsZW4kTWljcm9jaGFyY29hbA0KY2hhcjwtKGxvZzEwKGNoYXIpKSoxMA0KY2hhclshaXMuZmluaXRlKGNoYXIpXSA8LSAwDQpBVUc4cG9sbGVuJE1pY3JvY2hhcmNvYWw8LWNoYXINCnBvbGxlbkRhdGE8LUFVRzhwb2xsZW5bLDY6bmNvbChBVUc4cG9sbGVuKV0NCmNzPC1jb2xTdW1zKHBvbGxlbkRhdGEpDQpwb2xsZW5EYXRhW2lzLm5hKHBvbGxlbkRhdGEpXSA8LSAwDQpwb2xsZW5EYXRhMiA8LSBwb2xsZW5EYXRhDQpleDwtYyhGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxUUlVFLEZBTFNFLEZBTFNFLEZBTFNFLEZBTFNFKQ0KbXVsdDwtNQ0KDQoNCiNwbmcoIkFVOHBvbGxlbl9GSU5BTC5wbmciLHdpZHRoPTYsaGVpZ2h0PTMsdW5pdHM9ImluIixyZXM9NjAwKQ0KcGFyKG9tYT1jKDQsNCwzLDQpKQ0Kc3RyYXQucGxvdChwb2xsZW5EYXRhMix5dmFyPXJvdW5kKGFnZUJQKSwgeS50a3M9YygwLDEwMDAwLDIwMDAwLDMwMDAwKSwgeS5yZXY9VFJVRSwgcGxvdC5saW5lPUZBTFNFLCBwbG90LnBvbHk9VFJVRSwgcGxvdC5iYXI9RkFMU0UsIGNvbC5wb2x5PWNvbG9ycywgY29sLnBvbHkubGluZT0id2hpdGUiLCBsd2QucG9seT0wLjAxLCBzY2FsZS5wZXJjZW50PVRSVUUsIHhTcGFjZT0wLjAyLHgucGMuaW5jPWMoMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTApLCB4LnBjLm9taXQwPVRSVUUsY2V4LnhsYWJlbD0wLjUsY2V4LmF4aXM9MC41LGNleC55YXhpcz0wLjUseWxhYmVsID0gIlllYXJzIEJQIixjZXgueWxhYmVsPTAuNSxzcnQueGxhYmVsPTQ1LGxhcz0zLG1ncD1jKDMsMC42LDApLGV4YWc9ZXgsZXhhZy5tdWx0ID0gNSxjb2wuZXhhZz0iYXV0byIsZXhhZy5hbHBoYSA9IDAuNSx4Lm5hbWVzPW5hbWVzKQ0KI2Rldi5vZmYoKQ0KDQoNCg0KI1RFUlJBQ0VTDQoNCm5hbWVzPC1jKCdQb2FjZWFlJywnQXN0ZXJhY2VhZScsJ1Bpbm9wc2lkYSB1bmRpZmYuJyxleHByZXNzaW9uKGl0YWxpYygnQnJhY2h5c3RlZ2lhJykpLCdDeXBlcmFjZWFlJywiRmVybnMiLCJUcmVlcyIsIkhlcmJzIixleHByZXNzaW9uKHBhc3RlKCJNaWNyb2NoYXJjb2FsICgjL2NtIl4zLCIpIikpKQ0KDQpwb2x5Y29sb3JzPC1jKCJnb2xkIiwiZ29sZCIsImZvcmVzdGdyZWVuIiwiZm9yZXN0Z3JlZW4iLCJkYXJrY3lhbiIsImRhcmtncmV5IiwiZm9yZXN0Z3JlZW4iLCJnb2xkIiwiYmxhY2siKQ0KDQoNClRFUlJBQ0VEQVRFUzwtcmVhZC5jc3YoIlRFUlJBQ0VkYXRhLmNzdiIsc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQ0KYWdlS0E8LVRFUlJBQ0VEQVRFUyRhZ2VLQQ0KVEVSUkFDRXBvbGxlbjwtcmVhZC5jc3YoIlRlcnJhY2VQb2xsZW5fMS5jc3YiLHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSkNCmNoYXI8LVRFUlJBQ0Vwb2xsZW4kTWljcm9jaGFyY29hbA0KY2hhcjwtKGxvZzEwKGNoYXIpKSoxMA0KY2hhclshaXMuZmluaXRlKGNoYXIpXSA8LSAwDQpURVJSQUNFcG9sbGVuJE1pY3JvY2hhcmNvYWw8LWNoYXINCnBvbGxlbkRhdGE8LVRFUlJBQ0Vwb2xsZW5bLDY6bmNvbChURVJSQUNFcG9sbGVuKV0NCmNzPC1jb2xTdW1zKHBvbGxlbkRhdGEpDQpwb2xsZW5EYXRhW2lzLm5hKHBvbGxlbkRhdGEpXSA8LSAwDQpwb2xsZW5EYXRhMiA8LSBwb2xsZW5EYXRhI1ssIHdoaWNoKGNzID4gNSldDQpleDwtYyhGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSxGQUxTRSkNCm11bHQ8LTUNCg0KDQojcG5nKCJURVJSQUNFcG9sbGVuX0ZJTkFMMi5wbmciLHdpZHRoPTYsaGVpZ2h0PTMsdW5pdHM9ImluIixyZXM9NjAwKQ0KcGFyKG9tYT1jKDQsNCwzLDQpKQ0Kc3RyYXQucGxvdChwb2xsZW5EYXRhMix5dmFyPXJvdW5kKGFnZUtBKSwgeS50a3M9YygwLDIwMDAsNDAwMCw2MDAwLDgwMDApLCB5LnJldj1UUlVFLCBwbG90LmxpbmU9RkFMU0UsIHBsb3QucG9seT1UUlVFLCBwbG90LmJhcj1GQUxTRSwgY29sLnBvbHk9cG9seWNvbG9ycywgY29sLnBvbHkubGluZT0id2hpdGUiLCBsd2QucG9seT0wLjAxLCBzY2FsZS5wZXJjZW50PVRSVUUsIHhTcGFjZT0wLjAyLHgucGMuaW5jPWMoMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTAsMTApLCB4LnBjLm9taXQwPVRSVUUsY2V4LnhsYWJlbD0wLjUsY2V4LmF4aXM9MC41LGNleC55YXhpcz0wLjUseWxhYmVsID0gIiIsY2V4LnlsYWJlbD0wLjUsc3J0LnhsYWJlbD00NSxsYXM9MyxtZ3A9YygzLDAuNiwwKSxleGFnPWV4LGV4YWcubXVsdCA9IDUsY29sLmV4YWc9ImF1dG8iLGV4YWcuYWxwaGEgPSAwLjUseC5uYW1lcz1uYW1lcykNCiNkZXYub2ZmKCkNCg0KDQpgYGA=