Runs SPLASH at a single location for one day

run_one_day(lat, elv, n, y, wn, sf, tc, pn, kCw = 1.05, kWm = 150)

Arguments

lat

double, decimal degrees.

elv

double, elevation, m A.S.L. Default: \(0\).

n

double, day of year.

y

double, year. Default: \(0\).

wn

double, daily soil moisture content, mm (wn).

sf

double, fraction of sunshine hours. Default: \(1\).

tc

double, mean daily air temperature, degrees C. Default: \(23.0\).

pn

double, daily precipitation, mm/day.

kCw

double, supply constant, mm/hr. Default: \(1.05\) (Federer, 1982)

kWm

double, soil moisture capacity, mm. Default: \(150\) (Cramer-Prentice, 1988)

Value

List with the following components:

  • ho .......... daily solar irradiation, J/m2

  • hn .......... daily net radiation, J/m2

  • ppfd ........ daily PPFD, mol/m2

  • cond ........ daily condensation water, mm

  • eet ......... daily equilibrium ET, mm

  • pet ......... daily potential ET, mm

  • aet ......... daily actual ET, mm

  • wn .......... daily soil moisture, mm

  • ro .......... daily runoff, mm

References

Cramer, W. and Prentice, I.C., 1988. Simulation of regional soil moisture deficits on a European scale. Norsk Geografisk Tidsskrift - Norwegian Journal of Geography, 42(2-3), pp.149–151. doi: 10.1080/00291958808552193

Federer, C.A., 1982. Transpirational supply and demand: plant, soil, and atmospheric effects evaluated by simulation. Water Resources Research, 18(2), pp.355-362. doi: 10.1029/WR018i002p00355

Examples

soil <- run_one_day(lat = 37.7, elv = 142, n = 172, y = 2000, wn = 75, sf = 1, tc = 23, pn = 5)
#> [1] "in run_one_day: pn = 5"
cat(sprintf("Soil moisture (run one day):\n"))
#> Soil moisture (run one day):
cat(sprintf(" Ho: %0.6f J/m2\n", soil$ho))
#> Ho: 41646763.360234 J/m2
cat(sprintf(" Hn: %0.6f J/m2\n", soil$hn))
#> Hn: 21774952.713925 J/m2
cat(sprintf(" PPFD: %0.6f mol/m2\n", soil$ppfd))
#> PPFD: 62.042300 mol/m2
cat(sprintf(" EET: %0.6f mm/d\n", soil$eet))
#> EET: 6.405468 mm/d
cat(sprintf(" PET: %0.6f mm/d\n", soil$pet))
#> PET: 8.070889 mm/d
cat(sprintf(" AET: %0.6f mm/d\n", soil$aet))
#> AET: 5.748034 mm/d
cat(sprintf(" Cn: %0.6f mm/d\n", soil$cond))
#> Cn: 0.885192 mm/d
cat(sprintf(" Wn: %0.6f mm\n", soil$wn))
#> Wn: 75.137158 mm
cat(sprintf(" RO: %0.6f mm\n", soil$ro))
#> RO: 0.000000 mm