# ------------------------------------------------------------------
# LOAD DATA TO PRODUCE ALL FIGURES
# ------------------------------------------------------------------

# Sebastian Sippel
# 03.09.2021

library(raster)
library(ncdf4)
library(hydroGOF)
library(fields)
library(matrixStats)
library(bigmemory)
library(weights)

# ADJUST working directory to /figures directory:
setwd("/net/h2o/climphys1/sippels/_projects/low_freq_anchor_v2/figures/")


## 0.a. Load data to produce all figures:
## ------------------------------------------------------------------------------------------

# i. Load Anchor regression fingerprints and projections on model data:
load("../data/_cmip_anchor_LFIV/02_scen_cs.RData")
load("../data/_cmip_anchor_LFIV/02_piC_cs.RData")
load("../data/_cmip_anchor_LFIV/02_scen_all.RData")
# load("../data/_cmip_anchor_LFIV/02_piC_all.RData")

# ii. Load Observations projected onto anchor regression fingerprints:
load("../data/_final_processed/cv.OBS.RData")

# iii. Load model-specific trend distributions for each detection metrics:
load("../data/_final_processed/cv.piC.dist.RData")

# v. Load "maximum contribution of internal variability" at extremely likely level:
load("../data/_final_processed/frac.internal.95percent.RData")



## 0.b Load functions for figures:
## ------------------------------------------------------------------------------------------

source("../code/_convenience/frenchcolormap.R")
source("../code/_convenience/project_raster.R")
source("../code/_convenience/_plot_projected_worldmap_v2.R")