
   #@ [2d-serial]
   #@ exe = mechanics
   #@ dim = 2
   #@ nprocs = 1
   #@ check = true
   #@ timeout = 120
   #@
   #@ [2d-parallel]
   #@ exe = mechanics
   #@ dim = 2
   #@ nprocs = 2
   #@ check = true
   #@ timeout = 120
   #@
   
   alamo.program = mechanics
   alamo.program.mechanics.model = finite.neohookean
   
   type = static
   timestep = 0.1
   stop_time = 0.1
   
   amr.plot_int = 1
   amr.max_level = 1
   amr.n_cell = 32 32 8
   amr.blocking_factor = 4
   amr.max_grid_size = 16
   amr.regrid_int = -1
   amr.grid_eff = 1.0
   amr.node.all = 1
   amr.cell.all = 1
   
   explicitmesh.on = 1
   explicitmesh.lo1 = 8 8 0
   explicitmesh.hi1 = 23 23 0
   
   geometry.prob_lo = 0 0 0
   geometry.prob_hi = 1 1 1
   geometry.is_periodic = 0 0 0
   
   ic.type = ellipse
   ic.ellipse.a = 0.22 0.22 0.22
   ic.ellipse.x0 = 0.5 0.5 0.5
   ic.ellipse.eps = 0.04
   
   nmodels = 2
   model1.mu = 1
   model1.kappa = 1
   model2.mu = 26000
   model2.kappa = 70000
   
   rhs.type = expression
   rhs.expression.region0 = "-(x-0.5)/0.04*exp(-((sqrt((x-0.5)^2+(y-0.5)^2)-0.22)/0.04)^2)"
   rhs.expression.region1 = "-(y-0.5)/0.04*exp(-((sqrt((x-0.5)^2+(y-0.5)^2)-0.22)/0.04)^2)"
   
   solver.conservative_face_flux = 1
   solver.line_search = 1
   solver.resync_coeffs = 1
   solver.nriters = 3
   solver.nrtolerance = 0
   solver.max_iter = 1500
   solver.pre_smooth = 2
   solver.post_smooth = 2
   solver.bottom_solver = bicgstab
   solver.verbose = 3
   
   max_coarsening_level = 2
   zero_out_displacement = 1
   tol_rel = 2.0e-5
   tol_abs = 1.0e-8
   ref_threshold = 100
   
   bc.type = constant
   bc.constant.type.xlo = disp disp
   bc.constant.type.xhi = disp trac
   bc.constant.type.ylo = trac disp
   bc.constant.type.yhi = disp disp
   bc.constant.type.xloylo = disp disp
   bc.constant.type.xloyhi = disp disp
   bc.constant.type.xhiylo = disp disp
   bc.constant.type.xhiyhi = disp disp