.onAttach <- function(libname, pkgname) {
  packageStartupMessage("This information is preliminary or provisional and is subject to revision. It is being provided to meet the need for timely best science. The information has not received final approval by the U.S. Geological Survey (USGS) and is provided on the condition that neither the USGS nor the U.S. Government shall be held liable for any damages resulting from the authorized or unauthorized use of the information.")
}

> fabric <- webdata(url="http://cida.usgs.gov/thredds/dodsC/mows/sr")
> stencil <- webgeom("state::CO")
> job <- geoknife(stencil, fabric) #ignore the checkAttrNamespaces warning for now - it's a known issue
Warning message:
In checkAttrNamespaces(getEffectiveNamespaces(node), .attrs, suppressNamespaceWarning) :
  missing namespace definitions for prefix(es) xlink

> check(job)
Error in match.arg(state, c("none", "ProcessStarted", "Process successful",  : 
  'arg' should be one of “none”, “ProcessStarted”, “Process successful”, “ProcessFailed”, “unknown”
> traceback()
5: stop(gettextf("'arg' should be one of %s", paste(dQuote(choices), 
       collapse = ", ")), domain = NA)
4: match.arg(state, c("none", "ProcessStarted", "Process successful", 
       "ProcessFailed", "unknown"))
3: setJobState(process$status)
2: check(job)
1: check(job)

Error in start(geojob) : 
  Cannot start a new geojob until a previous one is completed or is error. See "check(geojob)"

setWFS(gk) <- 'https://www.sciencebase.gov/catalogMaps/mapping/ows/50d35261e4b062c7914ebd14'
getShapefiles(gk)

# geoknife v1.0

# -- new geoknife classes -- 
geojob
webgeom
webdata
webprocess

# -creation methods for webdata
data <- webdata(list) # create and populate slots
data <- webdata(algorithm) # empty webdata w/ slots based on algorithm
data <- webdata(geojob) # get dataset object from previous job

# -webdata methods
query(webdata,'variables')
query(webdata) # get all metadata from dataset
query(webdata,'times') # get start and stop times of dataset
times(webdata)<- list #start_time = POSIXct,...





...
# -creation methodd for webgeom
webgeom <- webgeom(WFS = "cida...")
webgeom <- webgeom(geojob)

# -webgeom methods
query(webgeom,'shapes')
query(webgeom,'attributes')
query(webgeom,'IDs')
IDs(webgeom) <- 'character'
IDs(webgeom)
attribute(webgeom) <- 'character'
attribute(webgeom)
shape(webgeom) <- 'character'
shape(webgeom)


# geoknife would be the verb for the process starting ("slice this!"), creates geojob object
geojob <- geoknife(geom = "numeric", data = "character", excecute=TRUE, 
    knife = "Area Grid Statistics (weighted)"...) # as a quickstart
geojob <- geoknife(stencil = "sp::SpatialPolygonsDataFrame", fabric = dataset, excecute=TRUE, ...) # as a quickstart
geojob <- geoknife(stencil = "webgeom", fabric = "webdata", excecute=TRUE, ...) # as a quickstart
geojob <- geoknife(geom = "geojob", fabric = "character", excecute=TRUE, ...) # as a quickstart
geojob <- geoknife(stencil = "geojob", fabric = "character", knife = "webprocess", excecute = TRUE, ...) # as a quickstart
geojob <- geoknife(geojob, excecute = TRUE) # which just re-runs an old job or a job that was built w/ excecute=F

# -- supported input types that are handled w/ appropriate dispatch --
geom <- XML::XML
geom <- sp::SpatialPointsDataFrame
geom <- sp::SpatialPolygonsDataFrame
geom <- webgeom
geom <- numeric # linear ring of lat/lon
geom <- numeric # single lat/lon point
geom <- geojob # get geom from previous job
data <- "character" # for shortname quickstart for datasets lists that are built into the package (e.g., prism)



# -- methods for geojob -- 
algorithm(geojob)
algorithm(geojob)<-
geom(geojob)
geom(geojob)<-
webdata(geojob)
webdata(geojob)<-
summarize(geojob)
print(geojob)
check(geojob)
isRunning(geojob)
isError(geojob)
data.frame(geojob) # was loadOutput()
execute(geojob)

# -- functions in the geoknife package
getAlgorithms() # can specify WPS = "cida..."
getAlgorithms(webprocess)

job <- geoknife(stencil = c(-89,42), fabric = 'prism')
check(job)

#this job only takes a few seconds
sleep(10)
job <- geoknife(stencil = c(-89,42), fabric = 'prism')
check(job)

stationNodes = lapply(stationTriplets, xmlNode, name="stationTriplets")
  cmd = addChildren(cmd, kids=stationNodes)

wg <- webgeom()