subroutine cgrnxy(stcprm, x,y, enx,eny,enz)
      real stcprm(15)
      real map(3),pole(3)
      call xy_map(stcprm, x,y, map)
      do k=1,3
        pole(k) = stcprm(3*k - 1)
      enddo
      call proj_3d(stcprm, map, pole, enx,eny,enz)
      return
      end