This document is used to produce some tables from the Kampala test model output. The tables can be compared with other model revisions.
Print the attribution of sanitation types to the surface water emissions.
fpath <- testthat::test_path("output","kla","human_sources_rotavirus_kla_water.csv")
n_digits <- 5
df_human_sources_water <- read.csv(fpath)
df_print <- df_human_sources_water %>%
dplyr::mutate_if(is.numeric, format, digits = n_digits, scientific = TRUE) %>%
dplyr::rename_with(~ stringr::str_trunc(.x, 10))
df_print %>% dplyr::select(iso, 2:8)
| iso | bucketL… | compost… | contain… | flushOpen | flushPit | flushSe… | flushSewer |
|---|---|---|---|---|---|---|---|
| 1 | 0 | 1.0824e+10 | 0 | 5.8597e+12 | 6.7309e+10 | 2.3323e+11 | 1.1893e+13 |
| 2 | 0 | 3.7942e+09 | 0 | 4.2429e+12 | 4.5091e+11 | 1.1137e+12 | 0.0000e+00 |
| 3 | 0 | 1.8002e+10 | 0 | 1.1170e+13 | 7.2373e+11 | 2.2734e+12 | 0.0000e+00 |
| 4 | 0 | 6.2749e+09 | 0 | 6.5522e+12 | 5.0427e+11 | 2.1585e+12 | 0.0000e+00 |
| 5 | 0 | 7.4791e+09 | 0 | 5.2296e+12 | 7.8373e+11 | 1.4512e+12 | 1.8953e+13 |
df_print %>% dplyr::select(iso, 9:14)
| iso | flushUn… | hanging… | openDef… | other | pitNoSlab | pitSlab |
|---|---|---|---|---|---|---|
| 1 | 3.7176e+12 | 0.0000e+00 | 5.5511e+12 | 1.4867e+12 | 1.5559e+11 | 2.4404e+12 |
| 2 | 1.3258e+13 | 2.6803e+11 | 1.1617e+13 | 1.8836e+12 | 6.2408e+11 | 2.8806e+13 |
| 3 | 5.1082e+13 | 1.0084e+11 | 1.7778e+13 | 5.3607e+12 | 1.1778e+12 | 2.5170e+13 |
| 4 | 3.8535e+13 | 1.6317e+12 | 1.6672e+13 | 1.9984e+12 | 1.1613e+12 | 1.9031e+13 |
| 5 | 3.1370e+13 | 5.8429e+10 | 1.0857e+13 | 1.2768e+13 | 1.5000e+12 | 3.0162e+13 |
Print the attribution from various pathways to the surface water emissions.
fpath <- testthat::test_path("output","kla","human_emissions_rotavirus_kla_surface_water.csv")
df_pathways_water <- read.csv(fpath)
df_pathways_water %>%
dplyr::mutate_if(is.numeric, format, digits = n_digits, scientific = TRUE)
| iso | humans | land | wwtp |
|---|---|---|---|
| 1 | 1.9207e+13 | 7.1999e+09 | 1.2201e+13 |
| 2 | 6.2183e+13 | 8.5869e+10 | 0.0000e+00 |
| 3 | 1.1477e+14 | 8.1338e+10 | 0.0000e+00 |
| 4 | 8.8188e+13 | 6.3329e+10 | 0.0000e+00 |
| 5 | 9.3604e+13 | 9.2557e+10 | 1.9444e+13 |
fpath <- testthat::test_path("output","kla","pathways.tif")
rast_pathways <- terra::trim(terra::rast(fpath), padding = 0)
terra::plot(log10(rast_pathways), col = hcl.colors(50, palette = "Geyser"), buffer = TRUE)