C     SccsID = "@(#)azmeqd.f	1.1     07/23/03"
      SUBROUTINE AZMEQD (ICODE,FLAT,FLON,X,Y,CC,S)
      IMPLICIT REAL*8 (A-H,O-Z)
      CHARACTER*4 IGP,IPC,ICODE
      DIMENSION AED(18)
      DATA AED/5.0D4,3.092241724D1,5.2109550254D5,6.7686580D-3,
     1  3.087002482D1,4.82632837578D4,1.052893943D3,4.483386D0,2.3559D-2
     2 ,1.27564128D-1,4.82632837702D4,3.239388390D-2,1.87770D0,
     3  1.047546691D3,6.193011D0,5.0699D-2,5.2109550254D5,3.092241724D1/
      DATA RHOSEC/2.062648062471D05/
      DATA IPC/'PCGP'/
      IF (ICODE.EQ.IPC) GO TO 20
C
      SINP = DSIN(FLAT)
      COSP = DCOS(FLAT)
      TANP = DTAN(FLAT)
      SQ = DSQRT(1.0 - AED(4)*SINP**2)
      X = AED(1) + AED(2)*(FLON*RHOSEC-AED(3))*COSP/SQ
      X1 = (X-AED(1))**2/1.0D8
      T1 = FLAT*RHOSEC - AED(6) - (AED(7)-AED(8)*COSP**2+AED(9)*COSP**4)
     1  *SINP*COSP
      T2 = X1*TANP*SQ/AED(10)
      Y = AED(1) + AED(5)*T1 + T2
      CC = 0.
      S = 1.
      RETURN
C
   20 X1 = (X-AED(1))**2/1.0D8
      DO 30 I=1,3
        IF (I.GT.1) GO TO 25
        T2 = X1*AED(13)
        GO TO 28
   25   T2 = X1*TANP*SQ/AED(10)
   28   OMEGA = AED(11) + AED(12)*(Y-AED(1)-T2)
        SINM = DSIN(OMEGA/RHOSEC)
        COSM = DCOS(OMEGA/RHOSEC)
        PSEC = OMEGA + (AED(14) + AED(15)*COSM**2 + AED(16)*COSM**4)
     1    *SINM*COSM
        FLAT = PSEC/RHOSEC
        SINP = DSIN(FLAT)
        TANP = DTAN(FLAT)
        SQ = DSQRT(1.0 - AED(4)*SINP**2)
   30   CONTINUE
      COSP = DCOS(FLAT)
      ESEC = AED(17) + (X-AED(1))*SQ/AED(18)/COSP
      FLON = ESEC/RHOSEC
      RETURN
      END