-
Notifications
You must be signed in to change notification settings - Fork 0
/
get_wscstation.R
63 lines (55 loc) · 1.96 KB
/
get_wscstation.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
#' get_wscstation
#'
#' Retrieves station information for an individual Water Survey of Canada site,
#' adds a text string at position 21 that combines key elements for a title.
#'
#' @param {stnID} {A Water Survey of Canada station number}
#' @param {stn} {default is HYDAT_list - station information from ECDataExplorer}
#'
#' @author Paul Whitfield <[email protected]>
#'
#' @return {array} with 21 elements
#' \itemize{
#' \item {Station} {StationID}
#' \item {StationName} {Station Name}
#' \item {HYDStatus} {Active or Discontinued}
#' \item {Prov} {Province}
#' \item {Latitude}{}
#' \item {Longitude}{}
#' \item {DrainageArea} {km2}
#' \item {Years} {# of years with data}
#' \item {From} {Start Year}
#' \item {To} {End Year}
#' \item {Reg.} {Regulated }
#' \item {Flow} {if TRUE/Yes}
#' \item {Level} {if TRUE/Yes}
#' \item {Sed} {if TRUE/Yes}
#' \item {OperSched} {Continuous or Seasonal}
#' \item {RealTime} {if TRUE/Yes}
#' \item {RHBN} {if TRUE/Yes is in the reference hydrologic basin network}
#' \item {Region} {if TRUE/Yes is in the reference hydrologic basin network}
#' \item {Datum} {if TRUE/Yes is in the reference hydrologic basin network}
#' \item {Operator} {if TRUE/Yes is in the reference hydrologic basin network}
#' \item {Station_lname} {Added field combines ID,Name,Province and if RHBN an * is added}
#' }
#'
#' @export
#'
#' @examples
#' data(HYDAT_list)
#' s_info <- get_wscstation("05BB001")
#' title <- s_info[21]
#' print(title)
#'
get_wscstation <- function (stnID, stn=HYDAT_list) {
rhbn <-NULL
stninfo <-stn[stn$Station==stnID,]
if(length(stninfo[,1])==0) {
print(paste("WSC Station ",stnID," not found"))
return(stnID)
}
if(stninfo$RHBN==TRUE) (rhbn="*" )
stninfo[21] <-paste(stninfo$Station," - ",stninfo$StationName, " - ", stninfo$Prov, rhbn, sep="")
names (stninfo) [21] <-"Station_lname"
return (stninfo)
}