demo_discr3 Program

     program demo_discr3
     use M_datapac, only : discr3
     implicit none
     ! call discr3(x,y)
     end program demo_discr3