def METHOD_NAME(self):
    """Create geometry.
    """
    gmsh.initialize()
    self.box = gmsh.model.occ.add_box(-0.5*self.DOMAIN_X, -0.5*self.DOMAIN_Y, -self.DOMAIN_Z,
                                self.DOMAIN_X, self.DOMAIN_Y, self.DOMAIN_Z)
    self.s_xneg = 1
    self.s_xpos = 2
    self.s_yneg = 3
    self.s_ypos = 4
    self.s_zneg = 5
    self.s_zpos = 6
    p_faulttrace_yneg = gmsh.model.occ.add_point(0.0, -0.5*self.FAULT_LENGTH, 0.0)
    p_faulttrace_ypos = gmsh.model.occ.add_point(0.0, +0.5*self.FAULT_LENGTH, 0.0)
    p_faultbot_yneg = gmsh.model.occ.add_point(0.0, -0.5*self.FAULT_LENGTH, -self.FAULT_WIDTH)
    p_faultbot_ypos = gmsh.model.occ.add_point(0.0, +0.5*self.FAULT_LENGTH, -self.FAULT_WIDTH)
    self.c_faulttrace = gmsh.model.occ.add_line(p_faulttrace_yneg, p_faulttrace_ypos)
    self.c_fault_yneg = gmsh.model.occ.add_line(p_faulttrace_yneg, p_faultbot_yneg)
    self.c_faultbot = gmsh.model.occ.add_line(p_faultbot_yneg, p_faultbot_ypos)
    self.c_fault_ypos = gmsh.model.occ.add_line(p_faulttrace_ypos, p_faultbot_ypos)
    loop_fault = gmsh.model.occ.add_curve_loop([-self.c_faulttrace, self.c_fault_yneg, self.c_faultbot, -self.c_fault_ypos])
    self.s_fault = gmsh.model.occ.add_plane_surface([loop_fault])
    gmsh.model.occ.synchronize()
    gmsh.model.mesh.embed(2, [self.s_fault], 3, self.box)
    gmsh.model.mesh.embed(1, [self.c_faulttrace], 2, self.s_zpos)