We consider a Josephson junction with F 1 F 2 F 3 ferromagnetic trilayers in the ballistic regime, where the magnetization in each ferromagnet F i (i = 1, 2, 3), can have arbitrary orientations and magnetization strengths. The trilayers are sandwiched between two s-wave superconductors with a macroscopic phase difference ∆ϕ. With our generalised theoretical and numerical techniques, we are able to study the planar spatial profiles and ∆ϕdependencies of the charge supercurrents, spin supercurrents, spin torques, and density of states for complex systems that are finite in two dimensions. A broad range of magnetization strengths of the central F 2 layer are considered, from an unpolarized normal metal (N) to a half-metallic phase, supporting only one spin species. Our results reveal that when the magnetization configuration in F 1 F 2 F 3 has three orthogonal components, a supercurrent can flow at ∆ϕ = 0, and a strong second harmonic in the current-phase relation appears. Remarkably, upon increasing the magnetization strength in the central ferromagnet layer up to the half-metallic limit, the self-biased current and second harmonic component become dramatically enhanced, and the critical supercurrent reaches its maximum value. The higher harmonics in the current-phase relations can be controlled by the relative magnetization orientations, with negligible current damping compared to the corresponding F 1 N F 3 counterparts. Additionally, for a broad range of exchange field strengths in the central ferromagnet F 2 , the ground state of the system can be tuned to an arbitrary phase difference ϕ 0 , e.g., by rotating the magnetization in the outer ferromagnet F 3 . For intermediate exchange field strengths in F 2 , a ϕ 0 state can arise that creates a superconducting diode effect, whereby ∆ϕ can be tuned to create a one-way dissipationless current flow. The spatial maps of the spin currents and effective magnetic moments reveal a long-ranged spin torque in the halfmetallic phase. Moreover, the density of states studies demonstrate the emergence of zero energy peaks for the mutually orthogonal magnetization configurations, which is strongest in the half-metallic phase. Our results suggest that this simple trilayer Josephson junction can be an excellent candidate for producing experimentally accessible signatures for long-ranged self-biased supercurrents, and supercurrent diode-effects.