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