demo_geocdf Program

     program demo_geocdf
     use M_datapac, only : geocdf
     implicit none
     ! call geocdf(x,y)
     end program demo_geocdf