Saving a subset of a file

[nicholasbalasus]

Not a bug but a question... Is there a simple way to save a subset of a file? For example, how to save ds_subset as a new netCDF here (or modify the existing in-place):

using NCDatasets

src = "https://www.unidata.ucar.edu/software/netcdf/examples/tos_O1_2001-2002.nc"
dst = "tos_O1_2001-2002.nc"
download(src, dst)

ds = NCDataset("tos_O1_2001-2002.nc", "r")
ds_subset = NCDatasets.@select(ds, -10 <= lon <= 10 && -10 <= lat <= 10)

[nicholasbalasus]

I can't get it to work with the example above, but the following works for my specific netCDF file so I'll close this for now (based on this).

using NCDatasets

src = "http://geoschemdata.wustl.edu/ExtData/GEOS_0.5x0.625/MERRA2/2018/05/MERRA2.20180501.A1.05x0625.nc4"
dst = "MERRA2.20180501.A1.05x0625.nc4"
download(src, dst)

ds_all = NCDataset(dst, "r")
ds_afr = NCDataset(replace(dst, ".nc4" => ".crop.nc4"), "c")
write(ds_afr, NCDatasets.@select(ds_all,
      -21.25 <= lon <= 55.0 && -38.5 <= lat <= 42.0))
close(ds_afr)
close(ds_all)