Add com.f90 which calculates the center of mass motion