demo_discre Program

     program demo_discre
     use M_datapac, only : discre
     implicit none
     ! call discre(x,y)
     end program demo_discre