Testing Constrained Spinup

library('hector')
library('dplyr')

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

library('ggplot2')
library('ggthemes')

Test 1: RCP45 unconstrained

hc <- newcore('../inst/input/hector_rcp45.ini')
run(hc)

Hector core:    unnamed hector core
Start date: 1745
End date:   2300
Current date:   2300
Input file: ../inst/input/hector_rcp45.ini

fetchvars(hc, NA, PREINDUSTRIAL_CO2())

fetchvars(hc, 1745:1750, ATMOSPHERIC_CO2())

spinup_constrain <- fetchvars(hc, 1745:2300, c(ATMOSPHERIC_CO2(), GLOBAL_TEMP()),
                              'spinup-constrain')

core_hector <- read.csv('sample_outputstream_rcp45.csv', 
                        stringsAsFactors = F,
                        comment.char = '#')
core_hector$scenario <- 'core-hector'
core_hector <- 
    filter(core_hector, variable %in% c(ATMOSPHERIC_CO2(), GLOBAL_TEMP()),
           year >= 1745) %>%
    select(scenario, year, variable, value, units)

pltdata <- bind_rows(spinup_constrain, core_hector)
ggplot(data=pltdata, aes(x=year, y=value, color=scenario)) + 
    geom_line(size=1) +
    facet_grid(variable~., scales='free_y') +
    ggtitle('RCP 4.5 with carbon cycle') +
    scale_color_solarized() + 
    theme_solarized_2(light=F)

Test 2: RCP60 with CO2 concentration constrained

hc <- newcore('../inst/input/hector_rcp60_constrained.ini')
run(hc)

Hector core:    unnamed hector core
Start date: 1745
End date:   2300
Current date:   2300
Input file: ../inst/input/hector_rcp60_constrained.ini

fetchvars(hc, NA, PREINDUSTRIAL_CO2())

fetchvars(hc, 1745:1750, ATMOSPHERIC_CO2())

spinup_constrain60 <- fetchvars(hc, 1745:2300, c(ATMOSPHERIC_CO2(), GLOBAL_TEMP()),
                              'spinup-constrain')

core_hector60 <- read.csv('outputstream_rcp60_co2constrain.csv', 
                        stringsAsFactors = F,
                        comment.char = '#')
core_hector60$scenario <- 'core-hector'
core_hector60 <- 
    filter(core_hector60, variable %in% c(ATMOSPHERIC_CO2(), GLOBAL_TEMP()),
           year >= 1745) %>%
    select(scenario, year, variable, value, units)

pltdata60 <- bind_rows(spinup_constrain60, core_hector60)
ggplot(data=pltdata60, aes(x=year, y=value, color=scenario)) + 
    geom_line(size=1) +
    facet_grid(variable~., scales='free_y') +
    ggtitle('RCP 4.5 with constrained CO2 concentration') +
    scale_color_solarized() + 
    theme_solarized_2(light=F)