Why magrobs_[iObs][iSurf][i] = mag(robs_[iObs][iSurf][i]) is not used, instead caculated the |r| = (-U_0*x/c_0+R)/(1-(U_0/c_0)^2).