do ib=1,10
         ibb=ib/2
         if (2*ibb.lt.ib) ibb=ibb+1
         print *,ib,ibb
       enddo
       stop
       end