Pore Pressure oscillations and instability

[WaveHello]

There are pore pressures oscillations, due to the weakly compressible nature of the solution (Pore pressures are calculated assuming that the bulk modulus of water isn't as high as the actual value).

[WaveHello]

Added a smoothing procedure (orginally implemented by Abdel). The smoothing happens in the LagrangianPhase.For file. In essence what it does is average the pressures by mapping the pressures to the nodes, and then using the nodes to map the pressures back to the MPs.

    if (CalParams%ApplyPorePressureSmoothing) then
        ! Pressure Smoothing through shape functions - Following Alsardi et al 2024??
        call GetNodalPressure(NodalWaterPressure)
        call MapNodalPressureBackToMPs
    endif

This doesn't completely fix the problem. Honestly, I don't know if it helps. Testing is required for the undrained and the fully saturated conditions.