program demo_nbppf use M_datapac, only : nbppf implicit none ! call nbppf(x,y) end program demo_nbppf