subroutine xe_xy(stcprm, xi, eta, x, y)
      parameter (REARTH=6371.2)
      dimension stcprm(15)
        x = stcprm(11) + REARTH / stcprm(15) *
     A  (stcprm(13) * xi + stcprm(14) * eta)
        y = stcprm(12) + REARTH / stcprm(15) *
     B  (stcprm(13) * eta - stcprm(14) * xi)
        return
      end