         scalar_surface => extract_surface_field(tke, 'tke_bottom_boundary', "value")         call remap_field(bottom_surface_values, scalar_surface)         scalar_surface => extract_surface_field(tke, 'tke_top_boundary', "value")         call remap_field(top_surface_values, scalar_surface)     end if       