I don't know what those constants are about. For the 1989 fortran code, apart from c, G, h, Pi , some values of the masses from other grpups are fed in from a table for comparison (e.g. the CODATA 1998 data): the output table then lists the program values alongside the CODATA values or also differences if that option is used.) E.g. searching for 4 in all the code gave:
grep 4 *.f
Hprog.f:C (Powerstation 4.0 under NT4.0 SP6 using an AMD K7 processor)
Hprog.f:C character*20 t1,t2,t3,t4,t5,t6,t7,t8
Hprog.f: integer*4 NN1,NN2 ! limits for NN
Hprog.f: integer*4 iL
Hprog.f: real*4 fma,GG4,SS4,FF4,FI4,sum14,sum24! result
Hprog.f: iprint = 0 ! no output of parts in eq.4
Hprog.f: Rg = 376.730313461 D0 ! CODATA'98 # const
Hprog.f: beta = 1.D0/1.00001411D0 ! #const
Hprog.f: call fmass(fma,GG4,SS4,FF4,FI4,WN0,sum14,sum24 )
Hprog.f:301 format(1x, a, 1pE14.7 , 11I4 )
Hprog.f: write(6,*) ' sum1=',sum14,' sum2=',sum24
Hprog.f: write(6,*)' G =',GG4, ' S=',SS4,' F=',FF4
Hprog.f: write(6,*)' FIFI =',FI4
Ibin.f: integer*4 Function Ibin(n,k)
Ibin.f: integer*4 N,K, ibi, ilo, ih ,iden ,ibinom ,i
Ibin.f:3 write(6,*) ' n over k is restricted to n <= 17 using integer*4'
WN0fu.f: integer*4 IBIN
WN0fu.f: integer*4 B,H,ieq, Lex
WN0fu.f: real*8 z3,z4,z5,z6,z7,z8,z9,z10,z11,z12,z13
WN0fu.f: real*8 teil1,teil2,teil3,teil4
WN0fu.f: + -dfloat((3*q-1)*(k-1)) + 0.5d0*dfloat((PP-QQ)*(4+(B+1)*(1-q))
WN0fu.f: z4 = dfloat(PP)*( 0.5d0*dfloat(B)+2.d0 + dfloat(q) )*dfloat(2-k)
WN0fu.f: z5=-dfloat(QQ)*(0.5d0*dfloat(B)+dfloat(1-4*(1+4*q)))*dfloat(2-k)
WN0fu.f: teil2 = z4 + z5
WN0fu.f: teil4 = ! # KLAMMER ergaenzt
WN0fu.f: + -0.25d0* dfloat(q)*dfloat((1-q)*(B-4))-0.25d0*dfloat(B-2)
WN0fu.f: a2 = teil1 -(1-r)*( teil2 + teil3 +teil4 )
WN0fu.f: z2 = (PP-QQ)*(4+(B+1)*(1-q))/2.
WN0fu.f: z4 = PP*(B/2.+2+q)*(2-k)
WN0fu.f: z5 = -QQ*(B/2.+1-4*(1+4*q))*(2-k)
WN0fu.f: z6 = (B-2)*(1+3*(PP-QQ)/2.)/4.
WN0fu.f: z9 = -(B+2)*(1-q)/4.
WN0fu.f: z11 = -q*(1-q)*(B-4)/4.
WN0fu.f: z12 = -(B-2)/4.
WN0fu.f: teil2 = z4+z5
WN0fu.f: teil4= -ibin(PP,3)*( (2*(1+ieq)+z10 -q)+z11+z12+z13 )
WN0fu.f: a2 = teil1 -(1-r)*( teil2 + teil3 +teil4 )
WN0fu.f: write(6,*) ' z3 =', z3 , ' z4 =', z4
WN0fu.f:C write(6,*) ' teil3=',teil3, ' teil4=',teil4
WN0fu.f: goto 400
WN0fu.f:c zw =(wet/dfloat(k))*(4.d0*(2.d0-wet)
WN0fu.f:c & *dfloat(4*B+PP+QQ))
WN0fu.f:cc zw =((wet/dfloat(k))*(4.d0*(2.d0-wet)
WN0fu.f:cc & *dfloat(4*B+PP+QQ))
WN0fu.f:c z4= ( dfloat((PP-QQ)*(H+2))+dfloat(PP )*( dfloat(5*B*(1+q)*QQ) +
WN0fu.f:c + ( z3 + z1 + z4 + z5 )
WN0fu.f:c write(6,*) ' z4 = ' ,z4
WN0fu.f:400 continue ! neuer code fuer y
WN0fu.f: zw = ( wet/k)*( 4*(2-wet)-pi*ebn*(1-eta)*wet)*(k+ebn*wet*(k-1))
WN0fu.f: + + 5*(1-q)*(4*B+PP+QQ)/(2*k+(-1)**k)
WN0fu.f: + -q*(1+ieq)*( k*(PP*PP+1)*(B+2)+(PP*PP+PP+1)/4.)
WN0fu.f: z4 = ( (PP-QQ)*(H+2)
WN0fu.f: zw1= z2+ibin(PP,2)*(1-ibin(QQ,3))* ( z3+z1 )+z4+z5!)
WN0fu.f: y = r*zw + (1-r)*zw1 ! KLAMMER versetzt 4.5.03
WN0fu.f:c write(6,*) ' z4 = ' ,z4
WN0fu.f: a3 = dfloat(4*B)*y/(1.d0+y)-1.d0/dfloat(4+B)
WN0fu.f: eque = a3/dfloat(4*B) ! alt1
compini.f:C integer*4 igam ,ialfpm ,iparm ,itab
compini.f: pi = 3.14159265358979 D0
compini.f: hq = 1.054571596 D-34 ! CODATA'98 (+-82)
compini.f: c = 2.99792458 D8
compini.f: ! # const: from table page 54 :
compini.f: qn(2) =24
compini.f: qp(2) =34
compini.f: HH(2)=104
compini.f: oc = 4.D0/3.D0 ! # const
compini.f:C Rg = 376.730313461 D0 ! CODATA'98 # const c.f. Hprog.for
compini.f:C beta = 1.D0/1.00001411D0 ! #const "
compini.f: case (4)
compini.f:C fakMeV = 0.05609545 D31 (Sch)
compini.f: eta = pi/Dsqrt(Dsqrt(4.D0+ pi*pi*pi*pi))
compini.f: eq = 3.D0/(4.D0*pi*pi)*Dsqrt(2.D0*theta*hq/Rg)
compini.f: beta = 1.D0/1.00001411D0
compini.f:C eq.4a page 13 ( anno 1985 ? )
compini.f: write(6,101) 'c.f. Anhang B , pg. 54 '
compini.f: goto 400
compini.f:C (4-8) Parametermatrix: entfaellt
compini.f:400 continue
detailini.f: integer*4 jq,jk
detailini.f: do jq = 0,2; n4tab(jk,jq)= -99999.; enddo
detailout.f: integer*4 jq,jk
detailout.f: write(6,*) ' N4(k,q) ,k =', jk, ' q = 0,1,2 :'
detailout.f: write(6,*)( n4tab(jk,jq), jq = 0,2)
etaqk.f:C pi = 3.14159265358979 D0 I am Common/CONST/
etaqk.f: integer*4 kq,k
etaqk.f: zw1 = Dsqrt(kq*kq*kq*kq*(4.d0+k)+pi*pi*pi*pi)
fmass.f: Subroutine Fmass(fma,GG4,SS4,FF4,FI4,WN04,sum14,sum24 )
fmass.f:C + , N1,N2,N3,N4,EQK,EQ1,E1K
fmass.f: + ,zw1,zw2,zw3 ,zw4,zw5,AAA ,UU ,sum1,sum2
fmass.f: real*4 fma, WN04,GG4,SS4,FF4,FI4 ,sum14,sum24
fmass.f: integer*4 ibin
fmass.f: zw4 = - zw3*zw1* oc/u
fmass.f:C write(6,*) 'zw5,zw4 =',zw5,zw4
fmass.f: zw4 = zw5 + zw4
fmass.f:C write(6,*) ' ln(0.5*k*N3) = ' ,zw4 , ' q,k=',q,k
fmass.f: N3 = dexp( zw4) * 2.d0/dfloat(k)
fmass.f: N4 = dfloat( 4*( 1 + q*(k-1)) / k )
fmass.f: zw4 = (1.d0 -dsqrt(EQK))/(1.d0+dsqrt(EQK))
fmass.f: zw4 = zw4 * zw4
fmass.f: N5 = AAA*(1.d0+dfloat(k*(k-1)*2**(k*k+3))*AAA*fuNk(k)*zw4)
fmass.f: zw4 = 4.d0*(1.d0 - dsqrt(eta))/(1.d0 + dsqrt(eta))
fmass.f: zw4 = zw4*zw4
fmass.f: zw3 = eta*(1.d0 -alfm/alfp)*zw4*dfloat(qs(k))
fmass.f: N6 = 4.d0*dfloat(k)*
fmass.f: + + dfloat(4*r*BB(k)*(1-QQ)) /dfloat(3-2*q)
fmass.f: + -dfloat((PP-QQ)*(1-q))*4.d0*pi/dsqrt(dsqrt(2.d0)) ) )
fmass.f: fi=N4*dfloat(p*p/(1+p*p))*(dfloat(s+qs(k))/dsqrt(dfloat(1+s*s)))*
fmass.f: + ( dsqrt(dsqrt(2.d0)) -4.d0*UU*dfloat(BB(k))/WN0f )
fmass.f: + + dfloat(p*(p+1))*N3 + dfloat(4*s)
fmass.f: + + dfloat( qp(k)*(1+qp(k)))*N3 + dfloat(4*qs(k))
fmass.f: sum2 = amu* 4.d0* dfloat(q) * alfm
fmass.f:C Fma = mue*((GG + SS + FF + FI)*alfp + 4.*q * alfm )
fmass.f: GG4 = GG ; SS4 = SS; FF4 = FF; FI4 = FIFI
fmass.f: sum14 =sum1; sum24 =sum2 ; Wn04 = Wn0f
fmass.f: n4tab(k,q) = N4 ;n5tab(k,q) = N5 ;n6tab(k,q) = N6
fmass.f: integer*4 k , nk , nkfu
fmass.f: integer*4 function NSk (k) ! eq. 8f1 pg.15
fmass.f: integer*4 k , nk
thetaqk.f: integer*4 q ,k