demo_sd Program

     program demo_sd
     use M_datapac, only : sd
     implicit none
     ! call sd(x,y)
     end program demo_sd