scripts/step1_cal_HW_1con_movingSPL.R

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

library(ncdf4)
library(data.table)
library(lubridate)
library(tidyr)
library(slider)
library(purrr)

ncfile <- nc_open("/mnt/n/Research/cmip5/heatwave/data-raw/CN05.1_Tmin_1961_2018_daily_025x025.nc")
tmin_arr <- ncvar_get(ncfile, "tmin")
time <- (ncvar_get(ncfile, 'time')/24)
date <- as_date(time, origin = '1961-01-01')
nc_close(ncfile)

ncfile <- nc_open("/mnt/n/Research/cmip5/heatwave/data-raw/CN05.1_Tmax_1961_2018_daily_025x025.nc")
tmax_arr <- ncvar_get(ncfile, "tmax")
nc_close(ncfile)


HW8 <- HW_1con_movingSPL(tmax_arr, date = date, prob = 0.95)
HW9 <- HW_1con_movingSPL(tmin_arr, date = date, prob = 0.95)
save(HW8, file = 'OUTPUT/HW_indexes/nonHI/HW8.rda', compress = 'xz')
save(HW9, file = 'OUTPUT/HW_indexes/nonHI/HW9.rda', compress = 'xz')