def flat_to_3d(c_flat, n_species, n):
    """
    Convert flat array of concentrations to 3D array.

    Parameters
    ----------
    c_flat : array_like, shape (n_species*n[0]*n[1])
        Flattened array of concentrations
    n_species :
    """
