We numerically demonstrate heat rectification for linear chains of ions in trap lattices with graded trapping frequencies, in contact with thermal baths implemented by optical molasses. To calculate the local temperatures and heat currents we find the stationary state by solving a system of algebraic equations. This approach is much faster than the usual method that integrates the dynamical equations of the system and averages over noise realizations.