program demo_hfnran use M_datapac, only : hfnran implicit none ! call hfnran(x,y) end program demo_hfnran