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