julia subset data

subset.jl

using NCDataFrame
using NCDatasets
using Pkg
Pkg.add("IntervalSets")
using IntervalSets
using PlotlyJS

ds = NCDataset("ECMWF_ERA-40_subset.nc","r")


temp = ds["p2t"][:]


function selectnc(latmin,latmax,lonmin,lonmax,ncfile)
    nc = NCDataset(ncfile,"r")
    lats,lons = nc["latitude"][:],nc["longitude"][:]
    lonindices = findall((x -> x>(lonmin) && x<(lonmax)),lons)
    latindices = findall((y -> y>(latmin) && y<(latmax)),lats)
    subdata = nc["p2t"][lonindices,latindices,:]
    return subdata
end
 
dd = selectnc(4,25,77,99,"ECMWF_ERA-40_subset.nc")#


#plot(contour(x=ds["longitude"][:],y=ds["latitude"][:],z=temp[:,:,1]))
plot(contour(z=dd[:,:,1]))