def zonal_mean_column(num_lat=<NUM_LIT>, num_lev=<NUM_LIT:30>, water_depth=<NUM_LIT>, lat=None,<EOL>lev=None, **kwargs):
if lat is None:<EOL><INDENT>latax = Axis(axis_type='<STR_LIT>', num_points=num_lat)<EOL><DEDENT>elif isinstance(lat, Axis):<EOL><INDENT>latax = lat<EOL><DEDENT>else:<EOL><INDENT>try:<EOL><INDENT>latax = Axis(axis_type='<STR_LIT>', points=lat)<EOL><DEDENT>except:<EOL><INDENT>raise ValueError('<STR_LIT>')<EOL><DEDENT><DEDENT>if lev is None:<EOL><INDENT>levax = Axis(axis_type='<STR_LIT>', num_points=num_lev)<EOL><DEDENT>elif isinstance(lev, Axis):<EOL><INDENT>levax = lev<EOL><DEDENT>else:<EOL><INDENT>try:<EOL><INDENT>levax = Axis(axis_type='<STR_LIT>', points=lev)<EOL><DEDENT>except:<EOL><INDENT>raise ValueError('<STR_LIT>')<EOL><DEDENT><DEDENT>depthax = Axis(axis_type='<STR_LIT>', bounds=[water_depth, <NUM_LIT:0.>])<EOL>slab = SlabOcean(axes={'<STR_LIT>':latax, '<STR_LIT>':depthax}, **kwargs)<EOL>atm = Atmosphere(axes={'<STR_LIT>':latax, '<STR_LIT>':levax}, **kwargs)<EOL>return slab, atm<EOL>
Creates two Domains with one water cell, a latitude axis and
    a level/height axis.

        * SlabOcean:    one water cell and a latitude axis above
          (similar to :func:`zonal_mean_surface`)
        * Atmosphere: a latitude axis and a level/height axis (two dimensional)


    **Function-call argument** \n

    :param int num_lat:         number of latitude points on the axis
                                [default: 90]
    :param int num_lev:         number of pressure levels
                                (evenly spaced from surface to TOA) [default: 30]
    :param float water_depth:   depth of the water cell (slab ocean) [default: 10.]
    :param lat:                 specification for latitude axis (optional)
    :type lat:                  :class:`~climlab.domain.axis.Axis` or latitude array
    :param lev:                 specification for height axis (optional)
    :type lev:                  :class:`~climlab.domain.axis.Axis` or pressure array
    :raises: :exc:`ValueError`  if `lat` is given but neither Axis nor latitude array.
    :raises: :exc:`ValueError`  if `lev` is given but neither Axis nor pressure array.
    :returns:                   a list of 2 Domain objects (slab ocean, atmosphere)
    :rtype:                     :py:class:`list` of :class:`SlabOcean`, :class:`Atmosphere`

    :Example:

        ::

            >>> from climlab import domain
            >>> sfc, atm = domain.zonal_mean_column(num_lat=36,num_lev=10)

            >>> print sfc
            climlab Domain object with domain_type=ocean and shape=(36, 1)

            >>> print atm
            climlab Domain object with domain_type=atm and shape=(36, 10)