View source: R/WATCHcreateWFDobs.R
WATCHcreateWFDobs | R Documentation |
Extracts data from WATCH WFD netCDF files and builds a CRHM .obs file of 3-hour data containing t
, ea
, u10
, and p
values. The values of t
, ea
and u
are 6-hourly, with NA
values inserted. The data are output at MST. The output values can be interpolated to hourly values using the function HourlyWATCHObs
. The windspeeds are at 10m, so they are denoted as u10
. Air temperatures are at 2m. The values for ea
are computed from the atmospheric pressure (at 10m) and the absolute humidity (at 2m).
WATCHcreateWFDobs( nc.location = "", startyear = 1901, endyear = 2001, lon = 0, lat = 0, houroffset = 0, obsFileName = "", quiet = TRUE, logfile = "" )
nc.location |
Required. A character string of the directory holding the WATCH WFD netCDf files. This is a file path WITHOUT a terminal slash, e.g. z:\WATCH\WFD. |
startyear |
Optional. Year to begin. Must be in the range |
endyear |
Optional. Year to end. Must be in the range |
lon |
Required. Decimal longitude to extract for. |
lat |
Required. Decimal latitude to extract for. |
houroffset |
Required. Number of hours that the local location is offset from UTC (GMT). Must be negative in the western hemisphere. For Mountain Standard Time, the offset is |
obsFileName |
Required. Name of the .obs file to be created. |
quiet |
Optional. Suppresses display of messages, except for errors. If you are calling this function in an R script, you will usually leave |
logfile |
Optional. Name of the file to be used for logging the action. Normally not used. |
If successful, returns the value TRUE
and writes the specified .obs file. Each month's data is written as it is created. If unsuccessful, returns the value FALSE
.
Kevin Shook
R code for conversion of air pressure and absolute humidity was taken from project PEcAn The Predictive Ecosystem Analyzer http://pecanproject.github.io. The source code is available at https://github.com/PecanProject/pecan/blob/master/modules/data.atmosphere/R/metutils.R.
WATCHcreateWFDEIobs
## Not run: location <- 'z:\data\WATCH\WFD' obsFileName <- 'VermilionWATCH_WFD.obs' WATCHcreateWFDobs(location, 1979, 2001, -111.9, 53.2, -7, obsFileName) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.