demo_subse1 Program

     program demo_subse1
     use M_datapac, only : subse1
     implicit none
     ! call subse1(x,y)
     end program demo_subse1