GetCRN: Get the US Climate Reference Network data from the NCEP ftp...

In NCAR/rwrfhydro: R tools for the WRF Hydro Model

Description

GetCRN Given the information provided to the function will download the data from NCEP ftp serverr and then subset the data to meet the criteria defined such as being in a domian or being in specific time period. Return a dataframe containing.

Usage

GetCRN(gageNames = NULL, siteIds = NULL, elements = NULL,
  timestep = "hourly", domain = NULL, startDate = Sys.time() - 10 *
  24 * 3600, endDate = Sys.time(), realTime = FALSE,
  parallel = FALSE)

Arguments

gageNames	Character vector of the gage names which is used in the file names (or any part of the file name)
siteIds	Characcter vector of the stations IDs so called WBANNO.
elements	Character vector containing all elements required from the data.
timestep	Character (DEFAULT = "hourly"). It can be subhourly, hourly, daily and monthly
domain	data.frame, containing four columns named lonmax, latmin, latmax, latmin. It will be used to subset the data to only those gauges falling into the rectangle domain
startDate	POSIXct, ending date
endDate	POSIXct, strating date
realTime	Logical (DEFAULT = FALSE)
parallel	Logical (DEFAULT = FALSE)

Details

The contiguous U.S. network of 114 stations was completed in 2008. There are two USCRN stations in Hawaii and deployment of a network of 29 stations in Alaska continues. The vision of the USCRN program is to maintain a sustainable high-quality climate observation network.

These data are provided in 4 different temporal resolution (subhourly, hourly, daily and monthly), and depending on the temporal resolution, the variables provided changes. Below is a summary of all the variables in each temporal resolution and their units.

subhourly data contains:

• chech this README.txt ftp://ftp.ncdc.noaa.gov/pub/data/uscrn/products/subhourly01/README.txt

• AIR_TEMPERATURE Celsius

• PRECIPITATION (mm)

• SOLAR_RADIATION (W/m^2)

• SR_FLAG (X)

• SURFACE_TEMPERATURE (Celsius)

• ST_TYPE (X)

• ST_FLAG (X)

• RELATIVE_HUMIDITY (%)

• RH_FLAG (X)

• SOIL_MOISTURE_5 (m^3/m^3)

• SOIL_TEMPERATURE_5 (Celsius)

• WETNESS (Ohms)

• WET_FLAG (X)

• WIND_1_5 (m/s)

• WIND_FLAG (X)

hourly data contains:

• chech this README.txt ftp://ftp.ncdc.noaa.gov/pub/data/uscrn/products/hourly02/README.txt

• T_CALC (Celsius)

• T_HR_AVG (Celsius)

• T_MAX (Celsius)

• T_MIN (Celsius)

• P_CALC (mm)

• SOLARAD (W/m^2)

• SOLARAD_FLAG (X)

• SOLARAD_MAX (W/m^2)

• SOLARAD_MAX_FLAG (X)

• SOLARAD_MIN (W/m^2)

• SOLARAD_MIN_FLAG (X)

• SUR_TEMP_TYPE (X)

• SUR_TEMP (Celsius)

• SUR_TEMP_FLAG (X)

• SUR_TEMP_MAX (Celsius)

• SUR_TEMP_MAX_FLAG (X)

• SUR_TEMP_MIN (Celsius)

• SUR_TEMP_MIN_FLAG (X)

• RH_HR_AVG (%)

• RH_HR_AVG_FLAG (X)

• SOIL_MOISTURE_5 (m^3/m^3)

• SOIL_MOISTURE_10 (m^3/m^3)

• SOIL_MOISTURE_20 (m^3/m^3)

• SOIL_MOISTURE_50 (m^3/m^3)

• SOIL_MOISTURE_100 (m^3/m^3)

• SOIL_TEMP_5 (Celsius)

• SOIL_TEMP_10 (Celsius)

• SOIL_TEMP_20 (Celsius)

• SOIL_TEMP_50 (Celsius)

• SOIL_TEMP_100 (Celsius)

daily data contains:

• chech this README.txt ftp://ftp.ncdc.noaa.gov/pub/data/uscrn/products/daily01/README.txt

