IMPLICIT REAL*8(A-H,O-Z) 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 IS=1,2 DO IQ=1,2 DO IT=1,3 IF(IT.EQ.1) THEN TQM= 170.D0 ELSE IF(IT.EQ.2) THEN TQM= 175.D0 ELSE IF(IT.EQ.3) THEN TQM= 180.D0 ENDIF * OEXT= 'P' CALL TINIT(NRS,OEXT) IF(IS.EQ.2) THEN OVAL= 'A' OF3= 'OU9' CALL TCFLAG(OF3,OVAL) ENDIF IF(IQ.EQ.2) THEN OVAL= 'N' OF4= 'OU10' CALL TCFLAG(OF4,OVAL) ENDIF * 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)) IF(IS.EQ.1) THEN PRINT 4 ELSE PRINT 5 ENDIF IF(IQ.EQ.1) THEN PRINT 6 ELSE PRINT 7 ENDIF PRINT 3,TQM PRINT 1,OTPPO(1),OTPPO(35),OTPPO(36) PRINT 2,QW PRINT*,'---------------------------------------' ENDDO ENDDO ENDDO * 1 FORMAT(' W MASS = ',E20.6/ # ' C_1U = ',E20.6/ # ' C_1D = ',E20.6) 2 FORMAT(' Q_W = ',E20.6) 3 FORMAT(' M_T = ',E20.6) 4 FORMAT(' FACTORIZED ') 5 FORMAT(' ADDITIVE ') 6 FORMAT(' PERTURBATIVE ') 7 FORMAT(' NON-PERTURBATIVE ') * STOP END