subroutine stvmrc(stcprm, r_lat, r_long)
      real stcprm(15)
C*
C*  Set Map Parameters for a TransVerse MeRCator Projection
C*  Inputs: r_lat,r_longit - latitude and longitude of the reference
C*            point on the meridian tangent to the cylinder,
C*            and 180 degrees from the cut.
C*  Outputs: stcprm - map parameters
C*/
      real refpnt(3),anchrp(3),projpl(3)
      call ll_geo(r_lat, r_long, refpnt)
      call ll_geo(r_lat - 90.,r_long,anchrp)
      anorm = x_prod(anchrp, refpnt, projpl)
      call geo_ll(projpl, p_lat, p_long)
      call mpstrt(stcprm, 0., p_lat,p_long, r_lat, r_long)
      return
      end