#------------------------------------------------------------------
#	 Aufteilung des Bodens in Flüssig-, Gas- und Festphase 
#------------------------------------------------------------------
library(DAAG)
rm(list=ls()) # löscht alle Objekte im workspace
#
# Eingabedaten
#
rho_s <- 2.64 # Substanzdichte (Dezimalzahlen nie mit ",")
rho_h2o <- 1  # Dichte von Wasser
SZ_feu <- 215 # Masse des feuchten Bodens
SZ_tro <- 185 # Masse des trockenen Bodens
SZ_leer <- 75 # Masse des leeren Stechzylinders
V_sz <- 100   # Volumen des Stechzylinders
Tiefe <- 30   # Bodentiefe in cm
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*100)/V_sz
# Luftgehalt:
lg <- por - theta_vol
# Porenziffer:
pz <- por/(1-por)
# Wasserspeicherung in 0-30 cm Tiefe:
WSV <- Tiefe*10*theta_vol
#---------------------------
# Ausgabe der Ergebnisse
#---------------------------
cat("Ergebnisse\n")
cat("==========\n")
cat("Porosität  :",signif(por*100,dig=3)," % \n")
cat("Vol. Wassergehalt:",signif(theta_vol*100,dig=3)," % \n")
cat("Vol. Luftgehalt  :",signif(lg*100,3)," % \n")
cat("Vol. Festsubstanz:",signif(fs_vol*100,3)," % \n")
cat("Lagerungsdichte  :",signif(rho_b,3)," g cm^-3\n")
cat("Porenziffer  :",signif(pz,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---