https://devhubby.com/thread/what-is-the-correct-way-to-do-subtraction-in-fortran