sPuReMD: fix code to remove several compiler warnings. Remove unused code. Disable compilation of several debug routines. Other code clean-up.

sPuReMD: add support for user-defined molecular charge constraints with EEM (using MOLCHARGE keyword in BGF geometry files). Fix issue with far neighbor list not being initialized for certain edge cases.

sPuReMD: enable setting the number of OpenMP threads via control file keyword or OpenMP environment variable (OMP_NUM_THREADS).

sPuReMD: make zlib dependency optional (trajectory file compression). PG-PuReMD: remove erroneous zlib header file inclusion.

sPuReMD: fix OpenMP compilation issue with specifications prior to 4.0 (behavior change regarding implicit shared attribute assignment to const variables).

QMMM interface related: ignore QM-MM bonds while forming bond lists, update charge matrix so that the total charge constraint only applies to QM atoms

sPuReMD: standardize QM/MM interface. Add test codes for QM/MM interface for AVE dataset in Python and Fortran.

sPuReMD: merge changes from branch 'new_setup'. Rework interface for geometry initialization without file I/O. Other code clean-up.

sPuReMD: fix BGF parser to correctly remap atom coordinates to internal simulation box boundaries. Fix reax_list to handle case when total interactions is zero instead of terminating program.

sPuReMD: rework Fortran F90 interfaces for Amber specifics and declare this interface for such usage.

sPuReMD: fix Coulomb energy calculations for QM/MM mode. Indicate QEq model not supported with QM/MM mode. Add Fortran interface for QM/MM mode.

sPuReMD: add QMMM library interface for external codes (e.g., Amber). Rework file I/O parsing to jive with multiple consecutive simulation code changes. Improve file parser error checking (numeric types). Other code clean-up.

sPuReMD: fixes for multiple consecutive simulation codepath. Reorganization to this codepath (some data structure now carry a flag to designation their allocation state). Fix bug in force field parser (atom type not truncated correctly). Charge sparse_matrix structs to be non-pointers to jive with new allocation state flag. Other code clean-up.

sPuReMD: finalize sPuReMD multiple consecutive simulation code (for integration with other codebases).

PG-PuReMD: fixes to NVT with Berendsen thermostat. Use updated constants from sPuReMD. Add BGF geometry file parser. Other code clean-up.

sPuReMD: fix memory leak with ICHOLT preconditioner (droptol for OpenMP only) for charge solver. Default to Jacobi preconditioner. Resolve several compiler warnings around C strings (snprintf, etc.). Other code clean-up.

sPuReMD: add functionality to run multiple MD simulations (for usage via integration with external codebases).

sPuReMD: resolve OpenMP build errors with GCC v10 (behavior changes for function arguments). Suppress possibly unitialized variable warnings for BiCGStab.

sPuReMD: finalize corrections for pressure calculations. Change output units from GPa to ATMs. Other formatting changes.