rm(list=ls())
library(DAAG)
#-------------------------------------------
#  Graphische Darstellung zu Profildaten
#-------------------------------------------
#
#  --- Dateneingabe ---
#
daten <- read.table("profile_erg.txt");daten
summary(daten)  # statistischer Datenüberblick
names(daten)    # zeigt die Spaltennamen
# Bem.: mit names(daten) <- c("v1","v2",usw) 
# können die Spaltentitel umbenannt werden.
attach(daten)   # macht die Spalten direkt über die Namen zugänglich
par(mfrow = c(2,2)) # Untereilt das Graphikfenster in 2 Reihen und zwei Spalten
#
#  --- Aufteilung des Porenraums ---
#
# Leere Grafik mit Achsen und Beschriftungen:
plot(Wco_mu,z,type="n",xlim=c(0,max(Wco_mu,Pco_mu+Pco_sd)),ylim=c(100,0),
	main="pore space",xlab="content (% vol)",ylab="depth (cm)",
	cex.main=2.0,cex.lab=1.5)
pause()
#
# Porenraum:
lines(Pco_mu,z,lwd=2);
arrows(Pco_mu-Pco_sd,z,Pco_mu+Pco_sd,z,code=3,angle=90,lwd=2,length=0.05)
points(Pco_mu,z,pch=15,cex=1.5,col="red");
#
# Wassergehalte
lines(Wco_mu,z,lwd=2);
arrows(Wco_mu-Wco_sd,z,Wco_mu+Wco_sd,z,code=3,angle=90,lwd=2,length=0.05)
points(Wco_mu,z,pch=15,cex=1.5,col="blue")
#
# Luftgehalte
lines(Aco_mu,z,lwd=2);
arrows(Aco_mu-Aco_sd,z,Aco_mu+Aco_sd,z,code=3,angle=90,lwd=2,length=0.05)
points(Aco_mu,z,pch=15,cex=1.5,col="green")
#
# Legende (Position fix bei (0,09)
legend(0,0,c("total","water", "air"),pch=c(15,15,15),col=c("red","blue","green"))
#
#  --- Aufteilung des Bodenvolumens ---
#
plot(Sco_mu,z,type="n",xlim=c(0,105),ylim=c(100,0),
	main="soil volume",xlab="content (% vol)",ylab="depth (cm)",
	cex.main=2.0,cex.lab=1.5)
#
# Bodenvolumen gesamt (= SV+WV+LV = 100 %):
# Die gefärbten Flächen werden mit dem polygon-Befehl erzeugt;
# anschließend werden mit points noch Datenpunkte angezeit
BV <- Sco_mu + Wco_mu + Aco_mu 
polygon(c(0,BV[1]-0.5*(BV[2]-BV[1]),BV,BV[10]+0.5*(BV[10]-BV[9]),0),
        c(0,0,z,100,100),col="lightgrey",lwd=1)
points(BV,z,pch=15,cex=0.8,col="black")
# 
# Substanz- + Wasservolumen (= SV+WV in %) :
WV <- Sco_mu + Wco_mu 
polygon(c(0,WV[1]-0.5*(WV[2]-WV[1]),WV,WV[10]+0.5*(WV[10]-WV[9]),0),
        c(0,0,z,100,100),col="lightblue",lwd=1)
points(WV,z,pch=15,cex=0.8,col="blue")
#
# Substanzvolumen (%):
SV <- Sco_mu 
polygon(c(0,SV[1]-0.5*(SV[2]-SV[1]),SV,SV[10]+0.5*(SV[10]-SV[9]),0),
        c(0,0,z,100,100),col="darkgrey",lwd=1)
points(SV,z,pch=15,cex=0.8,col="yellow")
#
legend(0,0,c("air","water","solids"),pch=c(15,15,15),col=c("black","blue","yellow"))
#
#  --- Lagerungsdichten ---
#
plot(Mbd_mu,z,type="n",xlim=c(0,max(Mbd_mu+Mbd_sd)),ylim=c(100,0),
	main="bulk density",xlab="density(g/cm^3)",ylab="depth (cm)",
	cex.main=2.0,cex.lab=1.5)
#
# Trockenraumdichte:
lines(Dbd_mu,z,lwd=2);
arrows(Dbd_mu-Dbd_sd,z,Dbd_mu+Dbd_sd,z,code=3,angle=90,lwd=2,length=0.05)
points(Dbd_mu,z,pch=15,cex=1.5,col="red");
# 
# Feuchte Raumdichte:
lines(Mbd_mu,z,lwd=2);
arrows(Mbd_mu-Mbd_sd,z,Mbd_mu+Mbd_sd,z,code=3,angle=90,lwd=2,length=0.05)
points(Mbd_mu,z,pch=15,cex=1.5,col="blue")
#
# Legende:
legend(0,0,c("dry bulk density","moist bulk density"),pch=c(15,15),col=c("red","blue"))
#
#  --- Auflast ---
#
plot(W_t,z,type="n",xlim=c(0,max(W_t)),ylim=c(100,0),
	main="normal stress",xlab="density(g/cm^3)",ylab="depth (cm)",
	cex.main=2.0,cex.lab=1.5)
# 
# Auflast durch trockenen Boden (W_d, in N/m^2):
lines(W_d,z,lwd=2);
points(W_d,z,pch=15,cex=1.5,col="red");
#
# Auflast durch feuchten Boden (W_t, in N/m^2):
lines(W_t,z,lwd=2);
points(W_t,z,pch=15,cex=1.5,col="blue")
#
legend(0,70,c("dry soil","moist soil"),pch=c(15,15),col=c("red","blue"))
#
detach(daten) # verdeckt die Spaltennamen 
#
# --- Ende ---