subroutine finish_theta_grid_eik use geometry, only: finish_geometry initialized = .false. call finish_geometry call theta_grid_eik_config%reset() end subroutine finish_theta_grid_eik