weach: Simulates the hourly conditions from daily

In serbinsh/biocro: Simulation of Energycane and Sugarcane

Description

Manipulates weather data in the format obtained from WARM (see link below) and returns the format and units needed for most functions in this package. This function should be used for one year at a time. It returns hourly (or sub-daily) weather information.

Usage

  weach(X, lati, ts = 1,
    temp.units = c("Farenheit", "Celsius"),
    rh.units = c("percent", "fraction"),
    ws.units = c("mph", "mps"), pp.units = c("in", "mm"),
    ...)

Arguments

X	a matrix (or data frame) containing weather information. The input format is strict but it is meant to be used with the data usually obtained from weather stations in Illinois. The data frame should have 11 columns (see details).
lati	latitude at the specific location
ts	timestep for the simulation of sub-daily data from daily. For example a value of 3 would return data every 3 hours. Only divisors of 24 work (i.e. 1,2,3,4, etc.).
temp.units	Option to specify the units in which the temperature is entered. Default is Farenheit.
rh.units	Option to specify the units in which the relative humidity is entered. Default is percent.
ws.units	Option to specify the units in which the wind speed is entered. Default is miles per hour.
pp.units	Option to specify the units in which the precipitation is entered. Default is inches.
list()	additional arguments to be passed to lightME

Details

This function was originally used to transform daily data to hourly data. Some flexibility has been added so that other units can be used. The input data used originally looked as follows.

\itemcol

1year \itemcol 2day of the year (1–365). Does not consider leap years. \itemcol 3total daily solar radiation (MJ/m^2). \itemcol 4maximum temperature (Fahrenheit). \itemcol 5minimum temperature (Fahrenheit). \itemcol 6average temperature (Fahrenheit). \itemcol 7maximum relative humidity (%). \itemcol 8minimum relative humidity (%). \itemcol 9average relative humidity (%). \itemcol 10average wind speed (miles per hour). \itemcol 11precipitation (inches).

All the units above are the defaults but they can be changed as part of the arguments.

Value

a matrix returning hourly (or sub-daily) weather data. Dimensions 8760 (if hourly) by 8.

Examples

data(cmi0506)
tmp1 <- cmi0506[cmi0506$year == 2005,]
wet05 <- weach(tmp1,40)

# Return data every 3 hours
wet05.3 <- weach(tmp1,40,ts=3)

serbinsh/biocro documentation built on May 29, 2019, 6:57 p.m.