swdi distances seem broken

[sckott]

via #239

get_cordinates_df <- function(swdi_output) {
  swdi_output$shape %>%
    dplyr::mutate(shape = sub(".*POINT.*\\(","",shape)) %>%
    dplyr::mutate(shape = sub(")","",shape)) %>%
    tidyr::separate(shape, c("lon", "lat"), sep = " ") %>%
    dplyr::mutate(lat = as.numeric(lat), lon = as.numeric(lon))
}

center <- c(-102.0, 32.7)
out <- swdi(dataset='nx3tvs', startdate='20060506', enddate='20060507', 
    radius = 5, center = center, verbose = TRUE)

coordinates <- get_cordinates_df(out)
coordinates$distance <- meteo_spherical_distance(lat1=center[2], long1=center[1], 
	lat2=coordinates[2][[1]], long2=coordinates[1][[1]]) 
head(coordinates)

lon	lat	distance
-101.12693	34.48964	214.79856
-98.64874	36.41111	514.18312
-102.72588	33.20819	88.20470
-101.41145	32.89070	58.95814
-101.96767	33.13551	48.52061
-101.38741	32.88153	60.71622

distance values are in km, but still converting to miles means they are all much greater than 15 miles away from the focal point.

[sckott]

have emailed NOAA folks

[sckott]

heard back. radius/etc. stuff is not operational - they'll update docs soon, ideally wait till those are updated to update here.

[sckott]

docs at https://www.ncdc.noaa.gov/swdiws/