scripts/spatiotemporal-cluster/s1_check_clusterIds.R

In CUG-hydro/heatwave: What the Package Does (One Line, Title Case)

source("scripts/main_pkgs.R")
devtools::load_all() # shortcut: ctrl+shift+L

# polyline
file_arc <- "data-raw/china/china_arc_sml_01deg.shp"
sp_arc_CH <- read_polyline(file_arc)

file = "OUTPUT/clusterId/HW7.nc"
time = nc_date(file)
lon = ncread(file, "lon")
lat = ncread(file, "lat")
info = nc_info(file)

year = 2010
period = c("2010-06-01", "2010-08-31")
arr <- ncread(file, 1, DatePeriod = period)$data[[1]]

range = c(69.75, 140.25, 14.75, 55.25)
write_fig(
    plot.cluster(arr, times = 1:92, origin = period[1], range = range,
                 main = "HW7",
                 sp.layout = sp_arc_CH),
    "HW7-2010.pdf", 10, 15
)