# using 
using RCall
using GeoArrays
using Dates
using NCDatasets

R"""
getwd()
# load("dTmax.rda")
"""

file = "data-raw/CN05.1_Rhu_1961_2017_daily_025x025.nc"
gdalinfo(file)

Date(2014,1,29):Dates.Month(1):Date(2014,07,29)
ds = NCDataset(file)
ds.attrib


dates = ds["time"][:]
years = year.(dates)

year_grps = unique(years)
arr = ds

for year in year_grps
    println(year)
end
