user parameters [
     # NAMES
     simulation_name      = LVFZ_dyn_N900_Ht1200.0_fctx1.5_fctz1.5_Hh800.0_sigi0_InfB1
     load_simulation_name = LVFZ_sts_N900_Ht1200.0_fctx1.5_fctz1.5_Hh800.0_sigi0
     # FOLDERS
     output_folder       = /home/ga288/raw-output/
     paraview_folder     = /home/ga288/raw-output/
     restart_load_folder = /home/ga288/input-data/restart/static_equilibrium/
     restart_dump_folder = /home/ga288/input-data/restart/
     # BOUNDARY CONDITIONS
     top_traction   = [ 2.75e+07 , -0.00e+00 , 0. ]
     bot_traction   = [-2.75e+07 ,  0.00e+00 , 0. ]
     left_traction  = [     0 , -2.75e+07 , 0. ]
     right_traction = [     0 ,  2.75e+07 , 0.]
     #block_x_bcs = slider_top,slider_left,slider_right,slider_back,slider_front
     #block_y_bcs = slider_top,slider_left,slider_right,slider_back,slider_front
     #block_z_bcs = slider_top,slider_left,slider_right,slider_back,slider_front
     
     #block_z_bcs  = slider_back

     partition_method = tagging
     nb_procs = [48,1,1]

     # INTERFACE
     vertical_normals = true
     
     # NUCLEATION
     is_TPV3 = true # instanteneous nucleation square
     nuc_x_min    = 43490.0
     nuc_x_max    = 46510.0
     nuc_z_min    = 20990.0
     nuc_z_max    = 24010.0
     tau_nuc      = 3.100e+07
     tau_0        = 2.75e+07

     # STOP
     stop_x_min = 14990.000
     stop_x_max = 75010.000
     stop_z_min = 7490.000
     stop_z_max = 37510.000

     # OFF FAULT HET 
     off_fault_het = 1
     off_fault_het_x_start_position = 0.000e+00
     off_fault_het_x_trans_position = 9.000e+04
     off_fault_het_y_start_position = 0.000e+00
     off_fault_het_y_trans_position = 8.000e+02
     off_fault_het_z_start_position = 0.000e+00
     off_fault_het_z_trans_position = 4.500e+04

     # MESH
     spatial_dimension = 3
     mesh              = /home/ga288/input-data/meshes/block_3d_L90000.0_Ht1200.0_Wt45000.0_msh900_hexa8.msh
     antisym_setup     = true
     
     # TIME
     simulation_time  = 10
     time_step_factor = 0.10

     # INFINITE BOUNDARY CONDITION
     inf_bc = 1
     inf_bc_neumann = 1
     inf_top_nb_nodes_x = 900
     inf_top_length_x   = 90000.0
     inf_top_nb_nodes_z = 450
     inf_top_length_z   = 45000.0
     inf_top_rho  = 2670.0
     inf_top_nu   = 0.25
     inf_top_E    = 8.01e+10
     
     inf_bot_nb_nodes_x = 900
     inf_bot_length_x   = 90000.0
     inf_bot_nb_nodes_z = 450
     inf_bot_length_z   = 45000.0
     inf_bot_rho  = 2670.0
     inf_bot_nu   = 0.25
     inf_bot_E    = 8.01e+10
     

     # DUMP
     #dump_heights          = [1200.0]
     #dump_height_tolerance = 5e-8  

     # NODE DUMP
     #node_dump_heights     = [-1200.0,1200.0]
     #node_dump_height_tolerance = 5e-8


     # DUMP AT A PRECISE TIME
     # precise_dump_time = 6.51586e-05

     # FULL DUMP
     # full_dump_positions = [ 0.1, 0.12, 0.14, 0.16, 0.18 ]
     
     # FREQUENCY THRESHOLD
     rel_slip_for_highf = 0.0 # 0.1
          
     # INTERFACE TEXT DUMPER 
     interface_dump_highf = 1.00e-02
     interface_dump_lowf  = 1.00e-02#3e-7
     interface_dump_fields = cohesive_shear,shear_gap_rate,displacement
     
     # AT DISTANCE TEXT DUMPER
     at_distance_dump_highf = 1.00e-02
     at_distance_dump_lowf  = 1.00e-02#3e-7

     # GLOBAL TEXT DUMPER
     global_dump_highf = 1.00e-02
     global_dump_lowf  = 1.00e-02#1e-7
     
     # FULL PARAVIEW DUMPER  # value cannot be given in 1e6 format
     nb_max_paraview_dumps = 0
     int_paraview_dumps    = 10000000
]


friction linear_shear_cohesive [
     G_c   = 5.00e+05
     tau_c = 3.00e+07
     tau_r = 2.50e+07
]



material elastic [
     name = slider

     rho  = 2670.0
     nu   = 0.25
     E    = 8.01e+10
]
material elastic [
     name = heterogeneity
     rho  = 2670.0
     nu   = 0.25
     E    = 5.13e+10
]
material elastic [
     name = base
     rho  = 2670.0
     nu   = 0.25
     E    = 8.01e+10
]