#------------------------------------------------------------------
#	 Aufteilung des Bodens in Flüssig-, Gas- und Festphase 
#------------------------------------------------------------------
library(DAAG)
rm(list=ls()) # löscht alle Objekte im workspace
#
# Eingabedaten
#
rho_s <- 2.65 # Substanzdichte (Dezimalzahlen nie mit ",")
rho_h2o <- 1  # Dichte von Wasser
SZ_feu <- 250 # Masse des feuchten Bodens
SZ_tro <- 200 # Masse des trockenen Bodens
SZ_leer <- 95 # Masse des leeren Stechzylinders
V_sz <- 100   # Volumen des Stechzylinders
pause()
#---------------------------
# Berechnung der drei Phasen
#---------------------------
# Vol. Wassergehalt
#
theta_vol <- (SZ_feu - SZ_tro)/rho_h2o/V_sz
#
# Trockenraumdichte
#
rho_b <- (SZ_tro - SZ_leer)/V_sz
#
# Porosität
#
por <- 1 - rho_b/rho_s
#
# Festphase
#
FS_vol <- (V_sz - por)/V_sz
#
# Luftgehalt
#
Lg <- por - theta_vol
#
#---------------------------
# Ausgabe der Ergebnisse
#---------------------------
cat("Ergebnisse\n")
cat("==========\n")
cat("Vol. Wassergehalt:",theta_vol," %\n")
cat("Vol. Luftgehalt  :",Lg," %\n")
cat("Vol. Festsubstanz:",FS_vol," %\n")
cat("Lagerungsdichte  :",rho_b," g cm^-3\n")
pause()
#---------------------------
# Pie-chart
#---------------------------
par(lwd=3, col="black", cex = 1.5)
lbs <- c("WV","LV","SV")
dat <- c(theta_vol,Lg,FS_vol)
cls <- c("blue1","lightgrey","sienna")
pie(dat,lbs,col=cls,main="Wasser-, Luft- und Festphase")
#
#---Ende---