def METHOD_NAME():
    """Test calculation of column-average."""
    cube = get_cube([[[[1.0]], [[2.0]]]], air_pressure_coord=True,
                    var_name='ch4', units='1')
    hus_cube = get_cube([[[[0.2]], [[0.2]]]], air_pressure_coord=True,
                        var_name='hus', units='1')
    zg_cube = get_cube([[[100.0]]], air_pressure_coord=False, var_name='zg',
                       units='m')
    ps_cube = get_cube([[[100000.0]]], air_pressure_coord=False, var_name='ps',
                       units='Pa')
    x_cube = column_average(cube, hus_cube, zg_cube, ps_cube)
    assert x_cube.shape == (1, 1, 1)
    assert x_cube.units == '1'
    assert x_cube.coords('time')
    assert x_cube.coords('air_pressure')
    assert x_cube.coords('latitude')
    assert x_cube.coords('longitude')
    np.testing.assert_allclose(x_cube.data, [[[1.85]]])
    np.testing.assert_allclose(x_cube.coord('time').points, [0.0])
    np.testing.assert_allclose(x_cube.coord('air_pressure').points, 85000.0)
    np.testing.assert_allclose(x_cube.coord('air_pressure').bounds,
                               [[80000.0, 90000.0]])
    np.testing.assert_allclose(x_cube.coord('latitude').points, [45.0])
    np.testing.assert_allclose(x_cube.coord('longitude').points, [10.0])