{
    "_comments": {
        "##Geometry Parameters##": {
            "xmin": "Minimum X coordinate for the domain (m)",
            "xmax": "Maximum X coordinate for the domain (m)",
            "ymin": "Minimum Y coordinate for the domain (m)",
            "ymax": "Maximum Y coordinate for the domain (m)",
            "zmin": "Minimum Z coordinate for the domain (m)",
            "zmax": "Maximum Z coordinate for the domain (m)",
            "lc": "Mesh size at far boundary (m)",
            "lc_fault": "Mesh size at near fault (m)"
        },
        "##Material Parameters##": {
            "lambda_o": "Lamé Constant λ (Pa)",
            "shear_modulus_o": "Shear Modulus μ (Pa)",
            "rho": "Density ρ (kg/m^3)"
        },
        "##Continuum Damage-Breakage Model Parameters##": {
            "xi_0": "Strain invariant ratio at onset of damage ξ_0",
            "xi_d": "Strain invariant ratio at onset of breakage ξ_d",
            "Cd_constant": "Damage accumulation rate C_d (1/s)",
            "CdCb_multiplier": "Breakage accumulation rate multiplier Cm (C_B (1/s) = Cm * C_d)",
            "CBH_constant": "Breakage healing rate C_BH (1/s)",
            "C_1": "Damage healing rate C_1 (1/s)",
            "C_2": "Damage healing rate C_2 (1/s)",
            "beta_width": "Width of transitional region β",
            "C_g": "Compliance of fluidity of the fine grain material C_g (1/(Pa s))",
            "m1": "Coefficient of power law index m_1",
            "m2": "Coefficient of power law index m_2",
            "chi": "Ratio of two energy states χ"
        },
        "##Initial Damage Parameters##": {
            "nucl_center_x": "Nucleation center x coordinate (m)",
            "nucl_center_y": "Nucleation center y coordinate (m)",
            "nucl_center_z": "Nucleation center z coordinate (m)",
            "len_fault_strike": "Length of initial damage strip along strike direction (m)",
            "len_fault_dip": "Length of initial damage strip along dip direction (m)",
            "len_fault_normal": "Length of initial damage strip along normal direction (m)",
            "nucl_distance": "Length of nucleation patch along strike/dip direction (m)",
            "nucl_thickness": "Length of nucleation patch along normal direction (m)",
            "nucl_damage": "Value of initial damage strip",
            "e_damage": "Value of additional nucleation damage",
            "e_sigma": "Value of initial damage exponential decay",
            "duration": "Simulation time (s)"
        },
        "##Boundary Condition##": {
            "normal_traction_x": "Normal boundary traction along x direction (Pa)",
            "normal_traction_z": "Normal boundary traction along z direction (Pa)",
            "normal_traction_y": "Normal boundary traction along y direction (Pa)",
            "shear_traction": "Shear boundary traction along x-z direction (Pa)"
        },
        "##Simulation Parameters##": {
            "dt": "Simulation time step (s)",
            "end_time": "Total simulation time (s)",
            "time_step_interval": "Result saved time interval"
        }
    },
    "xmin": -10000,
    "xmax": 10000,
    "ymin": -10000,
    "ymax": 10000,
    "zmin": -10000,
    "zmax": 10000,
    "lc": 5000,
    "lc_fault": 100,
    "lambda_o": 3e10,
    "shear_modulus_o": 3e10,
    "rho": 2700,
    "xi_0": -0.8,
    "xi_d": -0.9,
    "Cd_constant": 1e4,
    "CdCb_multiplier": 1000,
    "CBH_constant": 1e4,
    "C_1": 300,
    "C_2": 0.05,
    "beta_width": 0.03,
    "C_g": 1e-10,
    "m1": 10,
    "m2": 1,
    "chi": 0.8,
    "nucl_center_x": 0,
    "nucl_center_y": -5000,
    "nucl_center_z": 0,
    "len_fault_strike": 5000,
    "len_fault_dip": 5000,
    "len_fault_normal": 1000,
    "nucl_distance": 400,
    "nucl_thickness": 200,
    "nucl_damage": 0.7,
    "e_damage": 0.3,
    "e_sigma": 250,
    "duration": 0.1,
    "normal_traction_x": 135e6,
    "normal_traction_z": 120e6,
    "normal_traction_y": 127.5e6,
    "shear_traction": 55e6,
    "dt": 1e-4,
    "end_time": 1e-3,
    "time_step_interval": 1
}
