get_sph_holder Subroutine

    subroutine get_sph_holder(self, region_id, model_holder)
        class(type_material_manager), intent(inout) :: self
        integer(int32), intent(in) :: region_id
        type(holder_sphs), intent(inout) :: model_holder
        integer(int32) :: model_index

        model_index = self%region_id_map(region_id)

        if (model_index == 0) then
            print *, "Error: Invalid region_id in get_sph_holder:", region_id
            call exit(-1)
        end if

        model_holder = self%sph(model_index)

    end subroutine get_sph_holder