Calculate drift percentages and PECsw using the Rautmann formula

3 files changed, 53 insertions, 3 deletions

diff --git a/man/PEC_sw_drift.Rd b/man/PEC_sw_drift.Rd
index 9fd733e..47a1e41 100644
--- a/man/PEC_sw_drift.Rd
+++ b/man/PEC_sw_drift.Rd
@@ -10,7 +10,10 @@ PEC_sw_drift(
   water_depth = 30,
   drift_percentages = NULL,
   drift_data = c("JKI", "RF"),
-  crop = "Ackerbau",
+  crop_group_JKI = c("Ackerbau", "Obstbau frueh", "Obstbau spaet", "Weinbau frueh",
+    "Weinbau spaet", "Hopfenbau", "Flaechenkulturen > 900 l/ha", "Gleisanlagen"),
+  crop_group_focus = c("arable", "hops", "vines, late", "vines, early", "fruit, late",
+    "fruit, early", "aerial"),
   distances = c(1, 5, 10, 20),
   rate_units = "g/ha",
   PEC_units = "µg/L"
@@ -30,7 +33,10 @@ PEC_sw_drift(
 included in the package is used. If 'RF', the Rautmann formula is used, if
 implemented for the crop type and number of applications}
 
-\item{crop}{Crop name (use German names for JKI data), defaults to "Ackerbau"}
+\item{crop_group_JKI}{When using the 'JKI' drift data, one of the German names
+as used in \link{drift_parameters_focus}.}
+
+\item{crop_group_focus}{One of the crop groups as used in \link{drift_parameters_focus}}
 
 \item{distances}{The distances in m for which to get PEC values}
 
@@ -51,9 +57,13 @@ PEC_sw_drift(100)
 # Alternatively, we can use the formula for a single application to "Ackerbau" from the paper
 PEC_sw_drift(100, drift_data = "RF")
 
-# This makes it possible to also use different substances
+# This makes it possible to also use different distances
 PEC_sw_drift(100, distances = c(1, 3, 5, 6, 10, 20, 50, 100), drift_data = "RF")
 
+# or consider aerial application
+PEC_sw_drift(100, distances = c(1, 3, 5, 6, 10, 20, 50, 100), drift_data = "RF", 
+  crop_group_focus = "aerial")
+
 # Using custom drift percentages is also supported
 PEC_sw_drift(100, drift_percentages = c(2.77, 0.95, 0.57, 0.48, 0.29, 0.15, 0.06, 0.03))
 }
diff --git a/man/drift_parameters_focus.Rd b/man/drift_parameters_focus.Rd
index 2ba6a7a..61447ad 100644
--- a/man/drift_parameters_focus.Rd
+++ b/man/drift_parameters_focus.Rd
@@ -16,8 +16,15 @@ using the R code given in the file \code{data_generation/drift_parameters_focus.
 installed with this package. The appendix itself is not included in the package,
 as its licence is not clear.
 }
+\details{
+For the hinge distance, \code{Inf} was substituted for the cases where no hinge
+distance is given in the data, in this way parameters C and D are never
+used for any distance if A and B are used for the case that the distance
+is smaller than the hinge distance.
+}
 \examples{
 drift_parameters_focus
+unique(drift_parameters_focus$crop_group)
 }
 \references{
 FOCUS (2014) Generic guidance for Surface Water Scenarios (version 1.4).
diff --git a/man/drift_percentages_rautmann_formula.Rd b/man/drift_percentages_rautmann_formula.Rd
new file mode 100644
index 0000000..6823a49
--- /dev/null
+++ b/man/drift_percentages_rautmann_formula.Rd
@@ -0,0 +1,33 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PEC_sw_drift.R
+\name{drift_percentages_rautmann_formula}
+\alias{drift_percentages_rautmann_formula}
+\title{Calculate the drift percentages according to the Rautmann formula}
+\usage{
+drift_percentages_rautmann_formula(
+  distances,
+  applications = 1,
+  crop_group_focus = c("arable", "hops", "vines, late", "vines, early", "fruit, late",
+    "fruit, early", "aerial")
+)
+}
+\arguments{
+\item{distances}{The distances in m for which to get PEC values}
+
+\item{applications}{Number of applications for selection of drift percentile}
+
+\item{crop_group_focus}{One of the crop groups as used in \link{drift_parameters_focus}}
+}
+\description{
+Calculate the drift percentages according to the Rautmann formula
+}
+\examples{
+drift_data_JKI[[1]][as.character(c(1, 3, 5)), "Ackerbau"]
+drift_percentages_rautmann_formula(c(1, 3, 5))
+
+# One application to early or late fruit crops
+drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau frueh"]
+drift_percentages_rautmann_formula(c(3, 5, 20, 50), crop_group = "fruit, early")
+drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau spaet"]
+drift_percentages_rautmann_formula(c(3, 5, 20, 50), crop_group = "fruit, late")
+}