rmendels
commented
7 years ago Thanks for the feedback.
The reason I said you will be living dangerously is that there are so many possible variations when we go to any GRID file on any ERDDAP, that there is no way in the world I can possibly test everything. I have slowly been testing it on a variety of datasets, but there is a reason I warn people.
One possibility is the hycom went away for awhile from the ERDDAP, which means that all of a sudden you are reading non-existent data. The other is you may not be getting what you think because the function finds the "nearest" actual grid point, and there may be ties in the nearest.
In the example below, which is the value where it bombs? I will try to get into the code and reproduce the problem if you can narrow down the values that failed.
-Roy
PS - But getting 23,000 points ain't half bad. As I was saying to Lynn, what gets me is it is the same code just repeating as you feed in values.
On Jan 24, 2017, at 9:21 AM, Abram Fleishman notifications@github.com wrote:
I am getting an error when using the rxtracto(). We have successfully been using it to access the hycom_GLBa008_tyx dataset (we have extracted data for >23,000 points), however, halfway through we began to receive the following error:
Error in object[-omit, , drop = FALSE] : incorrect number of dimensions
Running through the function line by line I have traced it back to lines 208-215 of the rxtracto() source code. It is returning an NA when running line 209: paramdata <- ncdf4::ncvar_get(ncFile, parameter)
When looking at the output of line 208, it is saying that the date range and variable is coming from experiment: 91.2, however, when I look online it looks like it should be coming from experiment: 91.1 according to the date 2015-12-11 18:58:34 in example below. Any chance you have encountered this problem before? Thoughts?
library(rerddapXtracto) library(rerddap)
sample of data we are using
dataSub<-structure( list( Latitude = c( 18.411898, 18.412279, 18.413588, 18.4138915454545, 18.4141950909091, 18.4144986363636 ), Longitude = c( -104.210121, -104.21064, -104.211845, -104.212580272727, -104.213315545455, -104.214050818182 ), date = structure( c( 1449860314, 1449860254, 1449860194, 1449860134, 1449860074, 1449860014 ), class = c("POSIXct", "POSIXt"), tzone = "UTC" ), lon_hycom = c( 2270.8734875, 2270.867, 2270.8519375, 2270.84274659091, 2270.83355568182, 2270.82436477273 ), lat_hycom = c( 1898.32844792124, 1898.33590025503, 1898.36150420497, 1898.36744153318, 1898.37337886139, 1898.37931618959 ) ), .Names = c("Latitude", "Longitude", "date", "lon_hycom", "lat_hycom"), row.names = c(NA, -6L), class = "data.frame" )
input for rxtracto
xpos <- dataSub$lon_hycom ypos <- dataSub$lat_hycom tpos <- dataSub$date
sshInfo <- rerddap::info('hycom_GLBa008_tyx')
ssh <- rxtracto(dataInfo = sshInfo, parameter = 'ssh', xcoord = xpos, ycoord = ypos, tcoord = tpos, xlen = 0.1, ylen = 0.1,xName="X",yName = "Y", tName = "time",verbose = TRUE)
Error in object[-omit, , drop = FALSE] : incorrect number of dimensions
— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or mute the thread.
"The contents of this message do not reflect any position of the U.S. Government or NOAA."
Roy Mendelssohn Supervisory Operations Research Analyst NOAA/NMFS Environmental Research Division Southwest Fisheries Science Center Note new street address 110 McAllister Way Santa Cruz, CA 95060 Phone: (831)-420-3666 Fax: (831) 420-3980 e-mail: Roy.Mendelssohn@noaa.gov www: http://www.pfeg.noaa.gov/
"Old age and treachery will overcome youth and skill." "From those who have been given much, much will be expected" "the arc of the moral universe is long, but it bends toward justice" -MLK Jr.
abfleishman
I am getting an error when using the
rxtracto(). We have successfully been using it to access thehycom_GLBa008_tyxdataset (we have extracted data for >23,000 points), however, halfway through we began to receive the following error:Running through the function line by line I have traced it back to lines 208-215 of the
rxtracto()source code. It is returning anNAwhen running line 209:paramdata <- ncdf4::ncvar_get(ncFile, parameter)When looking at the output of line 208, it is saying that the date range and variable is coming from
experiment: 91.2, however, when I look online it looks like it should be coming fromexperiment: 91.1according to the date2015-12-11 18:58:34in example below.Any chance you have encountered this problem before? Thoughts?