run_cases_for_delivery_distributed.jl

using Distributed
addprocs(8)
##
@everywhere using Jutul, JutulDarcy, HYPRE, CSP11
cases_b = [
    # "nudge_b_400x60",
    # "nudge_b_819x117",
    "nudge_b_2640x380"
]

cases_c = [
    # "nudge_c_10x10x10",
    # "nudge_c_85x50x40",
    "nudge_c_170x100x50"
]

allcases = [cases_b, cases_c]
# allcases = [["b_400x60"]]
# allcases = [["c_10x10x10"]]

kgrad_to_run = [:tpfa, :avgmpfa]


all_cases = []
for cases_to_run in allcases
    for kg in kgrad_to_run
        for basename in cases_to_run
            push!(all_cases, (basename, kg))
        end
    end
end

@distributed for mycase in all_cases
    basename, kg = mycase
    specase = Symbol(basename[7])
    case, name = setup_spe11_case_from_mrst_grid(basename,
        case = specase,
        thermal = true,
        use_reporting_steps = true,
        kgrad = kg
    );
    domain = reservoir_domain(case.model)
    nc = number_of_cells(domain)
    println("Running case $name with $nc cells")
    pth = jutul_output_path(name, subfolder = "csp11_delivery_aug30")
    hook, csv_pth, f = CSP11.get_reporting_hook(pth, domain, specase = specase)
    simulate_reservoir(case,
        output_path = pth,
        max_timestep = 1*si_unit(:year),
        report_level = 0,
        tol_cnv = 0.01,
        nonlinear_tolerance = Inf, # Thermal, for AvgMPFA
        post_ministep_hook = hook,
        relaxation = SimpleRelaxation(),
        max_nonlinear_iterations = 20,
        max_timestep_cuts = 100,
        output_states = false,
        output_reports = false,
        info_level = 1
    );
    # Return nothing, results will have been written to disk
    nothing
end