Calculate drift percentages and PECsw using the Rautmann formula

4 files changed, 160 insertions, 15 deletions

diff --git a/docs/reference/PEC_sw_drift.html b/docs/reference/PEC_sw_drift.html
index d4cd191..8c5bc3a 100644
--- a/docs/reference/PEC_sw_drift.html
+++ b/docs/reference/PEC_sw_drift.html
@@ -60,7 +60,10 @@ with input via spray drift.</p>
 <span>  water_depth <span class="op">=</span> <span class="fl">30</span>,</span>
 <span>  drift_percentages <span class="op">=</span> <span class="cn">NULL</span>,</span>
 <span>  drift_data <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"JKI"</span>, <span class="st">"RF"</span><span class="op">)</span>,</span>
-<span>  crop <span class="op">=</span> <span class="st">"Ackerbau"</span>,</span>
+<span>  crop_group_JKI <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"Ackerbau"</span>, <span class="st">"Obstbau frueh"</span>, <span class="st">"Obstbau spaet"</span>, <span class="st">"Weinbau frueh"</span>,</span>
+<span>    <span class="st">"Weinbau spaet"</span>, <span class="st">"Hopfenbau"</span>, <span class="st">"Flaechenkulturen &gt; 900 l/ha"</span>, <span class="st">"Gleisanlagen"</span><span class="op">)</span>,</span>
+<span>  crop_group_focus <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"arable"</span>, <span class="st">"hops"</span>, <span class="st">"vines, late"</span>, <span class="st">"vines, early"</span>, <span class="st">"fruit, late"</span>,</span>
+<span>    <span class="st">"fruit, early"</span>, <span class="st">"aerial"</span><span class="op">)</span>,</span>
 <span>  distances <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">5</span>, <span class="fl">10</span>, <span class="fl">20</span><span class="op">)</span>,</span>
 <span>  rate_units <span class="op">=</span> <span class="st">"g/ha"</span>,</span>
 <span>  PEC_units <span class="op">=</span> <span class="st">"µg/L"</span></span>
@@ -92,8 +95,13 @@ included in the package is used. If 'RF', the Rautmann formula is used, if
 implemented for the crop type and number of applications</p></dd>
 
 
-<dt>crop</dt>
-<dd><p>Crop name (use German names for JKI data), defaults to "Ackerbau"</p></dd>
+<dt>crop_group_JKI</dt>
+<dd><p>When using the 'JKI' drift data, one of the German names
+as used in <a href="drift_parameters_focus.html">drift_parameters_focus</a>.</p></dd>
+
+
+<dt>crop_group_focus</dt>
+<dd><p>One of the crop groups as used in <a href="drift_parameters_focus.html">drift_parameters_focus</a></p></dd>
 
 
 <dt>distances</dt>
@@ -127,14 +135,22 @@ implemented for the crop type and number of applications</p></dd>
 <span class="r-in"><span><span class="co"># Alternatively, we can use the formula for a single application to "Ackerbau" from the paper</span></span></span>
 <span class="r-in"><span><span class="fu">PEC_sw_drift</span><span class="op">(</span><span class="fl">100</span>, drift_data <span class="op">=</span> <span class="st">"RF"</span><span class="op">)</span></span></span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>        1 m        5 m       10 m       20 m </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> 0.92350000 0.19114149 0.09699222 0.04921742 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 0.91976667 0.19064473 0.09680051 0.04915079 </span>
 <span class="r-in"><span></span></span>
