def METHOD_NAME(
        lat,
        cellWidthInSouth,
        cellWidthInNorth,
        latTransition,
        latWidthTransition):
    """
    Combine two cell width distributions using a ``tanh`` function. This is
    intended as part of the workflow to make an MPAS global mesh.
    Parameters
    ----------
    lat : ndarray
        vector of length n, with entries between -90 and 90, degrees
    cellWidthInSouth : ndarray
        vector of length n, first distribution
    cellWidthInNorth : ndarray
        vector of length n, second distribution
    latTransition : float
        lat to change from ``cellWidthInSouth`` to ``cellWidthInNorth`` in
        degrees
    latWidthTransition : float
        width of lat transition in degrees
    Returns
    -------
    cellWidthOut : ndarray
        vector of length n, entries are cell width as a function of lat
    """
    cellWidthOut = np.ones(lat.size)
    if latWidthTransition == 0:
        for j in range(lat.size):
            if lat[j] < latTransition:
                cellWidthOut[j] = cellWidthInSouth[j]
            else:
                cellWidthOut[j] = cellWidthInNorth[j]
    else:
        for j in range(lat.size):
            weightNorth = 0.5 * \
                (np.tanh((lat[j] - latTransition) / latWidthTransition) + 1.0)
            weightSouth = 1.0 - weightNorth
            cellWidthOut[j] = weightSouth * cellWidthInSouth[j] + \
                weightNorth * cellWidthInNorth[j]
    return cellWidthOut