ERAhourlyLongwave: Estimates hourly incoming longwave radiation from ERA-Interim...

In CentreForHydrology/Reanalysis: Reads reanalysis data for CRHM

Estimates hourly incoming longwave radiation from ERA-Interim 3-hourly values or ERA-40 6-hourly values

Description

This function is called after using the function getNearestERAtimeseries. This function deaccumulates the ERA-Interim 12-hour cumulative values,using the function ERAdeaccum and interpolates the 3-hour values to hourly values, based on the extra-terrestrial hourly radiation. ERA-40 6-hourly instantaneous values are interpolated directly to hourly values.

Usage

ERAhourlyLongwave(
  ERAstrd,
  strdColnum = 1,
  ERAt2m,
  t2mColnum = 1,
  method = "linear",
  quiet = TRUE,
  logfile = ""
)

Arguments

ERAstrd	Required. The CRHMr obs dataframe of ERA strd values. The values must not be deaccumulated, as the ERAdeaccum function is called by this function.
strdColnum	Optional. The column number containing the strd values, not including the datetime. Default is column 1.
ERAt2m	Required. The CRHMr obs dataframe of ERA t2m values. The values must not be deaccumulated, as the deaccumERA function is called by this function.
t2mColnum	Optional. The column number containing the t2m values, not including the datetime. Default is column 1.
method	Optional. The methods to be used for interpolation of the air temperature. Currently supported methods are linear and spline. The default is to use linear interpolation.
quiet	Optional. Suppresses display of messages, except for errors. If you are calling this function in an R script, you will usually leave quiet=TRUE (i.e. the default). If you are working interactively, you will probably want to set quiet=FALSE.
logfile	Optional. Name of the file to be used for logging the action. Normally not used

Value

If successful, returns an obs dataframe containing the interpolated hourly ERA incoming longwave radiation (Qli) in W/m^2, and the hourly 2m air temperatures (t) in C. If unsuccessful, returns the value FALSE.

Author(s)

Kevin Shook

References

The value of the Stefan-Boltzmann constant was obtained from http://physics.nist.gov/cgi-bin/cuu/Value?sigma

See Also

ERAgetNearestTimeseries ERAdeaccum ERAhourlyShortwave

Examples

## Not run: 
hourlyQli <- ERAhourlyLongwave(ERAstrd=strd, ERAt2m=t2m)
## End(Not run)

CentreForHydrology/Reanalysis documentation built on Nov. 23, 2022, 2:57 a.m.