LCPFCT solver–accessible from Python

shockwave animation of gun fire modeling
LCPFCT fast2d gun shockwave model.

The venerable NRL Flux-Corrected Transport algorithm for Solving Generalized Continuity Equations has been cited in hundreds of publications and is used by myself in my PhD thesis and general research.

I have updated the syntax of the original Fortran code slightly to fit with the FORTRAN 77 standard (not changing algorithm behavior, only syntactical correctness). As a byproduct, I have made demo plotting code available in Python that uses f2py to make the LCPFCT run from Python with 50x speedup overall since you don’t have to write output to disk and then read it back in via a custom parsing algorithm.