def fixedBoundary(eq, Jtor, psi):
    """
    Set psi=0 on all boundaries
    
    Inputs
    ------
    
    eq    Equilibrium object (not used)
    Jtor  2D array of toroidal current (not used)
    psi   2D array of psi values (modified by call)
    
    Returns
    -------
    None
    """
    psi[(0), :] = 0.0
    psi[:, (0)] = 0.0
    psi[(-1), :] = 0.0
    psi[:, (-1)] = 0.0
