program demo_geoppf use M_datapac, only : geoppf implicit none ! call geoppf(x,y) end program demo_geoppf