IMPLICIT REAL*8(A-H,O-Z)
*
*-----NOTE CHANGE IN NPO = 34 -> NPO = 36
*
      PARAMETER(NX=1,NY=2,NPO=36,NRO=26*NX,NFIT=7*NY)
      CHARACTER*1 OEXT,OVAL,ORAD,OCHAN
      CHARACTER*2 OV2
      CHARACTER*3 OF3
      CHARACTER*4 OF4
      CHARACTER*5 OF5
      CHARACTER*6 OF6
      CHARACTER*7 OF7
*
      COMMON/TK/IFK
      COMMON/TFITO/APO(9)
      COMMON/TFIT/FIT(NFIT)
*
      DIMENSION RS(NX)
      DIMENSION OTPPO(NPO),OTPRO(NRO),OTPROD(9,NX),OTPRODD(9,NX),
     #          OTPRODDD(9,NX),OTPROZZ(9,NX),OTPROWI(9,NX),
     #          OTPDRO(5,9,NX)
      DIMENSION SET(NY),AET(NY),SMUT(NY),AMUT(NY),STAUT(NY),
     #          ATAUT(NY),SHADT(NY),ESET(NY),EAET(NY),
     #          ESMUT(NY),EAMUT(NY),ESTAUT(NY),EATAUT(NY),
     #          ESHADT(NY),ACT(NY),ABT(NY),EACT(NY),EABT(NY),
     #          SBT(NY),ESBT(NY),SCT(NY),ESCT(NY),BRAT(NY),
     #          EBRAT(NY),CRAT(NY),ECRAT(NY),DECXSH(NX),DDECXSH(NX),
     #          DDDECXSH(NX),ZZECXSH(NX),WIECXSH(NX),HMI(5),
     #          FET(NY),BET(NY)
*
      IFK= 0
      NRS= 1
      NRST= 2
      ZM= 91.1875D0
      ALS= 0.119D0
      HM= 150.D0
*
      RS(1)= 92.D0
*
      DO IO=1,3
       IF(IO.EQ.1) THEN
        OVAL= 'L'
       ELSE IF(IO.EQ.2) THEN
        OVAL= 'C'
       ELSE IF(IO.EQ.3) THEN
        OVAL= 'R'
       ENDIF

      DO IT=1,11
       DO IH=1,11
        TQM= 170.D0+(IT-1)*1.D0
        HM= 100.D0+(IH-1)*10.D0
        PRINT 3,OVAL,TQM,HM 
*

      OEXT= 'P'
      CALL TINIT(NRS,OEXT)
      OF3= 'OU8'
      CALL TCFLAG(OF3,OVAL)
*
      SE= 1.D-1
      CALL TOPAZ0(SE,NRS,RS,ZM,TQM,HM,ALS,OTPPO,OTPRO,OTPDRO)
*
      Z= 55.D0
      RN= 78.D0
*
      QW= -2.D0*(OTPPO(35)*(2.D0*Z+RN)+OTPPO(36)*(Z+2.D0*RN))
      PRINT 1,(OTPPO(I),I=1,NPO)
      PRINT 2,QW
*
    1 FORMAT( '  W MASS   (GEV)        = ',E20.6/
     #        '  NU                    = ',E20.6/
     #        '  ELECTRON              = ',E20.6/
     #        '  MUON                  = ',E20.6/
     #        '  TAU                   = ',E20.6/
     #        '  UP                    = ',E20.6/
     #        '  DOWN(STRANGE)         = ',E20.6/
     #        '  CHARM                 = ',E20.6/
     #        '  BOTTOM                = ',E20.6/
     #        '  SIN^2(E)              = ',E20.6/
     #        '  SIN^2(B)              = ',E20.6/
     #        '  A_FB(L) EFF.          = ',E20.6/
     #        '  A_LR EFF.             = ',E20.6/
     #        '  TOTAL WIDTH (GEV)     = ',E20.6/
     #        '  G_H/G_E               = ',E20.6/
     #        '  SIGMA0_H  (NB)        = ',E20.6/
     #        '  G(B)/G(HAD)           = ',E20.6/
     #        '  A_FB(B)               = ',E20.6/
     #        '  HADRONIC WIDTH (GEV)  = ',E20.6/
     #        '  INVISIBLE             = ',E20.6/
     #        '  A_FB(C)               = ',E20.6/
     #        '  G(C)/G(HAD)           = ',E20.6/
     #        '  A_LR(B) EFF.          = ',E20.6/
     #        '  A_LR(C) EFF.          = ',E20.6/
     #        '  SIN^2(C)              = ',E20.6/
     #        '  RHO(E)                = ',E20.6/
     #        '  RHO(C)                = ',E20.6/
     #        '  RHO(B)                = ',E20.6/
     #        '  1/ALPHA^5             = ',E20.6/
     #        '  1/ALPHA               = ',E20.6/
     #        '  SIGMA0_L  (NB)        = ',E20.6/
     #        '  ALPHA_S(M_T)          = ',E20.6/
     #        '  G_H/G_M               = ',E20.6/
     #        '  G_H/G_T               = ',E20.6/
     #        '  C_1U                  = ',E20.6/
     #        '  C_1D                  = ',E20.6)
    2 FORMAT( '  Q_W(Cs)               = ',E20.6)
    3 FORMAT(1X,A1,1X,2E20.5)
*
        ENDDO
       ENDDO
      ENDDO
*
      STOP
      END