case.py

#!/usr/bin/env python3

import json

# Configuring case dictionary
print(json.dumps({
    # Logistics ================================================
    'case_dir'                     : '\'.\'',
    'run_time_info'                : 'T',
    # ==========================================================

    # Computational Domain Parameters ==========================
    'x_domain%beg'                 : 0.E+00,
    'x_domain%end'                 : 1.E+00,
    'y_domain%beg'                 : 0.E+00,
    'y_domain%end'                 : 1.E+00,
    'm'                            : 99,
    'n'                            : 99,
    'dt'                           : 5.E-07,
    't_step_start'                 : 0,
    't_step_stop'                  : 1000,
    't_step_save'                  : 100,
    # ==========================================================

    # Simulation Algorithm Parameters ==========================
    'num_patches'                  : 2,
    'num_fluids'                   : 2,
    'time_stepper'                 : 3,
    'weno_order'                   : 5,
    'weno_eps'                     : 1.E-16,
    'bc_x%beg'                     : -3,
    'bc_x%end'                     : -3,
    'bc_y%beg'                     : -3,
    'bc_y%end'                     : -3,
    # ==========================================================

    # Formatted Database Files Structure Parameters ============
    'format'                       : 1,
    'precision'                    : 2,
    'prim_vars_wrt'                :'T',
    'parallel_io'                  :'T',
    # ==========================================================
                                                            
    # Patch 1: Base ============================================
    'patch_icpp(1)%geometry'       : 3,
    'patch_icpp(1)%x_centroid'     : 0.5E+00,
    'patch_icpp(1)%y_centroid'     : 0.5E+00,
    'patch_icpp(1)%length_x'       : 1.E+00,
    'patch_icpp(1)%length_y'       : 1.E+00,
    'patch_icpp(1)%vel(1)'         : 100.E+00,
    'patch_icpp(1)%vel(2)'         : 100.E+00,
    'patch_icpp(1)%pres'           : 1.E+05,
    'patch_icpp(1)%alpha_rho(1)'   : 1000.E+00,
    'patch_icpp(1)%alpha_rho(2)'   : 1.,
    'patch_icpp(1)%alpha(1)'       : 1.E-12,
    'patch_icpp(1)%alpha(2)'       : 1. - 1.E-12,
    # ==========================================================

    # Patch 2: Density to transport ============================
    'patch_icpp(2)%geometry'       : 2,
    'patch_icpp(2)%smoothen'       : 'T',
    'patch_icpp(2)%smooth_patch_id' : 1,
    'patch_icpp(2)%smooth_coeff'   : 0.5E+00,
    'patch_icpp(2)%x_centroid'     : 0.1E+00,
    'patch_icpp(2)%y_centroid'     : 0.1E+00,
    'patch_icpp(2)%radius'         : 0.1E+00,
    'patch_icpp(2)%alter_patch(1)' : 'T',
    'patch_icpp(2)%vel(1)'         : 100.E+00,
    'patch_icpp(2)%vel(2)'         : 100.E+00,
    'patch_icpp(2)%pres'           : 1.E+05,
    'patch_icpp(2)%alpha_rho(1)'   : 1.,
    'patch_icpp(2)%alpha_rho(2)'   : 1.0,
    'patch_icpp(2)%alpha(1)'       : 0,
    'patch_icpp(2)%alpha(2)'       : 1.,
    # ==========================================================

    # Fluids Physical Parameters ===============================
    'fluid_pp(1)%gamma'            : 1.E+00/(2.35E+00-1.E+00),
    'fluid_pp(1)%pi_inf'           : 2.35E+00*1.E+09/(2.35E+00-1.E+00),
    'fluid_pp(2)%gamma'            : 1.E+00/(1.4E+00-1.E+00),
    'fluid_pp(2)%pi_inf'           : 0.E+00,
    # ==========================================================
}))
# ==============================================================================