subroutine basm2g(stcprm, xin, xout)
C  receives the vector xin, giving the components of a point in the map
C  system; returns xout, the components of the same point in the geo system.
      dimension stcprm(15),xin(3),xout(3)
      ifind(l,k) = -2 + 3*l + k
        do k=1,3
          xout(k) = 0.
        enddo
        do l=1,3
          do k=1,3
            xout(l) = xout(l) + xin(k) * stcprm(ifind(k,l))
          enddo
        enddo
        return
      end