subroutine cpolll(stcprm, alat,along, enx,eny,enz)
      real stcprm(15)
      real map(3),pole(3),geog(3)
      call ll_geo(alat,along, geog)
      call basg2m(stcprm, geog, map)
      do k=1,3
        pole(k) = stcprm(3*k + 1)
      enddo
      call proj_3d(stcprm, map, pole, enx,eny,enz)
      return
      end