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):
        M, N = zeta_star[isurf][0, :, :].shape
        angle = 0.
        for i in range(M):
            # Compute the step in azimuthal angle
            angle -= self.get_dphi(i, self.dphi1, self.ndphi1, self.r, self.dphimax)
            delta_t = -angle/np.linalg.norm(self.rotation_velocity)
            rot = algebra.rotation_matrix_around_axis(algebra.unit_vector(self.rotation_velocity), angle)
            for j in range(N):
                # Define the helicoidal
                aux_zeta_TE = zeta[isurf][:, -1, j] - (self.h_ref - self.h_corr)*self.shear_direction
                aux_zeta_TE = np.dot(rot, aux_zeta_TE) + (self.h_ref - self.h_corr)*self.shear_direction
                # Translate according to u_inf depending on the height
                h = np.dot(aux_zeta_TE, self.shear_direction) + self.h_corr
                zeta_star[isurf][:, i, j] = aux_zeta_TE + self.u_inf*self.u_inf_direction*(h/self.h_ref)**self.shear_exp*delta_t
                # zeta_star[isurf][:, i, j] = zeta[isurf][:, -1, j] + self.u_inf*self.u_inf_direction*self.dt*i
                # print(zeta_star[isurf][:, i, j])
        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]