SUBROUTINE TMSFMC(PHI,DLAM,SFO,A,B,SF,C)                          TMSFMC01
C***********************************************************************TMSFMC02
C*                                                                      TMSFMC03
C*       THIS ROUTINE COMPUTES THE POINT SCALE FACTOR AND MERIDIAN      TMSFMC04
C*     CONVERGENCE (FOR A POINT DEFINED BY PHI,DLAM) FOR A TRANSVERSE   TMSFMC05
C*     MERCATOR PROJECTION DEFINED BY THE SCALE FACTOR SFO AT THE       TMSFMC06
C*     CENTRAL MERIDIAN.                                                TMSFMC07
C*                                                                      TMSFMC08
C*                                                                      TMSFMC09
C*       INPUT:                                                         TMSFMC10
C*             PHI - ELLIPSOIDAL LATITUDE OF THE POINT, IN RADIANS.     TMSFMC11
C*             DLAM - ELLIPSOIDAL LONITUDE OF THE POINT MINUS THE       TMSFMC12
C*                    ELLIPSOIDAL LONGITUDE OF THE CENTRAL MERIDIAN OF  TMSFMC13
C*                    PROJECTION, (LONGITUDE POSITIVE EAST), IN RADIANS.TMSFMC14
C*             SFO - SCALE AT THE CENTRAL MERIDIAN.                     TMSFMC15
C*             A,B - SEMI-MAJOR AND SEMI-MINOR AXES OF THE REFERENCE    TMSFMC16
C*                   ELLIPSOID RESPECTIVELY, IN METRES.                 TMSFMC17
C*                                                                      TMSFMC18
C*       OUTPUT:                                                        TMSFMC19
C*             SF - POINT SCALE AT THE POINT.                           TMSFMC20
C*             C - MERIDIAN CONVERGENCE AT THE POINT, IN RADIANS.       TMSFMC21
C*                                                                      TMSFMC22
C*                                                                      TMSFMC23
C*                                     WRITTEN BY:                      TMSFMC24
C*                                         R.R. STEEVES, AUG., 1977     TMSFMC25
C*                                                                      TMSFMC26
C***********************************************************************TMSFMC27
      IMPLICIT REAL*8(A-Z)                                              TMSFMC28
      CP=DCOS(PHI)                                                      TMSFMC29
      T=DTAN(PHI)                                                       TMSFMC30
      ETA=DSQRT((A*A-B*B)/(B*B)*CP**2)                                  TMSFMC31
      C=DLAM*DSIN(PHI)*(1.D0+DLAM**2*CP**2/3.D0*(1.D0+3.D0*ETA**2+2.D0* TMSFMC32
     1  ETA**4)+DLAM**4*CP**4/15.D0*(2.D0-T**2))                        TMSFMC33
      SF=1.D0+DLAM**2*CP**2/2.D0*(1.D0+ETA**2)+DLAM**4*CP**4/24.D0*(5.D0TMSFMC34
     1   -4.D0*T**2)                                                    TMSFMC35
      SF=SF*SFO                                                         TMSFMC36
      RETURN                                                            TMSFMC37
      END                                                               TMSFMC38