khufkens

git checkout --orphan gh-pages

# preview files to be deleted
git rm -rf --dry-run .

# actually delete the files
git rm -rf .

khufkens

#' Format and run p-model drivers on Euler
#'
#' Takes site information for a single
#' or multiple sites and grabs all data
#' required for a p-model run in rsofun.
#'
#' Parameter settings are provided as
#' arguments, but could be altered after
#' the fact if desired.
#'

khufkens

library(terra)
library(stringr)

#---- static code ----

# read in the original data
r <- rast("LPX-Bern_SH1_gpppft.nc")

# first subset to only the required years
# 1978 - 2020

khufkens

find . -name '*.zip' -exec sh -c 'unzip {};'

khufkens

gdal_translate -a_ullr -180 90 180 -90 -a_srs EPSG:4326 world.topo.bathy.200408.3x5400x2700.jpg tmp.tif

khufkens

# Interpolation using a loess smoother
# and optimal span determination

library(MODISTools)
library(phenocamr)
library(tidyverse)

VI <- mt_subset(
  product = "MOD13Q1",
  lon = -72.1755,

khufkens

library(rvest)

# base html url
url <- "http://dapds00.nci.org.au/thredds/fileServer/ks32/CLEX_Data/PLUMBER2/v1-0/Flux/"

# catalogue
html_page = read_html("https://dap.nci.org.au/thredds/remoteCatalogService?catalog=http://dapds00.nci.org.au/thredds/catalog/ks32/CLEX_Data/PLUMBER2/v1-0/Flux/catalog.xml")

# grab files from the catalogue
files = html_page %>%

khufkens

# convert pngs to mp4
ffmpeg -r 2 -i tmp/%03d.png -y  out.mp4

# convert mp4 to twitter compatible format
ffmpeg -i out.mp4  -vf scale=1280:720 -b:v 5000K -minrate 5000K -maxrate 5000K -b:a 128K -r 30 -f mp4 -vcodec libx264 -profile:v high -pix_fmt yuv420p -strict -2  -ac 2 -acodec aac twitter.mp4

khufkens

# data wrangling
library(tidyverse)

# waffle plot
library(waffle)

# combining graphs
library(patchwork)

# text magic

khufkens

# Downsample HWSD data to the scale of 0.5 degree
# WFDEI / CRU data

library(raster)
library(hwsdr)
library(tidyverse)

# get mu_global at 0.5 degrees from the
# ORNL DAAC (which is easier than reading in
# that monstrous BIL file)