Add and clean some comments of correlation.f90

It is strange that corr_vec1() cannot be run on sigma, while it
is fine to run on iota and float.tk.