scripts/step1_cal_HW_3con.R

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

library(ncdf4)
library(data.table)
library(lubridate)
library(dplyr)
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)

HW10 <- HW_3con_matrix(tmax_arr, date = date, k = 3, p_high = 0.975, p_low = 0.81)
HW11 <- HW_3con_matrix(tmax_arr, date = date, k = 3, p_high = 0.9, p_low = 0.75)

HW12 <- HW_3con_matrix(tmin_arr, date = date, k = 3, p_high = 0.975, p_low = 0.81)
HW13 <- HW_3con_matrix(tmin_arr, date = date, k = 3, p_high = 0.9, p_low = 0.75)