# NOT RUN {
data(meteofrance)
## Choose a particular station
station_id <- unique(meteofrance$station_id)[3]
variables <- c("timestamp","DW","VW","TA","TD","RH","MFR_rr24",
"MFR_tn12","MFR_tn24","MFR_tx12","MFR_tx24","HS","HS_fresh")
header <- c("longitude","latitude","station_id" ,"altitude","location")
names(header) <- header
data <- meteofrance[meteofrance$station_id==station_id,c(header,variables)]
metaparam <- attr(meteofrance,"metaparam")
metaparam <- metaparam[metaparam$SMET_ID %in% names(data),]
header <- lapply(X=header,FUN=function(x,data) {data[1,x]},data=data)
data <- data[,variables]
attr(data,"header") <- header
attr(data,"metaparam") <- metaparam
sm <- as.smet(data)
# In case of multiple station, it return a list of SMET-class objects:
## This opton works for more than one station,
#in case of one station this works too only if it is FORCED!
sm_multi <- as.smet(meteofrance,variables=variables)
onestation <- meteofrance[meteofrance$station_id=="07499",]
sm_o <- as.smet(onestation,variables=variables,force.multistation=TRUE)
# }
Run the code above in your browser using DataLab