Extend drift percentage calcs with FOCUS method

5 files changed, 67 insertions, 41 deletions

diff --git a/man/PEC_sw_drift.Rd b/man/PEC_sw_drift.Rd
index fd88fc8..5d56c8a 100644
--- a/man/PEC_sw_drift.Rd
+++ b/man/PEC_sw_drift.Rd
@@ -68,7 +68,7 @@ PEC_sw_drift(100, distances = c(1, 3, 5, 6, 10, 20, 50, 100), drift_data = "RF",
 PEC_sw_drift(100, drift_percentages = c(2.77, 0.95, 0.57, 0.48, 0.29, 0.15, 0.06, 0.03))
 }
 \seealso{
-\link{drift_parameters_focus}, \link{drift_percentages_rautmann_formula}
+\link{drift_parameters_focus}, \link{drift_percentages_rautmann}
 }
 \author{
 Johannes Ranke
diff --git a/man/SSLRC_mobility_classification.Rd b/man/SSLRC_mobility_classification.Rd
index 43de457..daeebc2 100644
--- a/man/SSLRC_mobility_classification.Rd
+++ b/man/SSLRC_mobility_classification.Rd
@@ -15,7 +15,7 @@ compound transported per 10 mm drain water
 }
 \description{
 This implements the method specified in the UK data requirements handbook and was
-checked against the spreadsheet published on the CRC website
+checked against the spreadsheet published on the CRD website
 }
 \examples{
 SSLRC_mobility_classification(100)
diff --git a/man/drift_parameters_focus.Rd b/man/drift_parameters_focus.Rd
index 0320f7f..fe69784 100644
--- a/man/drift_parameters_focus.Rd
+++ b/man/drift_parameters_focus.Rd
@@ -41,6 +41,6 @@ the authorization procedure for plant protection products Mitt. Biol.
 Bundesanst. Land- Forstwirtsch. 383, 133-141
 }
 \seealso{
-\link{drift_percentages_rautmann_formula}, \link{PEC_sw_drift}
+\link{drift_percentages_rautmann}, \link{PEC_sw_drift}
 }
 \keyword{datasets}
diff --git a/man/drift_percentages_rautmann.Rd b/man/drift_percentages_rautmann.Rd
new file mode 100644
index 0000000..c3686ff
--- /dev/null
+++ b/man/drift_percentages_rautmann.Rd
@@ -0,0 +1,64 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PEC_sw_drift.R
+\name{drift_percentages_rautmann}
+\alias{drift_percentages_rautmann}
+\title{Calculate drift percentages based on Rautmann data}
+\usage{
+drift_percentages_rautmann(
+  distances,
+  applications = 1,
+  crop_group_focus = c("arable", "hops", "vines, late", "vines, early", "fruit, late",
+    "fruit, early", "aerial"),
+  formula = c("Rautmann", "FOCUS"),
+  widths = 1
+)
+}
+\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}}
+
+\item{formula}{By default, the original Rautmann formula is used. If you
+specify "FOCUS", mean drift input over the width of the water body is
+calculated as described in Chapter 5.4.5 of the FOCUS surface water guidance}
+
+\item{widths}{The widths of the water bodies (only used in the FOCUS formula)}
+}
+\description{
+Calculate drift percentages based on Rautmann data
+}
+\examples{
+# Compare JKI data with Rautmann formula
+# One application on field crops, for 1 m, 3 m and 5 m distance
+drift_data_JKI[[1]][as.character(c(1, 3, 5)), "Ackerbau"]
+drift_percentages_rautmann(c(1, 3, 5))
+drift_percentages_rautmann(c(1, 3, 5), formula = "FOCUS")
+
+# One application to early or late fruit crops
+drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau frueh"]
+drift_percentages_rautmann(c(3, 5, 20, 50), crop_group = "fruit, early")
+drift_percentages_rautmann(c(3, 5, 20, 50), crop_group = "fruit, early", formula = "FOCUS")
+drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau spaet"]
+drift_percentages_rautmann(c(3, 5, 20, 50), crop_group = "fruit, late")
+drift_percentages_rautmann(c(3, 5, 20, 50), crop_group = "fruit, late", formula = "FOCUS")
+
+# We get a continuum if the waterbody covers the hinge distance (11.4 m for 1 early app to fruit)
+drift_percentages_rautmann(seq(10, 13, by = 0.2), crop_group = "fruit, early", formula = "FOCUS")
+x <- seq(1, 30, by = 0.1)
+d <- drift_percentages_rautmann(x, crop_group = "fruit, early", formula = "FOCUS")
+plot(x, d, type = "l",
+  xlab = "Distance of near edge [m]", 
+  ylab = "Mean drift percentage over waterbody width",
+  main = "One application to fruit, early")
+abline(v = 11.4, lty = 2)
+}
+\references{
+FOCUS (2014) Generic guidance for Surface Water Scenarios (version 1.4).
+FOrum for the Co-ordination of pesticde fate models and their USe.
+\url{http://esdac.jrc.ec.europa.eu/public_path/projects_data/focus/sw/docs/Generic\%20FOCUS_SWS_vc1.4.pdf}
+}
+\seealso{
+\link{drift_parameters_focus}, \link{PEC_sw_drift}
+}
diff --git a/man/drift_percentages_rautmann_formula.Rd b/man/drift_percentages_rautmann_formula.Rd
deleted file mode 100644
index eaf5c51..0000000
--- a/man/drift_percentages_rautmann_formula.Rd
+++ /dev/null
@@ -1,38 +0,0 @@
-% 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{
-# Compare JKI data with Rautmann formula
-# One application on field crops, for 1 m, 3 m and 5 m distance
-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")
-}
-\seealso{
-\link{drift_parameters_focus}, \link{PEC_sw_drift}
-}