def sdss2decam(g_sdss, r_sdss, i_sdss, z_sdss):
    """
    Converts SDSS magnitudes to DECam magnitudes

    Args:
        [griz]_sdss: SDSS magnitudes (float or arrays of floats)

    Returns:
        g_decam, r_decam, z_decam

    Note: SDSS griz are inputs, but only grz (no i) are output
    """
    gr = g_sdss - r_sdss
    ri = r_sdss - i_sdss
    iz = i_sdss - z_sdss

    # - DESI-1788v1 equations 4-6
    g_decals = g_sdss + 0.01684 - 0.11169*gr
    r_decals = r_sdss - 0.03587 - 0.14114*ri
    z_decals = z_sdss - 0.00756 - 0.07692*iz

    return g_decals, r_decals, z_decals