Write force constants in GPUMD format

We should add functionality to write force constants to GPUMD format.