rm(list=ls())
#--------------------------------------------------------------------------
#   Berechnung von Sedimentationszeiten
#--------------------------------------------------------------------------
#
# Funktion zur Berechnung der Fallzeiten:
sed_time <- function(z,d) 
{
  sed_time <- 10*(z*18*visc)/(g_acc*(rho_s - rho_w)*d^2)
  # *10 wegen unterschiedl. Einheiten
}
#--------------------------------------------------------------------------
#   Grunddaten
#--------------------------------------------------------------------------
ndata_z <- 10
rho_s <- 2.65                                # g/cm^3
rho_w <- 1.00                                # g/cm^3
visc  <- 0.001002                            # kg/(m sec)
g_acc <- 9.81                                # m/sec^2
#
# Datenvektor mit Korndurchmessern
d <- c(0.2,0.063,0.02,0.0063,0.002)          # mm
#
# Datenvektor mit Tiefen:
z_min <- 5; z_max <- 50                      # cm
z <- seq(z_min,z_max,by=(z_max - z_min)/(ndata_z-1))
#
# Erzeugen der Datenmatrix mit den Funktionen "outer" und "sed_time"
sed_t <- outer(z,d,FUN="sed_time")
# 
A <- data.frame(round(sed_t[,1:5],1),round(z,1)) 
names(A) <- c(paste(d, c(" mm"), sep=""),"z")
cat("---------------------------------------------------------\n")
cat(" Fallzeiten (sec) für verschiedene Tiefen und Korngrößen:\n")
cat("---------------------------------------------------------\n")
print(A)
cat("---------------------------------------------------------\n")
#write.table(A,"Fallzeiten.txt",sep="\t")
#
#   --- Schluss ---
#