This project numerically modeled a simplified version of a fault system complete with a shear zone, similar to that of the fault zones in the Sierra Nevada, purely in Python to compare to that of one using ABAQUS. Modeling fluid flow in this system can help explain or predict where water is coming to the surface in places like the Sierra Nevada. The model took into account fluid flow based on Darcy’s law, which explains how fluid flows through a porous media, and applied it to the background diffusivities and pressures that it was given. This gave fluid velocities off initial pressure gradients and perturbations. After testing several initial pressure fields in the model against the previous models, it was found that, with Test 3, it was possible to match the “black box” of ABAQUS using the freeware of Python. This result opens up opportunities to explore further this system by refining the model and adding new parameters that were just not possible to do with ABAQUS.