In this work, we derived the L-MFE method to incorporate Lindblad jump operator dynamics into the meanfield Ehrenfest (MFE) approach. We took the density matrix equations of motion for Lindblad dynamics and mapped them onto an ensemble of pure state coefficients, using trajectory averages and expectation values of random variables. We then derived the L-MFE method to update the MFE coefficients at each time step which rigorously satisfies Lindblad jump operator dynamics. This established a method that exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase (on coefficients c 0,ξ (t) in Eq. 69b).