scripts/step1_cal_HW_1con_allSPL.R

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

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

ncfile <- nc_open("/mnt/n/Research/cmip5/heatwave/data-raw/CN05.1_Tm_1961_2018_daily_025x025.nc")
tmean_arr <- ncvar_get(ncfile, "tm")
nc_close(ncfile)

HW1 <- HW_1con_allSPL(tmean_arr, date = date, k = 2, prob = 0.9)
HW2 <- HW_1con_allSPL(tmean_arr, date = date, k = 2, prob = 0.95)
HW3 <- HW_1con_allSPL(tmax_arr, date = date, k = 2, prob = 0.9)
HW4 <- HW_1con_allSPL(tmin_arr, date = date, k = 2, prob = 0.9)
HW5 <- HW_1con_allSPL(tmean_arr, date = date, k = 2, prob = 0.98)
HW6 <- HW_1con_allSPL(tmean_arr, date = date, k = 2, prob = 0.99)
HW7 <- HW_1con_allSPL(tmax_arr, date = date, k = 3, THR = 35)