  Finalized PHENIX results are presented on $J/\psi$ modification from the RHIC
system-size scan, which include $p$+Al, $p$+Au, and $^{3}$He+Au (small systems)
collisions as well as $AA$ (large systems) collisions at forward and backward
rapidity ($1.2<|y|<2.2$). Comparisons with state-of-the-art theory calculations
and LHC data will be discussed. The forward rapidity data in small systems are
well described by model calculations for all three systems, however, the
backward rapidity $p$+Au and $^{3}$He+Au data are not reproduced by the
shadowing calculation alone. $J/\psi$ $v_2$ meausrements in Au+Au collisions at
forward rapidity is consistent with zero.
