diff --git a/R/readNWISunit.r b/R/readNWISunit.r
index e79e55ccb15f81683bfb98960f0161891f050823..7d2fa9cf505ed02228c13831d8b177e1a61dfedc 100644
--- a/R/readNWISunit.r
+++ b/R/readNWISunit.r
@@ -24,12 +24,13 @@
 #' summary(rawData)
 #' rawData2 <- readNWISunit(siteNumber,parameterCd,startDate,endDate,"tsv")
 #' summary(rawData2)
-#' timeZoneChange <- readNWISunit(siteNumber,parameterCd,
+#' timeZoneChange <- readNWISunit(c('04024430','04024000'),parameterCd,
 #'          "2013-11-03","2013-11-03","tsv")
+#' firstSite <- timeZoneChange[timeZoneChange$site_no == '04024430',]
 readNWISunit <- function (siteNumber,parameterCd,startDate,endDate,format="xml"){  
   
   url <- constructNWISURL(siteNumber,parameterCd,startDate,endDate,"uv",format=format)
-  if (format == "xml") {
+  if (format == "xml" | format == "wml1") {
     data <- importWaterML1(url,asDateTime=TRUE)
   } else {
     data <- importRDB1(url,asDateTime=TRUE)
diff --git a/man/readNWISunit.Rd b/man/readNWISunit.Rd
index 75db70c38e1cd9c7ae49633018aa9e17974dc025..0aa9778e0ac6d74ea951556c096dabbfb6ca5aff 100644
--- a/man/readNWISunit.Rd
+++ b/man/readNWISunit.Rd
@@ -36,8 +36,9 @@ rawData <- readNWISunit(siteNumber,parameterCd,startDate,endDate)
 summary(rawData)
 rawData2 <- readNWISunit(siteNumber,parameterCd,startDate,endDate,"tsv")
 summary(rawData2)
-timeZoneChange <- readNWISunit(siteNumber,parameterCd,
+timeZoneChange <- readNWISunit(c('04024430','04024000'),parameterCd,
          "2013-11-03","2013-11-03","tsv")
+firstSite <- timeZoneChange[timeZoneChange$site_no == '04024430',]
 }
 \keyword{USGS}
 \keyword{data}