A monolithic FORTRAN solver with 2D and axisymmetric formulation variants suitable for internal flows of solid rocket motor chambers and convergent-divergent nozzles