-<span class="r-in"><span><span class="co"># This makes it possible to also use different substances</span></span></span>
+<span class="r-in"><span><span class="co"># This makes it possible to also use different distances</span></span></span>
 <span class="r-in"><span><span class="fu">PEC_sw_drift</span><span class="op">(</span><span class="fl">100</span>, distances <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">6</span>, <span class="fl">10</span>, <span class="fl">20</span>, <span class="fl">50</span>, <span class="fl">100</span><span class="op">)</span>, drift_data <span class="op">=</span> <span class="st">"RF"</span><span class="op">)</span></span></span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>        1 m        3 m        5 m        6 m       10 m       20 m       50 m </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> 0.92350000 0.31512171 0.19114149 0.15990435 0.09699222 0.04921742 0.02007497 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 0.91976667 0.31415827 0.19064473 0.15951494 0.09680051 0.04915079 0.02006434 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>      100 m </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 0.01018774 </span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co"># or consider aerial application</span></span></span>
+<span class="r-in"><span><span class="fu">PEC_sw_drift</span><span class="op">(</span><span class="fl">100</span>, distances <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">6</span>, <span class="fl">10</span>, <span class="fl">20</span>, <span class="fl">50</span>, <span class="fl">100</span><span class="op">)</span>, drift_data <span class="op">=</span> <span class="st">"RF"</span>, </span></span>
+<span class="r-in"><span>  crop_group_focus <span class="op">=</span> <span class="st">"aerial"</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>        1 m        3 m        5 m        6 m       10 m       20 m       50 m </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 16.8233333 11.0585820  9.0986174  8.4866460  6.9825178  4.7004640  1.8820816 </span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>      100 m </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> 0.01018678 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>  0.9417586 </span>
 <span class="r-in"><span></span></span>
 <span class="r-in"><span><span class="co"># Using custom drift percentages is also supported</span></span></span>
 <span class="r-in"><span><span class="fu">PEC_sw_drift</span><span class="op">(</span><span class="fl">100</span>, drift_percentages <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">2.77</span>, <span class="fl">0.95</span>, <span class="fl">0.57</span>, <span class="fl">0.48</span>, <span class="fl">0.29</span>, <span class="fl">0.15</span>, <span class="fl">0.06</span>, <span class="fl">0.03</span><span class="op">)</span><span class="op">)</span></span></span>
diff --git a/docs/reference/drift_parameters_focus.html b/docs/reference/drift_parameters_focus.html
index 4c7539a..c3f6386 100644
--- a/docs/reference/drift_parameters_focus.html
+++ b/docs/reference/drift_parameters_focus.html
@@ -65,6 +65,13 @@ as its licence is not clear.</p>
     <p>A <a href="https://tibble.tidyverse.org/reference/tibble.html" class="external-link">tibble::tibble</a>.</p>
     </div>
     <div class="section level2">
+    <h2 id="details">Details<a class="anchor" aria-label="anchor" href="#details"></a></h2>
+    <p>For the hinge distance, <code>Inf</code> 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.</p>
+    </div>
+    <div class="section level2">
     <h2 id="references">References<a class="anchor" aria-label="anchor" href="#references"></a></h2>
     <p>FOCUS (2014) Generic guidance for Surface Water Scenarios (version 1.4).
 FOrum for the Co-ordination of pesticde fate models and their USe.
@@ -84,17 +91,20 @@ Bundesanst. Land- Forstwirtsch. 383, 133-141</p>
 <span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #949494;"># A tibble: 49 × 8</span></span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>    crop_group n_apps percentile     A      B     C     D hinge</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>    <span style="color: #949494; font-style: italic;">&lt;chr&gt;</span>       <span style="color: #949494; font-style: italic;">&lt;int&gt;</span>      <span style="color: #949494; font-style: italic;">&lt;int&gt;</span> <span style="color: #949494; font-style: italic;">&lt;dbl&gt;</span>  <span style="color: #949494; font-style: italic;">&lt;dbl&gt;</span> <span style="color: #949494; font-style: italic;">&lt;dbl&gt;</span> <span style="color: #949494; font-style: italic;">&lt;dbl&gt;</span> <span style="color: #949494; font-style: italic;">&lt;dbl&gt;</span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 1</span> arable          1         90  2.76 -<span style="color: #BB0000;">0.978</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 2</span> arable          2         82  2.44 -<span style="color: #BB0000;">1.01</span>    <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 3</span> arable          3         77  2.02 -<span style="color: #BB0000;">0.996</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 4</span> arable          4         74  1.86 -<span style="color: #BB0000;">0.986</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 5</span> arable          5         72  1.79 -<span style="color: #BB0000;">0.994</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 6</span> arable          6         70  1.63 -<span style="color: #BB0000;">0.986</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 7</span> arable          7         69  1.58 -<span style="color: #BB0000;">0.981</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 8</span> arable          8         67  1.51 -<span style="color: #BB0000;">0.983</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>     <span style="color: #BB0000;">NA</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 1</span> arable          1         90  2.76 -<span style="color: #BB0000;">0.978</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 2</span> arable          2         82  2.44 -<span style="color: #BB0000;">1.01</span>    <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 3</span> arable          3         77  2.02 -<span style="color: #BB0000;">0.996</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 4</span> arable          4         74  1.86 -<span style="color: #BB0000;">0.986</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 5</span> arable          5         72  1.79 -<span style="color: #BB0000;">0.994</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 6</span> arable          6         70  1.63 -<span style="color: #BB0000;">0.986</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 7</span> arable          7         69  1.58 -<span style="color: #BB0000;">0.981</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 8</span> arable          8         67  1.51 -<span style="color: #BB0000;">0.983</span>   <span style="color: #BB0000;">NA</span>  <span style="color: #BB0000;">NA</span>    <span style="color: #BB0000;">Inf</span>  </span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;"> 9</span> hops            1         90 58.2  -<span style="color: #BB0000;">1.00</span>  <span style="text-decoration: underline;">8</span>655. -<span style="color: #BB0000;">2.84</span>  15.3</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #BCBCBC;">10</span> hops            2         82 66.2  -<span style="color: #BB0000;">1.20</span>  <span style="text-decoration: underline;">5</span>555. -<span style="color: #BB0000;">2.82</span>  15.3</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> <span style="color: #949494;"># ℹ 39 more rows</span></span>
+<span class="r-in"><span><span class="fu"><a href="https://rdrr.io/r/base/unique.html" class="external-link">unique</a></span><span class="op">(</span><span class="va">drift_parameters_focus</span><span class="op">$</span><span class="va">crop_group</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] "arable"       "hops"         "vines, late"  "vines, early" "fruit, late" </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [6] "fruit, early" "aerial"      </span>
 </code></pre></div>
     </div>
   </main><aside class="col-md-3"><nav id="toc"><h2>On this page</h2>
diff --git a/docs/reference/drift_percentages_rautmann_formula.html b/docs/reference/drift_percentages_rautmann_formula.html
new file mode 100644
index 0000000..3d51fb0
--- /dev/null
+++ b/docs/reference/drift_percentages_rautmann_formula.html
@@ -0,0 +1,114 @@
+<!DOCTYPE html>
+<!-- Generated by pkgdown: do not edit by hand --><html lang="en"><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><meta charset="utf-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"><meta name="description" content="Calculate the drift percentages according to the Rautmann formula"><title>Calculate the drift percentages according to the Rautmann formula — drift_percentages_rautmann_formula • pfm</title><script src="../deps/jquery-3.6.0/jquery-3.6.0.min.js"></script><meta name="viewport" content="width=device-width, initial-scale=1, shrink-to-fit=no"><link href="../deps/bootstrap-5.3.1/bootstrap.min.css" rel="stylesheet"><script src="../deps/bootstrap-5.3.1/bootstrap.bundle.min.js"></script><!-- Font Awesome icons --><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/all.min.css" integrity="sha256-mmgLkCYLUQbXn0B1SRqzHar6dCnv9oZFPEC1g1cwlkk=" crossorigin="anonymous"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.12.1/css/v4-shims.min.css" integrity="sha256-wZjR52fzng1pJHwx4aV2AO3yyTOXrcDW7jBpJtTwVxw=" crossorigin="anonymous"><!-- bootstrap-toc --><script src="https://cdn.jsdelivr.net/gh/afeld/bootstrap-toc@v1.0.1/dist/bootstrap-toc.min.js" integrity="sha256-4veVQbu7//Lk5TSmc7YV48MxtMy98e26cf5MrgZYnwo=" crossorigin="anonymous"></script><!-- headroom.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/headroom.min.js" integrity="sha256-AsUX4SJE1+yuDu5+mAVzJbuYNPHj/WroHuZ8Ir/CkE0=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.11.0/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script><!-- clipboard.js --><script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.6/clipboard.min.js" integrity="sha256-inc5kl9MA1hkeYUt+EC3BhlIgyp/2jDIyBLS6k3UxPI=" crossorigin="anonymous"></script><!-- search --><script src="https://cdnjs.cloudflare.com/ajax/libs/fuse.js/6.4.6/fuse.js" integrity="sha512-zv6Ywkjyktsohkbp9bb45V6tEMoWhzFzXis+LrMehmJZZSys19Yxf1dopHx7WzIKxr5tK2dVcYmaCk2uqdjF4A==" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/autocomplete.js/0.38.0/autocomplete.jquery.min.js" integrity="sha512-GU9ayf+66Xx2TmpxqJpliWbT5PiGYxpaG8rfnBEk1LL8l1KGkRShhngwdXK1UgqhAzWpZHSiYPc09/NwDQIGyg==" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mark.js/8.11.1/mark.min.js" integrity="sha512-5CYOlHXGh6QpOFA/TeTylKLWfB3ftPsde7AnmhuitiTX4K5SqCLBeKro6sPS8ilsz1Q4NRx3v8Ko2IBiszzdww==" crossorigin="anonymous"></script><!-- pkgdown --><script src="../pkgdown.js"></script><meta property="og:title" content="Calculate the drift percentages according to the Rautmann formula — drift_percentages_rautmann_formula"><meta property="og:description" content="Calculate the drift percentages according to the Rautmann formula"><!-- mathjax --><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script><script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script><!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]--></head><body>
+    <a href="#main" class="visually-hidden-focusable">Skip to contents</a>
+    
+
+    <nav class="navbar fixed-top navbar-light navbar-expand-lg bg-light"><div class="container">
+    
+    <a class="navbar-brand me-2" href="../index.html">pfm</a>
+
+    <small class="nav-text text-muted me-auto" data-bs-toggle="tooltip" data-bs-placement="bottom" title="">0.6.1</small>
+
+    
+    <button class="navbar-toggler" type="button" data-bs-toggle="collapse" data-bs-target="#navbar" aria-controls="navbar" aria-expanded="false" aria-label="Toggle navigation">
+      <span class="navbar-toggler-icon"></span>
+    </button>
+
+    <div id="navbar" class="collapse navbar-collapse ms-3">
+      <ul class="navbar-nav me-auto"><li class="active nav-item">
+  <a class="nav-link" href="../reference/index.html">Reference</a>
+</li>
+      </ul><form class="form-inline my-2 my-lg-0" role="search">
+        <input type="search" class="form-control me-sm-2" aria-label="Toggle navigation" name="search-input" data-search-index="../search.json" id="search-input" placeholder="Search for" autocomplete="off"></form>
+
+      <ul class="navbar-nav"><li class="nav-item">
+  <a class="external-link nav-link" href="https://github.com/jranke/pfm/" aria-label="github">
+    <span class="fab fa fab fa-github fa-lg"></span>
+     
+  </a>
+</li>
+      </ul></div>
+
+    
+  </div>
+</nav><div class="container template-reference-topic">
+<div class="row">
+  <main id="main" class="col-md-9"><div class="page-header">
+      <img src="" class="logo" alt=""><h1>Calculate the drift percentages according to the Rautmann formula</h1>
+      <small class="dont-index">Source: <a href="https://github.com/jranke/pfm/blob/HEAD/R/PEC_sw_drift.R" class="external-link"><code>R/PEC_sw_drift.R</code></a></small>
+      <div class="d-none name"><code>drift_percentages_rautmann_formula.Rd</code></div>
+    </div>
+
+    <div class="ref-description section level2">
+    <p>Calculate the drift percentages according to the Rautmann formula</p>
+    </div>
+
+    <div class="section level2">
+    <h2 id="ref-usage">Usage<a class="anchor" aria-label="anchor" href="#ref-usage"></a></h2>
+    <div class="sourceCode"><pre class="sourceCode r"><code><span><span class="fu">drift_percentages_rautmann_formula</span><span class="op">(</span></span>
+<span>  <span class="va">distances</span>,</span>
+<span>  applications <span class="op">=</span> <span class="fl">1</span>,</span>
+<span>  crop_group_focus <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"arable"</span>, <span class="st">"hops"</span>, <span class="st">"vines, late"</span>, <span class="st">"vines, early"</span>, <span class="st">"fruit, late"</span>,</span>
+<span>    <span class="st">"fruit, early"</span>, <span class="st">"aerial"</span><span class="op">)</span></span>
+<span><span class="op">)</span></span></code></pre></div>
+    </div>
+
+    <div class="section level2">
+    <h2 id="arguments">Arguments<a class="anchor" aria-label="anchor" href="#arguments"></a></h2>
+    <dl><dt>distances</dt>
+<dd><p>The distances in m for which to get PEC values</p></dd>
+
+
+<dt>applications</dt>
+<dd><p>Number of applications for selection of drift percentile</p></dd>
+
+
+<dt>crop_group_focus</dt>
+<dd><p>One of the crop groups as used in <a href="drift_parameters_focus.html">drift_parameters_focus</a></p></dd>
+
+</dl></div>
+
+    <div class="section level2">
+    <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-examples"></a></h2>
+    <div class="sourceCode"><pre class="sourceCode r"><code><span class="r-in"><span><span class="va">drift_data_JKI</span><span class="op">[[</span><span class="fl">1</span><span class="op">]</span><span class="op">]</span><span class="op">[</span><span class="fu"><a href="https://rdrr.io/r/base/character.html" class="external-link">as.character</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">3</span>, <span class="fl">5</span><span class="op">)</span><span class="op">)</span>, <span class="st">"Ackerbau"</span><span class="op">]</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>    1    3    5 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 2.77 0.95 0.57 </span>
+<span class="r-in"><span><span class="fu">drift_percentages_rautmann_formula</span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">1</span>, <span class="fl">3</span>, <span class="fl">5</span><span class="op">)</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] 2.7593000 0.9424748 0.5719342</span>
+<span class="r-in"><span></span></span>
+<span class="r-in"><span><span class="co"># One application to early or late fruit crops</span></span></span>
+<span class="r-in"><span><span class="va">drift_data_JKI</span><span class="op">[[</span><span class="fl">1</span><span class="op">]</span><span class="op">]</span><span class="op">[</span><span class="fu"><a href="https://rdrr.io/r/base/character.html" class="external-link">as.character</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">20</span>, <span class="fl">50</span><span class="op">)</span><span class="op">)</span>, <span class="st">"Obstbau frueh"</span><span class="op">]</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>     3     5    20    50 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 29.20 19.89  2.77  0.30 </span>
+<span class="r-in"><span><span class="fu">drift_percentages_rautmann_formula</span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">20</span>, <span class="fl">50</span><span class="op">)</span>, crop_group <span class="op">=</span> <span class="st">"fruit, early"</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] 29.1973659 19.8844422  2.7618138  0.3012008</span>
+<span class="r-in"><span><span class="va">drift_data_JKI</span><span class="op">[[</span><span class="fl">1</span><span class="op">]</span><span class="op">]</span><span class="op">[</span><span class="fu"><a href="https://rdrr.io/r/base/character.html" class="external-link">as.character</a></span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">20</span>, <span class="fl">50</span><span class="op">)</span><span class="op">)</span>, <span class="st">"Obstbau spaet"</span><span class="op">]</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span>     3     5    20    50 </span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> 15.73  8.41  1.09  0.22 </span>
+<span class="r-in"><span><span class="fu">drift_percentages_rautmann_formula</span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="fl">3</span>, <span class="fl">5</span>, <span class="fl">20</span>, <span class="fl">50</span><span class="op">)</span>, crop_group <span class="op">=</span> <span class="st">"fruit, late"</span><span class="op">)</span></span></span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [1] 15.7246994  8.4108487  1.0813887  0.2155992</span>
+</code></pre></div>
+    </div>
+  </main><aside class="col-md-3"><nav id="toc"><h2>On this page</h2>
+    </nav></aside></div>
+
+
+    <footer><div class="pkgdown-footer-left">
+  <p></p><p>Developed by Johannes Ranke.</p>
+</div>
+
+<div class="pkgdown-footer-right">
+  <p></p><p>Site built with <a href="https://pkgdown.r-lib.org/" class="external-link">pkgdown</a> 2.0.7.</p>
+</div>
+
+    </footer></div>
+
+  
+
+  
+
+  </body></html>
+
diff --git a/docs/reference/index.html b/docs/reference/index.html
index 7f0b768..24f3012 100644
--- a/docs/reference/index.html
+++ b/docs/reference/index.html
@@ -197,6 +197,11 @@ published by the JKI</dd>
         <dd>Regression parameters for the Rautmann drift data</dd>
       </dl><dl><dt>
           
+          <code><a href="drift_percentages_rautmann_formula.html">drift_percentages_rautmann_formula()</a></code> 
+        </dt>
+        <dd>Calculate the drift percentages according to the Rautmann formula</dd>
+      </dl><dl><dt>
+          
           <code><a href="PEC_sw_drainage_UK.html">PEC_sw_drainage_UK()</a></code> 
         </dt>
         <dd>Calculate initial predicted environmental concentrations in surface water due to drainage using the UK method</dd>