Contains not work

[kongdd]

The example is not work at L23:

using GeoData, ArchGDAL, NCDatasets, Plots, Statistics

geturl(url, filename=splitdir(url)[2]) = begin
    isfile(filename) || download(url, filename)
    filename
end


# Load some layers from NetCDF #############################################

ncurl = "https://www.unidata.ucar.edu/software/netcdf/examples/tos_O1_2001-2002.nc"
ncfilename = geturl(ncurl, "tos_O1_2001-2002.nc")
stack = NCDstack(ncfilename)

# Load the sea surface temperature layer
A = stack[:tos]

# Plot the 1st, 4th, 7th and 10th months 
A[Ti(1:3:12)] |> plot
savefig("tos_4.png")

# Plot the Australia region
A[Ti(Contains(DateTime360Day(2001, 01, 17))), Lat(Between(0, -50)), Lon(Between(100, 160))] |> plot

I got the error:

A[Ti(Contains(DateTime360Day(2001, 01, 17))), Lat(Between(0, -50)), Lon(Between(100, 160))]
ERROR: LoadError: ArgumentError: `Contains` has no meaning with `Points`. Use `Near`

https://github.com/rafaqz/GeoData.jl/blob/master/examples/examples.jl#L23

[rafaqz]

Oh right I should change the example.

The netcdf time dimension is Points by default now, not Intervals. Points cant contain values.

Instead use Near:

A[Ti(Near(DateTime360Day(2001, 01, 17))), Lat(Between(0, -50)), Lon(Between(100, 160))] |> plot

[kongdd]

Thank you!