x<-c(1,1,1.15,1.15)
lito.1<-c(60,0,0,60)
lito.2<-c(60,63,63,60)
lito.3<-c(63,83,83,63)
lito.4<-c(83,200,200,83)
lito.5<-c(200,290,290,200)
lito.6<-c(290,310,310,290)
lito.7<-c(310,360,360,310)
#zonas
col<-c(1.25,1.25,1.40,1.40)
zona.1<-c(52.5,0,0,52.5)
zona.2<-c(52.5,58,58,52.5)
zona.3<-c(58,74,74,58)
zona.4<-c(58,82,82,58)
zona.5<-c(82,200,200,82)
zona.6<-c(200,290,290,200)
zona.7<-c(290,310,310,290)
zona.8<-c(310,360,360,310)
#########
par(mar = c(4, 3.5, 3.5, 0.5)) #c(bottom, left, top, right)
par(tcl=-0.3) # modifica o argumento sempre q for usado
par(mgp = c(1.6, -7, -8))
par(xpd=NA)
plot(1,1,xlim=c(0,3),ylim=c(360,0),type="l",col="black",lwd=20,
xlab="",ylab="",axe=F)
polygon(x,lito.1,density=1000,col="black")
polygon(x,lito.2,density=0)
polygon(x,lito.3,density=0)
polygon(x,lito.3,density=30,angle=180,lty=3)
polygon(x,lito.4,density=1000,col="black")
polygon(x,lito.5,density=20,angle=45)
polygon(x,lito.6,density=45,angle=90)
polygon(x,lito.7,density=20,angle=45)
polygon(col,zona.1,lty=3)
polygon(col,zona.2,lty=3)
polygon(col,zona.3,lty=3)
polygon(col,zona.4,lty=3)
polygon(col,zona.5,lty=3)
polygon(col,zona.6,lty=3)
polygon(col,zona.7,lty=3)
polygon(col,zona.8,lty=3)
#### Add idades
text(.67,14,"2.910±60",cex=.65)
la<-c(1,.9,.88)
lb<-c(52.5,52.5,48)
lines(la,lb,type="l")
text(.67,46,"6.590±150",cex=.65)
pa<-c(1,.9,.88)
pb<-c(56,56,56)
lines(pa,pb,type="l")
text(.67,55,"8.250±100",cex=.65)
da<-c(1,.9,.88)
db<-c(67.5,67.5,69)
lines(da,db,type="l")
text(.67,68.5,"18.740±140",cex=.65)
ua<-c(1,.9,.88)
ub<-c(74.5,74.5,75)
lines(ua,ub,type="l")
text(.67,75,"21.260±840",cex=.65)
va<-c(1,.9,.88)
vb<-c(82,82,82)
lines(va,vb,type="l")
text(.67,82,"27.320±500",cex=.65)
ha<-c(1,.9,.88)
hb<-c(84.5,84.5,89)
lines(ha,hb,type="l")
text(.67,89,"35.220±430",cex=.65)
sa<-c(1,.9,.88)
sb<-c(111,111,111)
lines(sa,sb,type="l")
text(.67,111,"41.690±430",cex=.65)
fa<-c(1,.9,.88)
fb<-c(117.5,117.5,124)
lines(fa,fb,type="l")
text(.65,124,"43.200±1320",cex=.65)
text(.67,157.5,"45.650±2300",cex=.65)
#Nomes das zonas
text(1.325,245,"1-A",cex=.65)
text(1.325,130,"1-B",cex=.65)
text(1.325,78,"2-A",cex=.65)
text(1.325,65,"2-B",cex=.65)
text(1.325,55,"3-A",cex=.65)
text(1.325,30,"3-B",cex=.65)
####Keys
K<-c(1,1,1.15,1.15)
Ka<-c(1.6,1.6,1.75,1.75)
Kb<-c(2.2,2.2,2.35,2.35)
Kc<-c(2.8,2.8,2.95,2.95)
K.1<-c(386,370,370,386)
polygon(K,K.1,density=1000,col="black")
polygon(Ka,K.1,density=30,angle=180,lty=3)
polygon(Ka,K.1)
polygon(Kb,K.1,density=20,angle=45)
polygon(Kc,K.1,density=50,angle=90)
text(1.07,395,"Black gyttja",cex=.65)
text(1.66,395,"Nodular yellowish clay",cex=.65)
text(2.32,395,"Brown silty clay",cex=.65)
text(2.9,395,"Dark brown silty clay",cex=.65)
##### Títulos
par(srt="45")
text(1.14,-22,"Lithology",cex=.65)
text(.8,-28,"ages (cal yr BP)",cex=.65)
text(1.67,-22,"Depth (cm)",cex=.65)
text(1.35,-19.5,"Zones",cex=.65)
