[To-Do]: streamline workflow by extracting f1 and f2 from equilibrium and import to softsliding.c

ability to extract facets for f1 and f2 from the equilibrium solution.

Steps:
1. There are two ways this can be done:
- [ ] output_facets() on the fly from inside the equilibrium.c (not recommended as it might be incompatible with the import function in softsliding.c
- [ ] write a python/jupyter script to extract f1, f2 from equilibriumSolution, convert to distance.h compatible format (see: https://github.com/VatsalSy/BurstingBubble_VE_coated/tree/main/initialconditions for examples) -- this method is preferred for customization.

2. Read the new interfaces into softsliding.c -- the code needs to be updated.