def get_rbounds(step):
if step.geom is not None:<EOL><INDENT>rcmb = step.geom.rcmb<EOL><DEDENT>else:<EOL><INDENT>rcmb = step.sdat.par['<STR_LIT>']['<STR_LIT>']<EOL>if step.sdat.par['<STR_LIT>']['<STR_LIT>'].lower() == '<STR_LIT>':<EOL><INDENT>rcmb = <NUM_LIT:0><EOL><DEDENT><DEDENT>rcmb = max(rcmb, <NUM_LIT:0>)<EOL>return rcmb, rcmb + <NUM_LIT:1><EOL>
Radial or vertical position of boundaries.

    Args:
        step (:class:`~stagpy.stagyydata._Step`): a step of a StagyyData
            instance.
    Returns:
        tuple of floats: radial or vertical positions of boundaries of the
        domain.