diff --git a/analysis/README.md b/analysis/README.md index e69de29..269419e 100644 --- a/analysis/README.md +++ b/analysis/README.md @@ -0,0 +1,5 @@ +# cwd_global: Apply CWD algorithm to global files + +## Workflow + +1. `make_tidy_cmip6.R`: Make original global files tidy diff --git a/analysis/make_tidy_cmip6.R b/analysis/make_tidy_cmip6.R index 09ccc90..80692da 100644 --- a/analysis/make_tidy_cmip6.R +++ b/analysis/make_tidy_cmip6.R @@ -3,10 +3,16 @@ library(dplyr) library(ggplot2) # list demo file path -path <- "/data/scratch/CMIP6ng/cmip6-ng/pr/day/native/" +path_cmip6 <- "~/data/cmip6-ng/" -## Precipitation --------------- -files <- list.files(path, pattern = "pr_day", full.names = TRUE) +## Evapotranspiration +varnam <- "evspsbl" +res <- "mon" +files <- list.files( + paste0(path_cmip6, varnam, "/", res, "/native/"), + pattern = ".nc", + full.names = TRUE + ) # load and convert df <- map2tidy( @@ -21,4 +27,9 @@ df <- map2tidy( fileprefix = "pr_day_CESM2_historical_r1i1p1f1_native" ) +## Precipitation --------------- +files <- list.files(path, pattern = "pr_day", full.names = TRUE) + + +