You can install the development version from GitHub with:
# install.packages("devtools")
devtools::install_github("joaobtj/sw")
This is a basic example:
library(sw)
swb_franca_sp <- tibble::tribble(
~m, ~et, ~p,
"Jan", 117, 275,
"Feb", 102, 218,
"Mar", 104, 180,
"Apr", 79, 60,
"May", 60, 25,
"Jun", 49, 20,
"Jul", 54, 15,
"Aug", 74, 12,
"Sep", 93, 48,
"Oct", 107, 113,
"Nov", 108, 180,
"Dec", 117, 245
)
normal_swb(month = "m", ET0 = "et", P = "p", data = swb_franca_sp, awc = 125, start_month = "Apr")
#> m et p pmet0 L A DA ETa D
#> 1 Apr 79 60 -19 19.00000 107.37354 -17.626465 77.62646 1.373535
#> 2 May 60 25 -35 54.00000 81.15117 -26.222363 51.22236 8.777637
#> 3 Jun 49 20 -29 83.00000 64.34851 -16.802665 36.80266 12.197335
#> 4 Jul 54 15 -39 122.00000 47.10192 -17.246589 32.24659 21.753411
#> 5 Aug 74 12 -62 184.00000 28.68326 -18.418657 30.41866 43.581343
#> 6 Sep 93 48 -45 229.00000 20.01163 -8.671629 56.67163 36.328371
#> 7 Oct 107 113 6 196.22124 26.01163 6.000000 107.00000 0.000000
#> 8 Nov 108 180 72 30.40345 98.01163 72.000000 108.00000 0.000000
#> 9 Dec 117 245 128 0.00000 125.00000 26.988367 117.00000 0.000000
#> 10 Jan 117 275 158 0.00000 125.00000 0.000000 117.00000 0.000000
#> 11 Feb 102 218 116 0.00000 125.00000 0.000000 102.00000 0.000000
#> 12 Mar 104 180 76 0.00000 125.00000 0.000000 104.00000 0.000000
#> E
#> 1 0.0000
#> 2 0.0000
#> 3 0.0000
#> 4 0.0000
#> 5 0.0000
#> 6 0.0000
#> 7 0.0000
#> 8 0.0000
#> 9 101.0116
#> 10 158.0000
#> 11 116.0000
#> 12 76.0000