A numerical method for computing initial conditions of Lagrangian invariant tori using the frequency map