subroutine ccrvll(stcprm, alat, along, gx, gy)
      parameter (REARTH = 6371.2)
      real stcprm(15)
      real map(3),geog(3)
      call ll_geo(alat,along, geog)
      call basg2m(stcprm, geog, map)
      if ( abs(map(3)) .ge. 1.) then
        if (stcprm(1) .eq. map(3)) then
          gx = 0.
          gy = 0.
          return
        else
          xi = 0.
          eta = sign(.5e21,map(3)) / REARTH
        endif
      else
        fact = (stcprm(1) - map(3)) /
     C    (1. - map(3)) / (1. + map(3)) / REARTH
        xi = - map(2) * fact
        eta = map(1) * fact
        if (abs(stcprm(1)) .lt. 1.) then
          glambda = (stcprm(1) - 1.) * atan2(map(2), map(1))
          clambda = cos(glambda)
          slambda = sin(glambda)
          fact = xi * clambda - eta * slambda
          eta = xi * slambda + eta * clambda
          xi = fact
        endif
      endif
      gx = xi * stcprm(13) + eta * stcprm(14)
      gy = eta * stcprm(13) - xi * stcprm(14)
      return
      end