demo_subset Program

     program demo_subset
     use M_datapac, only : subset
     implicit none
     ! call subset(x,y)
     end program demo_subset