def podaac_to_harmony_bbox(bbox):
    """
    Convert PO.DAAC bbox ((west, east), (south, north))
    to Harmony bbox (west, south, east, north)

    Parameters
    ----------
    bbox : np.array
        Podaac bbox

    Returns
    -------
    array, int or float
        Harmony bbox
    """
    return [bbox[0][0], bbox[1][0], bbox[0][1], bbox[1][1]]