FORTRAN implementation of the Finite Difference Frequency Domain method