def METHOD_NAME(self, params):
    # Renaming for convenience
    zeta = params['zeta']
    zeta_star = params['zeta_star']
    gamma = params['gamma']
    gamma_star = params['gamma_star']
    dist_to_orig = params['dist_to_orig']
    nsurf = len(zeta)
    for isurf in range(nsurf):
        r_surf, ndx1_surf, dx1_surf, dxmax_surf = self.get_all_surface_parameters(isurf)
        M, N = zeta_star[isurf][0, :, :].shape
        for j in range(N):
            zeta_star[isurf][:, 0, j] = zeta[isurf][:, -1, j]
            for i in range(1, M):
                deltax = self.get_deltax(i, dx1_surf, ndx1_surf, r_surf, dxmax_surf)
                zeta_star[isurf][:, i, j] = zeta_star[isurf][:, i - 1, j] + deltax*self.u_inf_direction
        gamma[isurf] *= 0.
        gamma_star[isurf] *= 0.
    for isurf in range(nsurf):
        M, N = zeta_star[isurf][0, :, :].shape
        dist_to_orig[isurf][0] = 0.
        for j in range(0, N):
            for i in range(1, M):
                dist_to_orig[isurf][i, j] = (dist_to_orig[isurf][i - 1, j] +
                                      np.linalg.norm(zeta_star[isurf][:, i, j] -
                                                     zeta_star[isurf][:, i - 1, j]))
            dist_to_orig[isurf][:, j] /= dist_to_orig[isurf][-1, j]