hourlysis: Calculate hourly radiation from daily radiation

In ilyamaclean/UKCP18adjust: Correct UKCP18 to make it consistent with historic data

Description

Calculate hourly radiation from daily radiation

Usage

hourlysis(dailysis, modfit, rad_gam = NA, Trace = T)

Arguments

dailysis	an object of class spatialarray with hourly radiation values
modfit	a list of semivariagram model coefficients as returned by radfit()
rad_gam	an optional gam object of correction coefficients to apply to radiation data as returned by gamcorrect(). The default is for no correction to be applied.
Trace	logical value indicating whether to plot progress

Details

Hourly radiation values are computed by deriving daily optical depths from daily radiation values, interpolating these to hourly and back-calculating radiation for each hour. To accunt for sub-daily variation in cloud patchiness spatially and temporally autocorrelated cloud patchiness is simulated using paramaters from modfit.

Value

an object of class spatialarray containing the following components:

arraydata

a three-dimensional array of hourly radiation values

times

POSIXlt object of times associated with arraydata

crs

crs object of coordinate reference system associated with arraydata

extent

extent object giving extent covered by arraydata

units

units of arraydata

description

character description of arraydata

See Also

radfit()

Examples

# Takes a few seconds to run
modfit <- radfit(sis_2005)
# Takes a few minutes to run
ukcpsishourly <- hourlysis(sis_ukcpdaily, modfit)

ilyamaclean/UKCP18adjust documentation built on Nov. 4, 2019, 2:08 p.m.