Implement `n_max_surf`, which is the maximum toroidal mode number that can be resolved by the flux surface. A variable exists in XGC called `m_max_surf`. `n_max_surf` can be obtained from `m_max_surf` by dividing by the safety factor `q`.