subroutine xy_xe(stcprm,  x,y, xi,eta)
      real stcprm(15)
      parameter (REARTH=6371.2)
        xtmp = x - stcprm(11)
        ytmp = y - stcprm(12)
        xi = stcprm(15) / REARTH *
     C    (stcprm(13) * xtmp - stcprm(14) * ytmp)
        eta = stcprm(15) / REARTH *
     C    (stcprm(13) * ytmp + stcprm(14) * xtmp)
        return
      end