program demo_discr2 use M_datapac, only : discr2 implicit none ! call discr2(x,y) end program demo_discr2