XML error when querying or downloading MACAv2MetData downscaled daily climate data via geoknife

[flakesw]

I'm looking to use geoknife to access the University of Idaho MACAv2MetData product. I've done so in the past without any issues, but I've recently encountered an error that, as far as I can tell, is only affecting the daily future MACAv2MetData product (https://cida.usgs.gov/thredds/dodsC/macav2metdata_daily_future). I was able to download it using the GeoData Portal, so I think it's not an issue with the data itself. The other MACAv2MetData products are working as expected (i.e., the daily historical product and monthly products), and so are all the other products I've tested.

> webdatasets = query('webdata')
> webdatasets[58]
An object of class "datagroup":
[1] Daily Future Multivariate Adaptive Constructed Analogs (MACA) CMIP5 Statistically Downscaled Data for Coterminous USA 
  url: https://cida.usgs.gov/thredds/dodsC/macav2metdata_daily_future 
> fabric <- webdata(webdatasets[58])
> query(fabric, 'variables')
Error in read_xml.raw(charToRaw(enc2utf8(x)), "UTF-8", ..., as_html = as_html,  : 
  Opening and ending tag mismatch: META line 2 and HEAD [76]

Here's the traceback:

13: read_xml.raw(charToRaw(enc2utf8(x)), "UTF-8", ..., as_html = as_html, 
        options = options)
12: read_xml.character(iconv(readBin(response$content, character()), 
        from = "UTF-8", to = "UTF-8"))
11: xml2::read_xml(iconv(readBin(response$content, character()), 
        from = "UTF-8", to = "UTF-8"))
10: gcontent(response)
9: xml2::xml_find_all(gcontent(response), "//gdp:shortname")
8: xml2::xml_text(xml2::xml_find_all(gcontent(response), "//gdp:shortname"), 
       )
7: variables_query(fabric, knife)
6: variables_query(fabric, knife)
5: variables_query(fabric = new("webdata", times = structure(c(NA_real_, 
   NA_real_), class = c("POSIXct", "POSIXt"), tzone = "UTC"), url = "https://cida.usgs.gov/thredds/dodsC/macav2metdata_daily_future", 
       variables = NA_character_, dataList = "gov.usgs.cida.gdp.wps.algorithm.discovery.ListOpendapGrids", 
       timeList = "gov.usgs.cida.gdp.wps.algorithm.discovery.GetGridTimeRange"))
4: variables_query(fabric = new("webdata", times = structure(c(NA_real_, 
   NA_real_), class = c("POSIXct", "POSIXt"), tzone = "UTC"), url = "https://cida.usgs.gov/thredds/dodsC/macav2metdata_daily_future", 
       variables = NA_character_, dataList = "gov.usgs.cida.gdp.wps.algorithm.discovery.ListOpendapGrids", 
       timeList = "gov.usgs.cida.gdp.wps.algorithm.discovery.GetGridTimeRange"))
3: do.call(paste0(field, "_query"), args)
2: query(fabric, "variables")
1: query(fabric, "variables")

I did recently update R, so perhaps it's an issue with a change in R or my local state, but I was able to replicate the error in both versions 4.2 and 4.3, and before and after cleaning my RStudio environment. I'm using geoknife 1.6.10. R version info below.

Thank for any advice on troubleshooting this error!

Sam

> R.Version()
$platform
[1] "x86_64-w64-mingw32"

$arch
[1] "x86_64"

$os
[1] "mingw32"

$crt
[1] "ucrt"

$system
[1] "x86_64, mingw32"

$status
[1] ""

$major
[1] "4"

$minor
[1] "3.0"

$year
[1] "2023"

$month
[1] "04"

$day
[1] "21"

$`svn rev`
[1] "84292"

$language
[1] "R"

$version.string
[1] "R version 4.3.0 (2023-04-21 ucrt)"

$nickname
[1] "Already Tomorrow"

[dblodgett-usgs]

Thanks for the report. I'm able to reproduce this.

Looking in a little further, something is timing out and failing. I'll need to do a bit more digging to see what might be causing the issue.

[flakesw]

Thanks for looking into it. No hurry on my end -- I was able to get everything I needed from the GeoData Portal.

[dblodgett-usgs]

Great. It may be the slow decline of the system, unfortunately. We are going to be retiring the GDP and migrating to other tools in the coming year. GeoKnife will have the same functionality but with a different back end. From: Sam Flake ***@***.***> Sent: Tuesday, June 6, 2023 7:11 AM To: DOI-USGS/geoknife ***@***.***> Cc: Blodgett, David L ***@***.***>; Comment ***@***.***> Subject: [EXTERNAL] Re: [DOI-USGS/geoknife] XML error when querying or downloading MACAv2MetData downscaled daily climate data via geoknife (Issue #434) This email has been received from outside of DOI - Use caution before clicking on links, opening attachments, or responding. Thanks for looking into it. No hurry on my end -- I was able to get everything I needed from the GeoData Portal. - Reply to this email directly, view it on GitHub<#434 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AALMOQ4HCNJSFHMUA2257B3XJ426DANCNFSM6AAAAAAY3B23LQ>. You are receiving this because you commented.Message ID: ***@***.******@***.***>>