grid.warm: Grid warm rain statistics

In jmuelmen/sciadv.aaz6433: Tools to process precipitation producs

Description

Grid warm rain statistics

Usage

grid.warm(df = get.cloud.summaries(), bins_lon = 180, bins_lat = 90,
  out.fname = "cloud-warm-grid.rds", precip.classes = 1:3)

Arguments

df	Data frame to grid, get.cloud.summaries() by default
bins_lon	Number of bins or vector of bin boundaries (longitude)
bins_lat	Number of bins or vector of bin boundaries (latitude)
out.fname	File name for output

Examples

## make the \code{warmrain.rda} file used by \code{data(\link{warmrain})}
library(devtools)
warmrain <- grid.warm(get.cloud.summaries(),
                      bins_lon = 180, bins_lat = 90,
                      NULL) 
str(warmrain)
devtools::use_data(warmrain, pkg = "/home/jmuelmen/lib64/R/library/preciptools")


## look at the Great Lakes region, if that's your thing
df <- grid.warm(get.cloud.summaries() %>%
               filter(lon > -95, lon < -75, lat > 40, lat < 50),
               seq(-95, -75, 0.1), seq(40, 50, 0.1),
               "great-lakes-grid.rds")
col.frac <- rev(RColorBrewer::brewer.pal(9, "RdYlBu"))[c(1:5,7,9)]
phase.breaks <- c(0, 0.05, 0.1, 0.2 * 1:5)
df %>%
    mutate(sum = ice + liq + mixed) %>%
    gather(phase, frac, c(ice:mixed, sum)) %>%
    discretize(frac, phase.breaks, TRUE) %>%
    ggplot(aes(x = lon, y = lat)) +
    geom_raster(aes(fill = frac)) +
    geom_world_polygon(highres = TRUE, lwd = 0.25) +
    xlim(c(-95, -75)) + ylim(c(40,50)) +
    scale_fill_manual(values = col.frac, drop = FALSE) +
    facet_grid(phase ~ denom) +
    coord_fixed(sqrt(2)) +
    theme_bw()

jmuelmen/sciadv.aaz6433 documentation built on March 29, 2020, 12:58 a.m.