add 3 functions

add calculate_s2m(), calculate_curvature_centrifugal(MatrixXd &res), interpolate_circulation()

move some initialization operation from construction function to flow_field_initialization()