• T_DAILY_MAX (Celsius)

• T_DAILY_MIN (Celsius)

• T_DAILY_MEAN (Celsius)

• T_DAILY_AVG (Celsius)

• P_DAILY_CALC (mm)

• SOLARAD_DAILY (MJ/m^2)

• SUR_TEMP_DAILY_TYPE (X)

• SUR_TEMP_DAILY_MAX (Celsius)

• SUR_TEMP_DAILY_MIN (Celsius)

• SUR_TEMP_DAILY_AVG (Celsius)

• RH_DAILY_MAX (%)

• RH_DAILY_MIN (%)

• RH_DAILY_AVG (%)

• SOIL_MOISTURE_5_DAILY (m^3/m^3)

• SOIL_MOISTURE_10_DAILY (m^3/m^3)

• SOIL_MOISTURE_20_DAILY (m^3/m^3)

• SOIL_MOISTURE_50_DAILY (m^3/m^3)

• SOIL_MOISTURE_100_DAILY (m^3/m^3)

• SOIL_TEMP_5_DAILY (Celsius)

• SOIL_TEMP_10_DAILY (Celsius)

• SOIL_TEMP_20_DAILY (Celsius)

• SOIL_TEMP_50_DAILY (Celsius)

• SOIL_TEMP_100_DAILY (Celsius)

monthly data contains:

• chech this README.txt ftp://ftp.ncdc.noaa.gov/pub/data/uscrn/products/monthly01/README.txt

• T_MONTHLY_MAX (Celsius)

• T_MONTHLY_MIN (Celsius)

• T_MONTHLY_MEAN (Celsius)

• T_MONTHLY_AVG (Celsius)

• P_MONTHLY_CALC (mm)

• SOLRAD_MONTHLY_AVG (MJ/m^2)

• SUR_TEMP_MONTHLY_TYPE (X)

• SUR_TEMP_MONTHLY_MAX (Celsius)

• SUR_TEMP_MONTHLY_MIN (Celsius)

• SUR_TEMP_MONTHLY_AVG (Celsius)

Value

A data.table containing the US CRN data.

Examples

## Not run: 
# example 1: getting the precipition data at the subhourly time steps, 
# for the first three months of 2016
startDate = as.POSIXct("2016-01-01 00:00:00", format= "%Y-%m-%d %H:%M:%S", tz = "UTC")
endDate   = as.POSIXct("2016-04-01 00:00:00", format= "%Y-%m-%d %H:%M:%S", tz = "UTC")
output <- GetCRN(elements = "P_CALC", timestep = "hourly", 
startDate = startDate, endDate= endDate)

# choose only gauge in the following domain
domain <- data.frame(minlon = -80, maxlon = -60, minlat = 30, maxlat = 50)
output <- GetCRN(elements = "P_CALC", timestep = "hourly", 
startDate = startDate, endDate= endDate, domain = domain)

# choose only gauges in the Colorado 
outputHourly <- GetCRN(gageNames = "CO", elements = "P_CALC", timestep = "hourly", 
startDate = startDate, endDate= endDate)

# get daily data only for basins in colorado
outputDaily <- GetCRN(gageNames = "CO", elements = "P_DAILY_CALC", timestep = "daily", 
startDate = startDate, endDate= endDate)

# get monthly data only for basins in colorado
outputMonthly <- GetCRN(gageNames = "CO", elements = "P_MONTHLY_CALC", timestep = "monthly", 
startDate = startDate, endDate= endDate)

# get subhourly data only for basins in colorado
outputsubhourly <- GetCRN(gageNames = "CO", elements = "PRECIPITATION", 
timestep = "subhourly", startDate = startDate, endDate= endDate)

# get hourly soik moisure data
outputHourly <- GetCRN(gageNames = "CO", elements = c("SOIL_MOISTURE_5", 
"SOIL_MOISTURE_10", "SOIL_MOISTURE_20", "SOIL_MOISTURE_50", "SOIL_MOISTURE_100"), 
timestep = "hourly")

## End(Not run)

NCAR/rwrfhydro documentation built on Feb. 28, 2021, 12:47 p.m.