diff --git a/src/conex_mod.F b/src/conex_mod.F
index 29decc9a048469d15691282aa462d9a5a30e132d..122d893d1268c9c20175df14fa74786affe31041 100644
--- a/src/conex_mod.F
+++ b/src/conex_mod.F
@@ -48987,8 +48987,8 @@ c external type declarations
c COMMONs
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER LDIFF
COMMON /S_CLDIF/ LDIFF
C The final particle output is contained in COMMON /S_PLIST/
@@ -49046,7 +49046,7 @@ c internal type declarations
DATA LL /5*0,7*2,2*1,12*0,2,6*0,6*1,19*0,2,2,10*0,
& 2,2,0,2,2,11*0,1,1,1,9*0,1/
- if(Ndebug.gt.0)then
+ if(Idebug.gt.0)then
WRITE(LUN,'(A42,I3,I3,1X,F10.2)')
& ' SIBYLL: called with (K_beam,IATARG,Ecm):',
& K_beam,IATARG,Ecm
@@ -49095,14 +49095,14 @@ C...Diffractive/non-diffractive interactions
ENDIF
IF (IREJ.NE.0) THEN
- if(Ndebug.gt.0) WRITE(LUN,'(A38,F10.2,I3,I3,I3)')
+ if(Idebug.gt.0) WRITE(LUN,'(A38,F10.2,I3,I3,I3)')
& ' SIBYLL: rejection (Ecm,Ncall,Nw,JDIF):',Ecm,Ncall,NW,JDIF(1)
GOTO 100
ENDIF
do J=1,NP
if (P(J,4).lt.0.D0 ) then
- if(Ndebug.gt.0)then
+ if(Idebug.gt.0)then
WRITE(LUN,*)' negative energy particle!' , P(J,4)
CALL SIB_LIST(LUN)
endif
@@ -49125,12 +49125,12 @@ c stop
goto 100
ENDIF
IF (ABS(PZsum+0.5D0*Ecm*DBLE(NW-1)) .GT. 0.1D0) THEN
- if(Ndebug.gt.0)THEN
+ if(Idebug.gt.0)THEN
WRITE(LUN,*) ' SIBYLL: momentum not conserved (L,call): ',
& L,Ncall
WRITE(LUN,*) ' pz_inp = ', 0., ' pz_out = ', pzsum
ENDIF
- IF(ndebug.gt.0)then
+ IF(idebug.gt.0)then
CALL PRNT_PRTN_STCK
CALL SIB_LIST(LUN)
WRITE(LUN,*) ' SIBYLL: event rejected'
@@ -49155,7 +49155,7 @@ c exchange pi0 with charged pions for meson projectiles
C...list final state particles
- if(Ndebug.gt.10) CALL SIB_LIST(LUN)
+ if(Idebug.gt.10) CALL SIB_LIST(LUN)
END
@@ -49838,8 +49838,8 @@ c commons
INTEGER IPAR
COMMON /S_CFLAFR/ PAR(NPAR_max), IPAR(NIPAR_max)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION STR_mass_val, STR_mass_val_hyp, STR_mass_sea
COMMON /S_CUTOFF/ STR_mass_val, STR_mass_val_hyp, STR_mass_sea
@@ -49867,17 +49867,17 @@ c commons
c set standard parameters (full set)
CALL PAR_INI
c read new parameters from file
- IF(ndebug.gt.0)WRITE(LUN,*)'reading parameter file: sibyll.par'
+ IF(idebug.gt.0)WRITE(LUN,*)'reading parameter file: sibyll.par'
DO WHILE (istat.ge.0)
READ(4,14,iostat=ISTAT) CNAME,NUMBER
IF(CNAME.eq.'IPAR ')THEN
READ(NUMBER,*) j, ival
- IF(ndebug.gt.1)write(LUN,15) 'IPAR(',j,')=', ival
+ IF(idebug.gt.1)write(LUN,15) 'IPAR(',j,')=', ival
IPAR(J) = iVAL
ELSEif(CNAME.eq.'PAR ')THEN
READ(NUMBER,*) j, val
PAR(J) = VAL
- IF(ndebug.gt.1)write(LUN,16) ' PAR(',j,')=', val
+ IF(idebug.gt.1)write(LUN,16) ' PAR(',j,')=', val
ELSE
WRITE(LUN,*)'wrong format in parameter file!'
WRITE(6,*)'wrong format in parameter file!'
@@ -49936,8 +49936,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -50031,7 +50031,7 @@ c last in use:93
C...Fragmentation of nuclei
DATA KODFRAG /0/
C...Debug label and event counter
- DATA Ndebug /0/
+ DATA Idebug /0/
DATA Ncall /0/
END
@@ -50046,7 +50046,7 @@ C-----------------------------------------------------------------------
COMMON /S_CZLEAD/ CLEAD, FLEAD
COMMON /S_CPSPL/ CCHIK(4,99)
- COMMON /S_DEBUG/ Ncall, Ndebug, Lunn
+ COMMON /S_DEBUG/ Ncall, Idebug, Lunn
INTEGER NW_max
PARAMETER (NW_max = 20)
C--------------------------------------------------------------------
@@ -50179,7 +50179,7 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- COMMON /S_DEBUG/ Ncall, Ndebug, Lunn
+ COMMON /S_DEBUG/ Ncall, Idebug, Lunn
INTEGER NW_max
PARAMETER (NW_max = 20)
C--------------------------------------------------------------------
@@ -50434,8 +50434,8 @@ C----------------------------------------------------------------
IMPLICIT INTEGER(I-N)
COMMON /S_PDG2PID/ IPID_PDG_LIST(99),ID_LIST(577)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION CBR
INTEGER KDEC,LBARP,IDB
@@ -50446,7 +50446,7 @@ C----------------------------------------------------------------
Nin = abs(Npdg)
if((Nin.gt.999999).or.(Nin.eq.0)) then
C invalid particle number
- if(ndebug.gt.5) write(6,'(1x,A,I10)')
+ if(idebug.gt.5) write(6,'(1x,A,I10)')
& ' ISIB_PDG2PID: invalid PDG ID number ',Npdg
ISIB_PDG2PID = 0
return
@@ -50462,7 +50462,7 @@ C use hash algorithm
C particle not in table
if(ID_list(Nout).Eq.0) then
- if(ndebug.gt.0) write(6,'(1x,A,I10)')
+ if(idebug.gt.0) write(6,'(1x,A,I10)')
& ' ISIB_PDG2PID: particle not in table ',Npdg
ISIB_PDG2PID = 0
return
@@ -50494,8 +50494,8 @@ C=======================================================================
C-----------------------------------------------------------------------
C PDG conversion blocks \FR'13
C----------------------------------------------------------------
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER ( ID_PDG_MAX = 260 )
COMMON /S_PDG2PID/ ID_PDG_LIST(99),ID_LIST(577)
SAVE
@@ -50510,7 +50510,7 @@ C----------------------------------------------------------------
& -423,0,443,4222,4212,4112,4232,4132,4122,-15, !90
& 15,-16,16,4224,4214,4114,4324,4314,4332/
- IF(Ndebug.gt.2)
+ IF(Idebug.gt.2)
& WRITE(LUN,*) ' INITIALIZING PDG TABLES..'
CALL SIB_CPCINI(ID_pdg_max,ID_pdg_list,ID_list)
@@ -50534,8 +50534,8 @@ C
C***********************************************************************
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
integer Number(*),List(*),Nrows
Integer Nin,Nout,Ip,I
SAVE
@@ -50562,7 +50562,7 @@ C Calculate a list number for this particle id number
If(Nout.Lt.0) Then
C Count the bad entries
- IF(Ndebug.gt.3) Write(LUN,'(1x,a,i10)')
+ IF(Idebug.gt.3) Write(LUN,'(1x,a,i10)')
& ' SIB_CPCINI: invalid particle ID',Nin
Go to 500
End If
@@ -50570,7 +50570,7 @@ C Count the bad entries
List(Nout) = Ip
Else
If(Nin.eq.Number(List(Nout))) Then
- IF(Ndebug.gt.3)Write(LUN,'(1x,a,i10)')
+ IF(Idebug.gt.3)Write(LUN,'(1x,a,i10)')
& ' SIB_CPCINI: double particle ID',Nin
End If
Nout = Nout + 5
@@ -50731,8 +50731,8 @@ Cf2py integer,intent(out) :: kf
INTEGER IPAR
COMMON /S_CFLAFR/ PAR(NPAR_max), IPAR(NIPAR_max)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -50767,7 +50767,7 @@ C--------------------------------------------------------------------
+ 246,256,271,0,0,1,3,6,10,15,21,60,61,64,70,292,307,328,356,
+ 0,1,3,6,10,15,21,28,16*0/
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*)' SIB_FLAV: input:',IFL1, IFL2_A, IRNK, IFL2, KF
c set rho0 / ( omega, phi ) ratio, i.e. I=1 to I=0
@@ -50911,7 +50911,7 @@ c suppress rank1 (leading) omega
ENDIF
ENDIF
c PRINT*,' I4FLAV returns :(IFL1,IFL2,LL)',IFL1,IFL2,KF
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*)' SIB_FLAV: output:',IFL1, IFL2_A, IRNK, IFL2, KF
RETURN
ENDIF
@@ -51006,7 +51006,7 @@ c i.e. forbid: Blead-Bbar-B combination
IFL2 = IFL2 + ISIGN(100,IFL2)
endif
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*)' SIB_FLAV: output:',IFL1, IFL2_A, IRNK, IFL2, KF
RETURN
@@ -51032,8 +51032,8 @@ Cf2py integer,intent(out) :: ifl
IMPLICIT NONE
DOUBLE PRECISION Q2
INTEGER IFL,IS
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
C--------------------------------------------------------------------
@@ -51055,7 +51055,7 @@ C--------------------------------------------------------------------
DOUBLE PRECISION XMS2,XMC2,P_S,P_C,S_RNDM,QMASS,FERMI
INTEGER IFL1,IS0
- IF( NDEBUG.gt.6 )
+ IF( IDEBUG.gt.6 )
& WRITE(LUN,*)' SIB_ICFLAV: input (Q2,IFL,IS):',Q2,IFL,IS
c quark or antiquark, sampled if input is zero
@@ -51082,7 +51082,7 @@ c P_c: 0 (s only) --> 1 (s,c) equal
& + PAR(155) * FERMI( Q2, XMC2, -PAR(152) )
P_C = PAR(156) * 0.5D0*FERMI( Q2, XMC2, -PAR(152) )
ENDIF
- IF(NDEBUG.gt.6)THEN
+ IF(IDEBUG.gt.6)THEN
WRITE(LUN,*)' SIB_ICFLAV: (4*M_S**2, P_S, kT):',
& xms2, P_s, PAR(151)
WRITE(LUN,*)' SIB_ICFLAV: (4*M_C**2, P_C, kT):',
@@ -51097,7 +51097,7 @@ c replace s with c
IFL = IS*IFL1
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*)' SIB_ICFLAV: output (Q2,IFL,IS):',Q2,IFL,IS
END
@@ -51117,8 +51117,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -51178,7 +51178,7 @@ c charmed mesons, charmed baryons
DATA (KK(I), I= 1,39) /5*0,3*2,4*3,2*1,6*0,6*2,3,6*2,6*4/
DATA (KK(I), I=40,99) /19*0,5,5,10*0,5,5,0,5,5,11*0,6,6,6,9*0,6/
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& WRITE(LUN,*)' SIB_DIFF: called with (L0,JDIF1,Ecm):',
& L0,JDIF1,Ecm
@@ -51207,7 +51207,7 @@ c add interactions
LL = L0
LA = IABS(L0)
XM2MAX = PAR(13)*Ecm*Ecm
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& WRITE(LUN,*)' SIB_DIFF: max diff. mass (M,Xi):',XM2MAX,PAR(13)
C...Double diffraction
@@ -51231,7 +51231,7 @@ c add diffractive system to parton stack
CALL ADD_PRTN_4VEC(P2,-3,0,0,Iref)
CALL ADD_INT_REF(Iref,1)
ENDIF
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& write(lun,*)' double-diff.: (kb,xmb,kt,xmt)',LL,xmb,13,xmt
CALL DIFDEC (LL, Irec, IDBAD, P1)
IF(IDBAD.eq.1)goto 20
@@ -51251,7 +51251,7 @@ C...Single diffraction
EM2 = AM2(13)
L = 13
ZD = -1.D0
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& write(lun,*)' single-diff. (beam): (kb)',LL
ELSE
K = 1
@@ -51260,7 +51260,7 @@ C...Single diffraction
L = LL
LL = 13
ZD = +1.D0
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& write(lun,*)' single-diff. (target): (kt)', LL
ENDIF
@@ -51326,8 +51326,8 @@ C=======================================================================
C----------------------------------------------------------------------
IMPLICIT NONE
INTEGER I,NPION
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION XM2MIN,ALXMIN,SLOP0,ASLOP,BSLOP,XMASS
COMMON /S_DIFMAss/ XM2MIN(6),ALXMIN(6),SLOP0,ASLOP,BSLOP,XMASS(2)
INTEGER NIPAR_max,NPAR_max
@@ -51365,7 +51365,7 @@ C charmed baryons (6), lowest mass: lambda_c
ALXMIN(6) = log(XM2MIN(6))
c debug output
- IF(NDEBUG.gt.1)THEN
+ IF(IDEBUG.gt.1)THEN
WRITE(LUN,*)'DIFF_INI: setting diff. mass parameters'
WRITE(LUN,*)' min mass: ', (XM2MIN(I), I=1,6)
WRITE(LUN,*)' log min mass: ', (ALXMIN(I), I=1,6)
@@ -52624,8 +52624,8 @@ C-----------------------------------------------------------------------
IMPLICIT NONE
INTEGER IDX,LBAD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -52701,13 +52701,13 @@ c read string 4momentum from stack
C kinematic variables
EE = PST(5) ! string mass
- IF(NDEBUG.gt.1) WRITE(LUN,*)' FRAG_VLNCE: IDX,EE,IFLB,IFLT',
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' FRAG_VLNCE: IDX,EE,IFLB,IFLT',
& IDX,EE,IFLB,IFLT
IF(IDX.ne.ISTH) then
c read merged string and add hadron to final particle stack..
CALL RD_PRTN_4VEC(ISTH,PstH,LSTH,IST)
- IF(NDEBUG.gt.1) WRITE(LUN,*)' FRAG_VLNCE: found merged string',
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' FRAG_VLNCE: found merged string',
& LSTH,(PSTH(I),I=1,5)
IF(IDX.ne.IST) then
write(lun,*) ' FRAG_VLNCE: reference loop broken!' , IDX
@@ -52739,7 +52739,7 @@ c charm fractions in different parameterizations
PAR(24) = PAR(25)*EXP(-PAR(26)/EE)
ENDIF
- IF(NDEBUG.gt.1)
+ IF(IDEBUG.gt.1)
& WRITE(LUN,*)' FRAG_VLNCE: parameters (CHM,DIQ,STR,VEC,POP)',
& PAR(24),PAR(1),PAR(2),PAR(5),PAR(8)
@@ -52768,7 +52768,7 @@ C... fragment strings in string restframe
PAR(1) = PAR1_def
KINT= 0
IF (IFBAD .EQ. 1) then
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& WRITE(LUN,*)' STRING_FRAG: rejection (Ncall):',Ncall
RETURN
ENDIF
@@ -52789,7 +52789,7 @@ c to calculate rotation angles in string center-of-mass
& P2TOT,p2(1),p2(2),p2(3),p2(4))
c should be back-to-back...
- IF(ndebug.gt.1)THEN
+ IF(idebug.gt.1)THEN
write(lun,*)
& ' FRAG_VLNCE: string c.m. momentum, parton 1 (Pabs,P(i)):' ,
& P1TOT, (P1(j),j=1,4)
@@ -52879,8 +52879,8 @@ C IAT : mass number of target
INTEGER KB,IAT,KT
COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -52978,8 +52978,8 @@ C... normalisation for Peterson/SLAC frag. func
DOUBLE PRECISION ZDMAX, EPSI
COMMON /S_CZDISc/ ZDMAX, EPSI
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
SAVE
c get the maximum zmefn value first for normalisation
@@ -52990,7 +52990,7 @@ c get the maximum zmefn value first for normalisation
z = dble(j)/dble(jmax+1)
zdmax = max(zdmax, zmefn(z,epsi))
enddo
- if (ndebug .gt. 0) WRITE(LUN,*)' ZDMAX,EPS:',zdmax, epsi
+ if (idebug .gt. 0) WRITE(LUN,*)' ZDMAX,EPS:',zdmax, epsi
RETURN
END
C=======================================================================
@@ -53964,8 +53964,8 @@ c DOUBLE PRECISION SQS,S,PTmin,XMIN,ZMIN
c INTEGER KB,IAT,KT
c COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
c
-c INTEGER NCALL, NDEBUG, LUN
-c COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+c INTEGER NCALL, IDEBUG, LUN
+c COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
c
c DOUBLE PRECISION PPT02
c COMMON /S_CQDIS2/ PPT02(44)
@@ -54058,8 +54058,8 @@ c DOUBLE PRECISION SQS,S,PTmin,XMIN,ZMIN
c INTEGER KB,IAT,KT
c COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
c
-c INTEGER NCALL, NDEBUG, LUN
-c COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+c INTEGER NCALL, IDEBUG, LUN
+c COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
c
c DOUBLE PRECISION PPT02
c COMMON /S_CQDIS2/ PPT02(44)
@@ -54146,8 +54146,8 @@ c DOUBLE PRECISION SQS,S,PTmin,XMIN,ZMIN
c INTEGER KB,IAT,KT
c COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
c
-c INTEGER NCALL, NDEBUG, LUN
-c COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+c INTEGER NCALL, IDEBUG, LUN
+c COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
c
c DOUBLE PRECISION PPT02
c COMMON /S_CQDIS2/ PPT02(44)
@@ -54239,8 +54239,8 @@ c DOUBLE PRECISION SQS,S,PTmin,XMIN,ZMIN
c INTEGER KB,IAT,KT
c COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
c
-c INTEGER NCALL, NDEBUG, LUN
-c COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+c INTEGER NCALL, IDEBUG, LUN
+c COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
c
c DOUBLE PRECISION PPT02
c COMMON /S_CQDIS2/ PPT02(44)
@@ -54314,8 +54314,8 @@ Cf2py double precision, intent(out) :: PT
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NW_max
PARAMETER (NW_max = 20)
C--------------------------------------------------------------------
@@ -54367,7 +54367,7 @@ C--------------------------------------------------------------------
ZSOF = 2.D0*LOG(STR_mass_min/SQS) ! minim. mass ~ x1 * x2
XMINA = MAX(EXP(ZSOF),EPS10)
AXMINA = 1.D0/XMINA
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*) ' SAMPLE_SOFT6: Mmin,ZSOF,XMINA,XMAX,SLOPE:',
& STR_mass_min,ZSOF,XMINA,XMAX,SLOPE
@@ -54375,7 +54375,7 @@ C--------------------------------------------------------------------
IF(NOSLOPE.eq.1) goto 200
XRNDM = S_RNDM(0)
XR = LOG(1.D0-X1)-LOG(1.D0-XMINA)
- IF(ndebug.gt.5)
+ IF(idebug.gt.5)
& write(lun,*) ' X1,XR,SLOPE*XR:',X1,XR,SLOPE*XR
if(SLOPE*XR.le.LOG(max(xrndm,eps10))) goto 100
@@ -54383,7 +54383,7 @@ C--------------------------------------------------------------------
IF(NOSLOPE.eq.1) goto 300
XRNDM = S_RNDM(1)
XR = log(1.D0-X2) - log(1.D0-XMINA)
- IF(ndebug.gt.5)
+ IF(idebug.gt.5)
& write(lun,*) ' X2,XR,SLOPE*XR:',X2,XR,SLOPE*XR
if(SLOPE*XR.le.log(max(xrndm,eps10))) goto 200
@@ -54393,7 +54393,7 @@ C--------------------------------------------------------------------
STR_mass2 = sqrt(X1*X2*S)/2.D0
PPTT = PPT02(10)
IF(IPAR(3).eq.8) PPTT = PPT02(20)
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*) ' SAMPLE_SOFT6: PPTT,Mmin2,PTmin:',
&PPTT,STR_mass2,PTmin
150 PT = PPTT*SQRT(-LOG(MAX(EPS10,S_RNDM(0))))
@@ -54405,7 +54405,7 @@ C--------------------------------------------------------------------
XMT2 = XMT**2
PT = SQRT(XMT2-XM2)
ENDIF
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*) ' XM,XMT2,PT:',XM,XMT2,PT
IF(PT.GT.PTmin) GOTO 150
IF(PT.GE.STR_mass2) GOTO 150
@@ -54434,8 +54434,8 @@ C.
C-----------------------------------------------------------------------
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
c external type declarations
INTEGER NW_max,JDIF,IA,L,IAFLG,NW
DOUBLE PRECISION SQS
@@ -54462,7 +54462,7 @@ c local type declarations
INTEGER K
SAVE
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
&WRITE(LUN,*)'SIB_START_EV:', SQS, L, IA, IAFLG, NW, JDIF
C...sample number of wounded nucleons
@@ -54498,7 +54498,7 @@ c sample number of wounded nucleons
ENDIF
NW = NA
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& WRITE(LUN,'(A50,2I3,1P,3E10.3)')
& ' START_EVT: IA, NW, SIGT, SLOPE, RHO:',IA,NW,SIGT,SLOPE,RHO
C...new treatment of diffraction
@@ -54576,8 +54576,8 @@ c COMMONs
INTEGER NS_max, NH_max
PARAMETER (NS_max = 20, NH_max = 80)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -54659,7 +54659,7 @@ c local types
DATA FOX /0.21522D0/ !atomic percentage of 'non-nitrogen' in air
- IF(NDEBUG.gt.0.and.IMOD.eq.1)
+ IF(IDEBUG.gt.0.and.IMOD.eq.1)
& WRITE(LUN,'(A50,F10.2,I4,I3,I3)')
& ' INI_EVENT: called with (ECM,KBEAM,IATARG,NCALL):',
& ECM,KBEAM,IATARG,NCALL
@@ -54712,7 +54712,7 @@ C... Generate an 'air' interaction by choosing Nitrogen or Oxygen
R = S_RNDM(0)
IATARG = 14
IF (R .LT. FOX) IATARG = 16
- if (NDEBUG.gt.0)
+ if (IDEBUG.gt.0)
* WRITE(lun,*)'fox,rndm,iatarg,eps:',fox,r,iatarg,eps8
ENDIF
DO JN=1,IATARG
@@ -54733,7 +54733,7 @@ c...EJA correction 2007.03.27
endif
XMIN = 4.D0*PTmin**2/S
ZMIN = LOG(XMIN)
- IF(ndebug.gt.0)then
+ IF(idebug.gt.0)then
write(lun,*) ' INI_EVENT: ncall:', ncall
write(lun,'(2X,A33,F10.2,1X,F16.2,F8.5,E10.3,F10.5)')
& 'INI_EVENT: (SQS,S,PTmin,Xmin,Zmin)',
@@ -54774,8 +54774,8 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -54796,7 +54796,7 @@ c level0 index
NPP0 = NPP0 + 1
ENDIF
IREFout = NPP
- IF(NDEBUG.gt.6)THEN
+ IF(IDEBUG.gt.6)THEN
WRITE(LUN,*) ' ADD_PRTN: (#,PID,LEVEL,REF)',
& NPP,LPID(NPP,1),LPID(NPP,2),LPID(NPP,3)
WRITE(LUN,*) ' 4momentum: ',(PP(NPP,JJ),JJ=1,5)
@@ -54813,8 +54813,8 @@ C----------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DIMENSION PIN(5)
SAVE
@@ -54834,8 +54834,8 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -54844,7 +54844,7 @@ C-------------------------------------------------------
c IF(LPID(IDX,3).ne.0) WRITE(LUN,*)
c & ' ADD_REF: warning particle already has defined reference,',
c & IDX,' overwritting..'
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
&WRITE(LUN,*) ' ADD_REF: (IDX,REFin)',IDX,Irefin
LPID(IDX,3) = Irefin
END
@@ -54860,15 +54860,15 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
Irefout = LPID(IDX,3)
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*) ' RD_ref: (IDX,REFout)',IDX,Irefout
END
@@ -54883,14 +54883,14 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*) ' ADD_INT_REF: (IDX,REFin)',IDX,Irefin
LPID(IDX,4) = Irefin
END
@@ -54906,8 +54906,8 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -54915,7 +54915,7 @@ C-------------------------------------------------------
Irefout = LPID(IDX,4)
IF(Irefout.ne.0) Iout = LPID(Irefout,1)
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*) ' RD_INT: (IDX,REFout,Iint)',IDX,Irefout,Iout
END
@@ -54929,14 +54929,14 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
- IF(NDEBUG.gt.6)THEN
+ IF(IDEBUG.gt.6)THEN
WRITE(LUN,*) ' EDT_PRTN: (#,PID,LEVEL,REF)',
& IDX,LPID(IDX,1),LPID(IDX,2),LPID(IDX,3)
WRITE(LUN,*) ' initial 4momentum:',(PP(IDX,JJ),JJ=1,5)
@@ -54948,7 +54948,7 @@ C-------------------------------------------------------
PP(IDX,5) = XMS
c return reference to other partons
IREFout = LPID(IDX,3)
- IF(NDEBUG.gt.6)
+ IF(IDEBUG.gt.6)
& WRITE(LUN,*) ' final 4momentum: ',(PP(IDX,JJ),JJ=1,5)
END
@@ -54962,14 +54962,14 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
- IF(NDEBUG.gt.6)THEN
+ IF(IDEBUG.gt.6)THEN
WRITE(LUN,*) ' RD_PRTN: (#,PID,LEVEL,REF)',
& IDX,LPID(IDX,1),LPID(IDX,2),LPID(IDX,3)
WRITE(LUN,*) ' 4momentum: ',(PP(IDX,JJ),JJ=1,5)
@@ -54994,8 +54994,8 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -55027,7 +55027,7 @@ C--------------------------------------------------------------------
IFL = LPID(IDX,1)
c return reference to other partons
IREFout = LPID(IDX,3)
- IF(NDEBUG.gt.6)THEN
+ IF(IDEBUG.gt.6)THEN
WRITE(LUN,*) ' RD_PRTN: (#,PID,LEVEL,REF)',
& IDX,IFL,LPID(IDX,2),IREFout
WRITE(LUN,*) ' 4momentum: ',(PIN(JJ),JJ=1,5)
@@ -55047,15 +55047,15 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
IDX = LVL0IDX(JJ)
- IF(ndebug.gt.6)
+ IF(idebug.gt.6)
& WRITE(LUN,*) ' ITR_LVL0_PRTN: JJ,IDX',JJ,IDX
LID = LPID(IDX,1)
IF(JJ+1.gt.NPP0) THEN
@@ -55076,20 +55076,20 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
SAVE
- IF(NDEBUG.gt.6) WRITE(LUN,*) ' PRTN_STCK: reset .. '
- IF(NDEBUG.gt.6) WRITE(LUN,*) ' old state: NPP,NPP0',NPP,NPP0
+ IF(IDEBUG.gt.6) WRITE(LUN,*) ' PRTN_STCK: reset .. '
+ IF(IDEBUG.gt.6) WRITE(LUN,*) ' old state: NPP,NPP0',NPP,NPP0
NPP = NOLD
NPP0 = N0OLD
- IF(NDEBUG.gt.6) WRITE(LUN,*) ' new state: NPP,NPP0',NPP,NPP0
+ IF(IDEBUG.gt.6) WRITE(LUN,*) ' new state: NPP,NPP0',NPP,NPP0
END
@@ -55101,8 +55101,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -55120,8 +55120,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -55148,8 +55148,8 @@ C-------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
PARAMETER (NPP_max = 1000, NPP0_max = 500)
COMMON /S_PRTNS/ PP(NPP_max,5), LPID(NPP_max,4), LVL0IDX(NPP0_max)
& ,NPP,NPP0
@@ -55447,8 +55447,8 @@ C--------------------------------------------------------------------
DIMENSION FACN(3:10)
COMMON /SIB_FAC/ FACN
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DIMENSION P1(5),P2(5),POUT(5)
SAVE
@@ -55459,7 +55459,7 @@ C--------------------------------------------------------------------
IF(XM2.LT.0)THEN
c virtual particle
POUT(5) = -1.D0
- IF(NDEBUG.gt.6)then
+ IF(IDEBUG.gt.6)then
WRITE(LUN,*)
& ' ADD_4VECS: resulting particle virtual!! (m**2):',XM2
WRITE(LUN,*) ' p**2' , POUT(1)**2+POUT(2)**2+POUT(3)**2
@@ -56009,11 +56009,11 @@ C-----------------------------------------------------------------------
INTEGER IPAR
COMMON /S_CFLAFR/ PAR(NPAR_max), IPAR(NIPAR_max)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
SAVE
- if ( ndebug .gt. 0 ) then
+ if ( idebug .gt. 0 ) then
write(lun,*)' -----------------------------------------'
write(lun,*)' SIBYLL DEC_INI: setting particle decays!'
write(lun,*)' to be used in stand-alone SIBYLL only ! '
@@ -56032,15 +56032,15 @@ c all particles with t<0.3e-10s are considered unstable
c i.e. all the mesons from K0s onwards(K0l is stable)
c----------------------------------------------------------
C K0s stable
- if (ndebug .gt. 0 ) write(lun,*)' making K0s stable..'
+ if (idebug .gt. 0 ) write(lun,*)' making K0s stable..'
IDB(12) = -abs(IDB(12))
C Lambda/Anti-lambda stable
- if (ndebug .gt. 0 ) write(lun,*)' making LAMBDA stable..'
+ if (idebug .gt. 0 ) write(lun,*)' making LAMBDA stable..'
IDB(39) = -abs(IDB(39))
c Sigmas stable
- if (ndebug .gt. 0 ) write(lun,*)' making SIGMAs stable..'
+ if (idebug .gt. 0 ) write(lun,*)' making SIGMAs stable..'
do i=34,36
IDB(i) = -abs(IDB(i))
enddo
@@ -56051,11 +56051,11 @@ cfr in reasonable contex eta is never stable !
cdh initializing the pythia routines is done in corsika/SIBINI
c IF(IPAR(44).eq.1)THEN
c use pythia decay routine
-c if (ndebug .gt. 0 ) write(LUN,*) ' using PYTHIA decay routine...'
+c if (idebug .gt. 0 ) write(LUN,*) ' using PYTHIA decay routine...'
c CALL PYDEC_INI
c endif
- if (ndebug .gt. 0 )
+ if (idebug .gt. 0 )
* write(lun,*)' ------------------------------------------'
end
C=======================================================================
@@ -56085,8 +56085,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -56171,7 +56171,7 @@ C--------------------------------------------------------------------
& 77,78,79,80,81,0,83/
DATA (LRES(I),I=84, 99) /94,95,96,97,98,89,4*0,94,95,96,97,98,99/
- IF(Ndebug.gt.3) THEN
+ IF(Idebug.gt.3) THEN
WRITE(LUN,*)
& ' STRING_FRAG_4FLV: called with ',
& '(E0,IFL1,IFL2,PX1,PY1,PX2,PY2,IVAL)',
@@ -56722,7 +56722,7 @@ C...Reset pT and flavor at ends of the string
C...Final two hadrons
400 IAFL1 = IABS(mod(IFL(JR),100))
IAFL2 = IABS(mod(IFL(3),100))
- IF(NDEBUG.gt.5)
+ IF(IDEBUG.gt.5)
& write(lun,*)'STRING_FRAG: final flavors:', IFL(JR), -IFL(3)
C.. check if flavor combination is allowed..
@@ -57005,7 +57005,7 @@ C...reorder particles along chain (in rank)
ENDDO
ENDIF
- if(Ndebug.gt.3)
+ if(Idebug.gt.3)
& WRITE(LUN,*)' STRING_FRAG_4FLV: NP after fragmentation:',NP
END
@@ -57024,8 +57024,8 @@ C...........................................................
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -57064,7 +57064,7 @@ C--------------------------------------------------------------------
DIMENSION WW(2,2),PTOT(4),PX(3),PY(3),IFL(3),PMQ(3)
SAVE
- if(Ndebug.gt.3) then
+ if(Idebug.gt.3) then
WRITE(LUN,*)
& ' GG_FRAG_4FLV: called with (E0)',
& E0
@@ -57292,7 +57292,7 @@ C...Final two hadrons
NP = I+1
c PAR(24) = PAR24def
- IF(Ndebug.gt.3) then
+ IF(Idebug.gt.3) then
WRITE(LUN,*)' GG_FRAG_4FLV: NP after fragmentation:',NP
ENDIF
RETURN
@@ -57319,8 +57319,8 @@ c external types
INTEGER NW_max
PARAMETER (NW_max = 20)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -57415,7 +57415,7 @@ c internal types
DATA PRES /0.7D0/ ! probability of forming a resonance
DATA NRJECT /0/
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,'(2X,A,1x,I2,1x,I2,/,5(2x,F10.3))')
& 'DIFDEC: (L0,Irec,P0):',L0,Irec,(P0(i),i=1,5)
@@ -57441,7 +57441,7 @@ c XI1=FERMI(DELTAE,EMIN,XSMR)
XI1=MAX((EMIN-DELTAE)/XSMR,0.D0)
XSMR=PAR(131)*EMIN2
XI2=MAX((EMIN2-DELTAE)/XSMR,0.D0)
- if(Ndebug.gt.2)
+ if(Idebug.gt.2)
& WRITE(LUN,'(1x,A29,2(2x,F5.2),2(2x,F8.3))')
& ' DIFDEC: EMIN1,EMIN2,XI1,XI2',
& EMIN,EMIN2,Xi1,Xi2
@@ -57449,7 +57449,7 @@ c XI1=FERMI(DELTAE,EMIN,XSMR)
C... pomeron-hadron scattering (pi0 is used instead of pomeron)
IF ((IPAR(10).gt.0).and.(Irec.gt.0).and.
& (DELTAE.gt.EMIN2.or.S_RNDM(LA).gt.XI2)) THEN
- if(Ndebug.gt.2)
+ if(Idebug.gt.2)
& WRITE(LUN,*)' DIFDEC: central (L0,DELTAE,NP,XI):',
& L0,DELTAE,NP,XI2
N1 = NP+1
@@ -57499,7 +57499,7 @@ c restore main event
C..."string-like" decay
ELSE IF (DELTAE .GT. EMIN .or. S_RNDM(LA).gt.XI1) THEN
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,'(2X,A,3(2x,F8.3))')
& 'DIFDEC: string-like, (DELTAE,E0,central prob.):',
& DELTAE,P0(5),1.D0-XI2
@@ -57532,7 +57532,7 @@ c randomize flavor orientation in string
CALL STRING_FRAG_4FLV
+ (P0(5), IFL1, IFL2, 0.D0,0.D0,0.D0,0.D0,IFBAD,-1)
IF (IFBAD .EQ. 1)then
- if(ndebug.gt.1)
+ if(idebug.gt.1)
& WRITE(lun,*)' SIB_DIFF: string-frag rejection! ',
& '(M,NCALL)',P0(5),NCALL
NTRYS = NTRYS + 1
@@ -57561,7 +57561,7 @@ c randomize flavor orientation in string
C...Phase space decay of the excited state
ELSEIF(DELTAE.GT.AM(7)+0.02D0)THEN
- if(Ndebug.gt.2)
+ if(Idebug.gt.2)
& WRITE(LUN,*)' DIFDEC: fireball, (DELTAE,string prob.):',
& DELTAE,1.D0-XI1
IF(IPAR(14).GT.0.and.IPAR(14).NE.7)THEN
@@ -57725,14 +57725,14 @@ c LL(K) = 27-INT(S_RNDM(0)*2.99999D0)
ELSEIF (IPAR(14).LE.-1) THEN
C... generalized fireball model
- IF(Ndebug.gt.2)
+ IF(Idebug.gt.2)
& WRITE(LUN,*)' DIFDEC: using generalized fireball!'
c set charge exchange probability,
c i.e. prob for p* -> n + pip
PCXG = PAR(61)
CALL FIREBALL_4FLV(L0,P0,PCXG,IFBAD)
IF(IFBAD.eq.1)THEN
- IF(ndebug.gt.0)THEN
+ IF(idebug.gt.0)THEN
IF(NRJECT.le.10)THEN
WRITE(LUN,*)
& ' DIFDEC: warning: fireball rejection! ',
@@ -57789,7 +57789,7 @@ c print *,'npi:',npi,'av',av,'p05',p0(5),am(la),deltae
ENDDO
ENDIF
ELSE
- if (ndebug .gt. 0) then
+ if (idebug .gt. 0) then
IF(NRJECT.le.10)THEN
WRITE(LUN,*) ' DIFDEC rejection! ',
& 'diff. mass to low to dissociate beam!'
@@ -57817,8 +57817,8 @@ C routine to produce massive excitations of beam and/or target \FR'14
C-----------------------------------------------------------------------
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION AM,AM2
COMMON /S_MASS1/ AM(99), AM2(99)
INTEGER NW_max
@@ -57927,11 +57927,11 @@ c DOUBLE PRECISION XDUMMY
c default return point, beam and target sampling
c IREJ = 1
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*) ' EXCT_RMNT: input (JW,KRMNT,IREJ)',
& JW,KRMNT,IREJ
- IF(NDEBUG.gt.3)THEN
+ IF(IDEBUG.gt.3)THEN
write(LUN,*)' beam remnant index: (lvl0,flv1,flv2) ',IBMRDX
write(LUN,*)' 1st central string index: (lvl0,bm,tg)',
& (ICSTDX(2*(JW-1)+1,ii),ii=1,3)
@@ -58001,7 +58001,7 @@ c no excited remnant case, jump straight to central strings..
ENDIF
- IF(NDEBUG.gt.3)then
+ IF(IDEBUG.gt.3)then
write(lun,*) ' beam parton1: ',IBM1
write(lun,*) ' beam parton2: ',IBM2
write(lun,*) ' target parton1:',ITG1
@@ -58013,7 +58013,7 @@ c save status of parton stack
10 ITRY(5) = ITRY(5) + 1
IF(ITRY(5).GT.NREJ(5))THEN
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*) ' EXCT_RMNT: no. of trials exceeded, ',
& NREJ(5), 'resample minijets ...' , IREJ
RETURN
@@ -58053,7 +58053,7 @@ c add beam and target side to get total 4momentum
ECM = PTT(5)
c catch virtual remnants
IF(PTT(5).LT.0.D0) THEN
- IF(NDEBUG.GT.2)THEN
+ IF(IDEBUG.GT.2)THEN
WRITE(LUN,*) ' EXCT_RMNT: too little mass left (Shat):',
& SHAT
WRITE(LUN,*) ' resample minijets...'
@@ -58063,8 +58063,8 @@ c catch virtual remnants
ENDIF
- IF(NDEBUG.GT.2) WRITE(LUN,*) ' EXCT_RMNT: try no.',ITRY(5)
- IF(NDEBUG.GT.3)THEN
+ IF(IDEBUG.GT.2) WRITE(LUN,*) ' EXCT_RMNT: try no.',ITRY(5)
+ IF(IDEBUG.GT.3)THEN
write(LUN,*) ' 4momenta before scattering:'
write(LUN,*) ' PBM1:' , (PBM1(jj),jj=1,5)
write(LUN,*) ' PBM2:' , (PBM2(jj),jj=1,5)
@@ -58077,7 +58077,7 @@ c catch virtual remnants
write(LUN,*) ' PTT:' , (PTT(jj),jj=1,5)
ENDIF
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: SHAT:',SHAT
XMFRAC = PAR(81)
@@ -58113,7 +58113,7 @@ C select exponent from COMMON
ALPHA = XRMEX(LL(IABS(KB)))
c sample beam mass
XMB2 = XM2DIS(XMSQMIN,XM2MAX,ALPHA)
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: XM2min,XM2max,ALPHA,XM2:',
& XMSQMIN,XM2MAX,ALPHA,XMB2
c check if resonance or massive hadron has to be formed
@@ -58131,7 +58131,7 @@ c target side mass
C select exponent from COMMON
ALPHA = XRMEX(LL(IABS(KT(JW))))
XMT2 = XM2DIS(XMSQMIN,XM2MAX,ALPHA)
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: XM2min,XM2max,ALPHA,XM2:',
& XMSQMIN,XM2MAX,ALPHA,XMT2
@@ -58154,7 +58154,7 @@ c remnant mass only increased by multiple interactions..
C select exponent from COMMON
ALPHA = XRMEX(LL(IABS(KB)))
XMB2 = XM2DIS(XMSQMIN,XM2MAX,ALPHA)
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: XM2min,XM2max,ALPHA,xm2:',
& XMSQMIN,XM2MAX,ALPHA,XMB2
@@ -58170,7 +58170,7 @@ c target always nucleon
C select exponent from COMMON
ALPHA = XRMEX(LL(IABS(KT(JW))))
XMT2 = XM2DIS(XMSQMIN,XM2MAX,ALPHA)
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: XM2min,XM2max,alpha,XM2:',
& XMSQMIN,XM2MAX,ALPHA,XMT2
@@ -58190,7 +58190,7 @@ c IF(SHAT.lt.XMB2+XMT2+0.3) GOTO 10
C transfer cm energy to mass of particle in parton-parton cm
CALL TRANSFONSHELL(ECM,XMB,XMT,XM2MAX,1,P1,P2,IBD)
IF(IBD.eq.1) THEN
- IF(NDEBUG.gt.2) WRITE(LUN,*) ' EXCT_RMNT: excitation rejected!'
+ IF(IDEBUG.gt.2) WRITE(LUN,*) ' EXCT_RMNT: excitation rejected!'
RETURN
ENDIF
@@ -58256,7 +58256,7 @@ c assign transverse momentum of partons according to random fraction
ENDDO
ENDIF
- IF(NDEBUG.GT.3)THEN
+ IF(IDEBUG.GT.3)THEN
write(LUN,*) ' parton 4momenta after scattering:'
write(LUN,*) ' PBM1:' , (PBM1(jj),jj=1,5)
write(LUN,*) ' PBM2:' , (PBM2(jj),jj=1,5)
@@ -58421,14 +58421,14 @@ c m2_max should be smaller than m2_min
c print *,'jw,jdif,nwd,l0,ifl1,ifl2,deltae,xmin,ee,xmax',
c & jw,jdif(jw),nwd,l0,ifl1,ifl2,deltae,xmmin,ee,par(13)*ee2
IF(DELTAE.lt.AM(7)+0.02D0) THEN
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& WRITE(lun,*) ' EXCT_RMNT: inchoherent diff. :',
& ' not enough mass left for excitation! (DELTAE,PION,',
& 'IREJ,NCALL)',DELTAE,AM(7)+0.02D0,IREJ,NCALL
RETURN
ENDIF
IF(PAR(13)*EE2.lt.XMMIN)THEN
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& WRITE(lun,*) ' EXCT_RMNT: inchoherent diff. :',
& ' not enough mass left for excitation! (min,max,',
& 'IREJ,NCALL)',PAR(13)*EE2,XMMIN,IREJ,NCALL
@@ -58441,7 +58441,7 @@ c reset references of partons
CALL ADD_REF(IBMST2,ITGST1)
CALL ADD_REF(ITGST1,ITGST2)
CALL ADD_REF(ITGST2,IREF)
- IF(ndebug.gt.2) THEN
+ IF(idebug.gt.2) THEN
WRITE(LUN,*) ' EXCT_RMNT: incoherent diff. ',
& '(IDX,IDX2,JDIF,ECM,L0)',IREF,IDHAD,JDIF(JW),PTT(5),L0
WRITE(LUN,*) ' EXCT_RMNT: DELTAE,XM2MAX:',DELTAE,PAR(13)*EE2
@@ -58475,9 +58475,9 @@ c available mass for string
c catch virtual strings
IF(PTT(5).lt.0.D0) THEN
IREJ = 1
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(LUN,*)' EXCT_RMNT: virt. string (M):',EE
- IF(ndebug.gt.3)then
+ IF(idebug.gt.3)then
CALL GET_IMASS2(IBMST,XM2)
write(LUN,*) ' PBM1:', (PBM1(j),j=1,5),XM2
CALL GET_IMASS2(ITGST,xm2)
@@ -58508,14 +58508,14 @@ c skip if two charm quarks need merging..
IF(EE.gt.AM(IABS(L0))) then
IMRG = IMRG + JJ
CALL ADD_PRTN_4VEC(PTT,L0,2,IBMST,ISTH)
- IF(ndebug.gt.2)then
+ IF(idebug.gt.2)then
write(lun,*)
& ' EXCT_RMNT: c.string mass too low! ',
& 'merge into hadron..',l0
ENDIF
ENDIF
ELSE
- IF(ndebug.gt.2)then
+ IF(idebug.gt.2)then
write(lun,*)
& ' EXCT_RMNT: c.string kinematic rejection!'
write(lun,*) ' EE,limit,XMT1,XMT2:',
@@ -58546,7 +58546,7 @@ c add reference to corresponding central string to target parton
c form single hadron from string if mass was too low ..
c need to put hadron on shell by exchanging energy with other string
IF(IMRG.eq.1.or.IMRG.eq.2)THEN
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: merging one string..',IMRG
c one string merged
c index of merged string and its last parton
@@ -58562,7 +58562,7 @@ c string two
CALL RD_PRTN_4VEC(ICST2,P2,IFL2,IDM)
c cm energy
CALL ADD_4VECS(P1,P2,PTT)
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
write(lun,*)' EXCT_RMNT: string A :',(P1(i),i=1,5)
write(lun,*)' EXCT_RMNT: string B :',(P2(i),i=1,5)
write(lun,*)' EXCT_RMNT: total :',(PTT(i),i=1,5)
@@ -58572,7 +58572,7 @@ c cm energy
XM2 = P2(5)
CALL TRANSFONSHELL(ECM,XM1,XM2,1.D0,3,P1N,P2N,LBD)
IF(LBD.eq.1) THEN
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: mass transfer failed!'
RETURN
ENDIF
@@ -58600,7 +58600,7 @@ c rotation factors
COF=COF/ANORF
SIF=SIF/ANORF
ENDIF
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
write(lun,*)' EXCT_RMNT: momentum in cm:',(PBM1(i),i=1,5)
write(lun,*)' EXCT_RMNT: rotation factors:',COD,SID,COF,SIF
write(lun,*)' EXCT_RMNT: rotation angles (theta,phi):',
@@ -58619,7 +58619,7 @@ c rotate parton momenta after interaction, still in parton-parton frame
P2N(1)=PX
P2N(2)=PY
P2N(3)=PZ
- IF(ndebug.gt.2) write(lun,*)' EXCT_RMNT: momentum*:',
+ IF(idebug.gt.2) write(lun,*)' EXCT_RMNT: momentum*:',
& sqrt(P1N(1)**2+P1N(2)**2+P1N(3)**2)
c boost back to hadron-hadron
@@ -58634,7 +58634,7 @@ c boost back to hadron-hadron
& P2N(1),P2N(2),P2N(3),P2N(4),
& P2TOT,P2(1),P2(2),P2(3),P2(4))
p2(5)=p2n(5)
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
write(lun,*)' EXCT_RMNT: momenta after scattering:'
write(lun,*)' EXCT_RMNT: hadron A :',(P1(i),i=1,5)
write(lun,*)' EXCT_RMNT: string B :',(P2(i),i=1,5)
@@ -58651,7 +58651,7 @@ c edit partons on stack
& (ICST2,P2(1),P2(2),P2(3),P2(4),P2(5),ICST21)
ELSEIF(IMRG.eq.3)THEN
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: merge both strings..'
c both strings merged
@@ -58673,7 +58673,7 @@ c cm energy
CALL ADD_4VECS(P1,P2,PTT)
ECM = PTT(5)
ETOT = PTT(4)
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
write(lun,*)' EXCT_RMNT: string A :',(P1(i),i=1,5)
write(lun,*)' EXCT_RMNT: string B :',(P2(i),i=1,5)
write(lun,*)' EXCT_RMNT: total :',(PTT(i),i=1,5)
@@ -58681,7 +58681,7 @@ c cm energy
CALL TRANSFONSHELL(ecm,xm1,xm2,1.D0,3,P1n,P2n,LBD)
IF(LBD.eq.1) THEN
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*)' EXCT_RMNT: mass transfer failed!'
RETURN
ENDIF
@@ -58728,7 +58728,7 @@ c boost massive hadrons back to hadron-hadron
& P2N(1),P2N(2),P2N(3),P2N(4),
& P2TOT,P2(1),P2(2),P2(3),P2(4))
P2(5)=P2N(5)
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
write(lun,*)' EXCT_RMNT: hadron A :',(P1(i),i=1,5)
write(lun,*)' EXCT_RMNT: hadron B :',(P2(i),i=1,5)
ENDIF
@@ -58764,8 +58764,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -58899,7 +58899,7 @@ c xmpsuppr = prob. accepting additional proton
ENDIF
XTEMPH=(AM(6)-AM(13))/dLOG(XMPSUPPR)
- IF(Ndebug.gt.3) THEN
+ IF(Idebug.gt.3) THEN
WRITE(LUN,*)' FIRBALL_4FLV: called with (L0,P0):',
& L0,P0
WRITE(LUN,*)' 2nd Proton rejection prob.:',XMPSUPPR
@@ -58929,7 +58929,7 @@ c RETURN
DELTAE = P0(5) - AM(LA)
IF(DELTAE.lt.AM(6)+0.02D0)THEN
IREJ = 1
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& WRITE(LUN,*)' FIRBALL_4FLV: too low mass!! aborting...',IREJ
c xa=-1.
c xa=log(xa)
@@ -58943,14 +58943,14 @@ c at least two
200 XRNDM = GASDEV(LA)
NPI = INT(AV*(1.D0+PAR(38)*XRNDM))
XMMIN = AM(LA)+DBLE(NPI-1)*AM(6)+0.02D0
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
& WRITE(LUN,*)' NPI,av,rndm,xmin,delta',
& NPI,av,XRNDM,xmmin,P0(5)-XMMIN
IF((NPI.LE.1).OR.(NPI.GT.9).OR.(P0(5).LT.XMMIN))THEN
GOTO 200
ENDIF
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
& WRITE(LUN,*)' FIRBALL_4FLV: No. of particles sampled. ',
& '(NPI,DELTAE,NTRY):',NPI,DELTAE,NTRY
@@ -58965,7 +58965,7 @@ cdh IF(NTRY.GT.5) GOTO 211
c split last hadron again to start hadron chain
CALL HSPLI (LL(I+1),IFL(1),IFL(2))
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
& WRITE(LUN,*)' FIRBALL_4FLV: Input hadron split. ',
& '(L0,IFL1,IFL2):',LL(I+1),IFL(1),IFL(2)
WREM = P0(5)
@@ -59002,7 +59002,7 @@ c add strangeness
& DELTAE.gt.XLIMLAM)THEN
LA1 = LSTR(LA)
c print *,'xlim<deltae?: ',xlimlam,deltae
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
&write(lun,*)' FIRBALL_4FLV: producing hyperon:',namp(LA),namp(LA1)
endif
IF(PCHEX.gt.S_RNDM(LA1))THEN
@@ -59032,7 +59032,7 @@ c split last hadron again to start hadron chain
cdh 211 CALL HSPLI (LL(I+1),IFL(1),IFL(2))
CALL HSPLI (LL(I+1),IFL(1),IFL(2))
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
& WRITE(LUN,*)' FIRBALL_4FLV: Input hadron split. ',
& '(L0,IFL1,IFL2):',LL(I+1),IFL(1),IFL(2)
WREM = P0(5)
@@ -59054,7 +59054,7 @@ cdh 211 CALL HSPLI (LL(I+1),IFL(1),IFL(2))
IF(NTRYS.gt.50) GOTO 210
NTRYS=NTRYS+1
W=dEXP(-AM(ABS(LL(I)))/XTEMPH)
- IF(Ndebug.gt.4)
+ IF(Idebug.gt.4)
& WRITE(LUN,*)' FIRBALL_4FLV: flavor added: ',
& '(I,NTRYS,LL(I),IFL3,W):',I,NTRYS,LL(I),IFL(3),W
IF(W.LT.S_RNDM(I).and.INONLEAD(JT).eq.1) GOTO 240
@@ -59062,12 +59062,12 @@ cdh 211 CALL HSPLI (LL(I+1),IFL(1),IFL(2))
c... kinematic limits...
WREM = WREM-AM(IABS(LL(I)))
WREM2_2=WREM2+2.D0*dSQRT(WREM2)*AM(IABS(LL(I)))+AM2(IABS(LL(I)))
- IF(Ndebug.gt.4)
+ IF(Idebug.gt.4)
& WRITE(LUN,*)' FIRBALL_4FLV: kinematic limits: ',
& '(I,NTRYS,P05**2,WREM2):',I,NTRYS,P0(5)**2,WREM2_2
IF(WREM2_2+0.2D0*S_RNDM(I+1).ge.P0(5)**2) GOTO 240
WREM2=WREM2_2
- IF(Ndebug.gt.3)
+ IF(Idebug.gt.3)
& WRITE(LUN,*)
& ' FIRBALL_4FLV: Hadron added: (KF,NAMP,I,NONlead,WRME2)',
& LL(I),NAMP(ABS(LL(I))),I,INONLEAD(JT),WREM2
@@ -59093,14 +59093,14 @@ c$$$ IF(IAFL1*IAFL2.GT.100) GOTO 100
IF(NTRYC.gt.10) GOTO 210
NTRYC=NTRYC+1
WREM2_2=WREM2+2.D0*dSQRT(WREM2)*AM(ABS(LL(I)))+AM2(ABS(LL(I)))
- IF(Ndebug.gt.5)
+ IF(Idebug.gt.5)
& WRITE(LUN,*)' FIRBALL_4FLV: closing List: (IFL1,IFL2,KF,',
& 'NAMP,I,NTRYC,WREM2)',
& IFL(JT),IFL(JR),LL(I),NAMP(ABS(LL(I))),I,NTRYC,WREM2_2
IF(WREM2_2+0.2D0*S_RNDM(I).ge.P0(5)**2) GOTO 250
- 300 IF(Ndebug.gt.3)
+ 300 IF(Idebug.gt.3)
& WRITE(LUN,*)
& ' FIRBALL_4FLV: flavors sampled. (NPI,LL,WREM,NTRYL):',
& NPI,(LL(ii),ii=1,NPI),WREM,NTRYL
@@ -59148,8 +59148,8 @@ c external types
DOUBLE PRECISION SQS,SIGT,SIGEL,SIGINEL,SLOPE,RHO!,SIGDIF
integer l,kt
c commons
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -59228,7 +59228,7 @@ c write(LUN,*) ' l,k,sig,rho:',l,k,sig,rho
c derive elastic and inelastic cross section
sigel = sigt**2*(1.D0+rho**2)/(16.D0*PI*slope*cmbarn)
siginel = sigt-sigel
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(LUN,*)
& ' SIG_RPP2014: L,KT,SQS,SIGT,SIGEL,SIGINEL,SLOPE,RHO',
& L,KT,SQS,SIGT,SIGEL,SIGINEL,SLOPE,RHO
@@ -59281,8 +59281,8 @@ C--------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -59329,10 +59329,10 @@ c parton stack index of incoming hadron
IKDH = 0
c if too low, fallback on beam
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*)' SEL_RES: input (XM2in,KDin,IRDX):',XM2,KD,IRDX
DELTAE = XM1-AM(ABS(KD))
- IF(ndebug.gt.1)then
+ IF(idebug.gt.1)then
write(lun,*)' SEL_RES: DELTAE,EMIN1,EMIN2',deltae,emin1,emin2
write(lun,*)' SEL_RES: XM,XM1,XM2',
& XM1,emin1+AM(ABS(KD)),emin2+AM(ABS(KD))
@@ -59356,7 +59356,7 @@ cdh & +INT((2.D0-EPS10)*S_RNDM(KD))
IF(KDD.eq.11.or.KDD.eq.12)KDD=21
& +INT(0.5D0+S_RNDM(KD))
IL = MRES(KDD,II)
- IF(ndebug.gt.2) then
+ IF(idebug.gt.2) then
write(lun,*) ' SEL_RES: res. select (KD,II,IL):',
& KD,II,IL
ENDif
@@ -59368,7 +59368,7 @@ c from the relativistic breit-wigner dist.
c scale widths to artificially increase or decrease resonance occurence
XWDTH = PAR(94)*AW2(IL)
PRES = BREIT_WIGNER(XM2,AM2(IL),XWDTH)
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*)
& ' SEL_RES: res. proposal (AM2,AW2,Prob.):',
& AM2(IL),XWDTH,PRES
@@ -59402,7 +59402,7 @@ c add new beam hadron to stack
CALL ADD_PRTN
& (0.D0,0.D0,0.D0,0.D0,XM1,KDH,2,IRDX,IKDH)
endif
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*)' SEL_RES: output (XM2in,KDin,KDH):',XM2,KD,KDH
RETURN
@@ -59450,8 +59450,8 @@ C-----------------------------------------------------------------------
IMPLICIT NONE
INTEGER IDX,IBAD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -59548,9 +59548,9 @@ c FI = ASIN(MIN((1.D0-EPS8),PBM(2)/PTJET))
c TH = PST(1)
c FI = PST(2)
- IF(NDEBUG.gt.1) WRITE(LUN,*)' FRAG_MINIJET: IDX,EE,IFLB,IFLT,PT',
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' FRAG_MINIJET: IDX,EE,IFLB,IFLT,PT',
& IDX,E0,IFLB,IFLT,PTJET,IBAD
- IF(NDEBUG.gt.1) WRITE(LUN,*)' FRAG_MINIJET: PTJET,TH,FI:',
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' FRAG_MINIJET: PTJET,TH,FI:',
& PTJET,TH,FI
C... parameter setup (string fragmentation)
@@ -59656,7 +59656,7 @@ c leading charm fraction
CALL STRING_FRAG_4FLV
& (E0,IFL1,-IFL1,0.D0,0.D0,0.D0,0.D0,IFBAD,0)
if(IFBAD.gt.0) then
- IF(ndebug.gt.1)
+ IF(idebug.gt.1)
& WRITE(LUN,*)
& ' JET_FRAG: rejection in STRING_FRAG (IFL,E0,NCALL):',
& IFL1,E0,NCALL
@@ -59678,7 +59678,7 @@ c DBETJ = (DX1J-DX2J)/(DX1J+DX2J)
DBETJ = PST(3)/PST(4)
CALL SIROBO (NOLD+1,NP,TH,FI,0.D0,0.D0,DBETJ)
- if(Ndebug.gt.1) WRITE(LUN,*)
+ if(Idebug.gt.1) WRITE(LUN,*)
& ' JET_FRAG: particles produced:',NP-NOLD
PAR(24) = PAR24_def
PAR(1) = PAR1_def
@@ -59771,8 +59771,8 @@ C----------------------------------------------------------
IMPLICIT NONE
character*16 text
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER II2,JJ2
DOUBLE PRECISION U2,C2,CD2,CM2
COMMON /SIB_RAND/ U2(97),C2,CD2,CM2,II2,JJ2
@@ -59828,9 +59828,9 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
-c COMMON /S_DEBUG/ Ncall, Ndebug, Lun
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
+c COMMON /S_DEBUG/ Ncall, Idebug, Lun
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -59869,7 +59869,7 @@ C check if tables initialized
xa = log(xa)
stop
ENDIF
- IF(NDEBUG.GT.1)
+ IF(IDEBUG.GT.1)
& WRITE(LUN,*) ' CUT_PRO: input: L, SQS, PTmin',L, SQS, PTmin
c choose nucleon or meson table
@@ -59923,7 +59923,7 @@ C...phase space limitation
NSOFR = I
NJETR = J
- if(Ndebug.gt.1)
+ if(Idebug.gt.1)
& write(lun,*)' CUT_PRO: (L,SQS,PTmin,Ns,Nh) ',K,SQS,PTmin,I,J
END
@@ -59961,8 +59961,8 @@ C-----------------------------------------------------------------------
COMMON /S_CCSIG3/ SSIG_SD1LM(61,3),SSIG_SD1HM(61,3),
& SSIG_SD2LM(61,3),SSIG_SD2HM(61,3),
& SSIG_DDLM(61,3),SSIG_DDHM(61,3)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DIMENSION Pjet(0:NS_max,0:NH_max)
DIMENSION SIG_df(3),SIG_df2(3,2),SIGDIF(3),SIGDIF_pi(3),
@@ -59994,7 +59994,7 @@ C...initialization of proton and pion tables
IF(LUN.ne.6) WRITE(6,*)' Calculating cross section tables...'
DO KK=1,2
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& WRITE(LUN,'(2(/,1X,A,A))')
& 'Table: J, sqs, PT_cut, SIG_tot, SIG_inel, B_el, ',
& 'rho, <n_s>, <n_h>, SIG_SD, SD1_lm, SD1_hm',
@@ -60056,7 +60056,7 @@ C...low-energy interpolation with data-parametrizations
PH_tab(J) = PH
PT_tab(J) = PTmin
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& WRITE(LUN,'(3X,I2,1P,E12.3,0P,4F8.2,6F8.3)')
& JINT,SQS,PTmin,SIG_tot,SIG_inel,B_el,RHO,PS,PH
& ,SIGDIF(1)+SIGDIF(2),SIG_df2(1,1),SIG_df2(1,2)
@@ -60068,7 +60068,7 @@ C...initialization of kaon tables
JINT = 3
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& WRITE(LUN,'(2(/,1X,A,A))')
& 'Table: J, sqs, PT_cut, SIG_tot, SIG_inel, B_el, ',
& 'rho, <n_s>, <n_h>',
@@ -60116,7 +60116,7 @@ C...low-energy interpolation with data-parametrizations
SSIG_B(J,3) = B_EL
SSIG_RHO(J,3) = RHO
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& WRITE(LUN,'(3X,I2,1P,E12.3,0P,4F8.2,3F8.3)')
& JINT,SQS,PTmin,SIG_tot,SIG_inel,B_el,RHO,PS,PH
@@ -60151,8 +60151,8 @@ C-----------------------------------------------------------------------
SAVE
DATA ATARGET /14.514D0/
- if ( ndebug .gt. 3 ) CALL PARAM_PRINT(LUN)
- if ( ndebug .gt. 0 ) THEN
+ if ( idebug .gt. 3 ) CALL PARAM_PRINT(LUN)
+ if ( idebug .gt. 0 ) THEN
WRITE (LUN, 10)
WRITE (LUN, 15)
WRITE (LUN, 16)
@@ -60226,8 +60226,8 @@ C-----------------------------------------------------------------------
& ALINT(61,3,3), ASQSMIN, ASQSMAX, DASQS, NSQS
COMMON /GLAUB_SCR/ XI_MAX , ALAM(61)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -60244,12 +60244,12 @@ C-----------------------------------------------------------------------
DATA TARGN /'air','nit','oxy'/
IF ( IPAR(12).GT.0 ) THEN
- if (ndebug.gt.0) then
+ if (idebug.gt.0) then
WRITE(LUN,*) ' SIG_AIR_INI:'
WRITE(LUN,*)' using Goulianos param. for res.coupling..'
endif
XI_MAX = 0.02D0
- if (ndebug.gt.0)WRITE(LUN,*)' low mass Xi_max: ' , XI_MAX
+ if (idebug.gt.0)WRITE(LUN,*)' low mass Xi_max: ' , XI_MAX
ENDIF
C...target loop (air, N, O)
@@ -60260,7 +60260,7 @@ C...target loop (air, N, O)
C...particle loop (p, pi, K)
DO K=1,3
- if (NDEBUG .gt. 0 ) then
+ if (IDEBUG .gt. 0 ) then
WRITE(6,'(/,1X,A,A)')
& 'Table: J, IK, sqs, SIGtot, SIGprod, SIG_SD,',
&' Lambda '
@@ -60294,7 +60294,7 @@ c fixed O-N mixture
& (IAT, SIGT, SLOPE, RHO, ALAM(J),
& SSIGT, SSIGEL, SSIGQE, SSIGSD, SSIGQSD)
ENDIF
- if (ndebug .gt. 0 ) WRITE(6,'(1X,I2,1P,5E12.3)')
+ if (idebug .gt. 0 ) WRITE(6,'(1X,I2,1P,5E12.3)')
& K,SQS,SSIGT,SSIGT-SSIGQE,SSIGQSD,ALAM(J)
C particle production cross section
SSIGN(J,K,IK) = SSIGT-SSIGQE
@@ -60311,7 +60311,7 @@ c interaction length
ENDDO
ENDDO
- if (ndebug .gt. 0 ) then
+ if (idebug .gt. 0 ) then
WRITE(6,'(/,1X,A)')
& ' SIG_AIR_INI: NUCLIB interaction lengths [g/cm**2]'
WRITE(6,*) 'target:', TARGN(IK)
@@ -60349,8 +60349,8 @@ c external types
DIMENSION KRMNT(NW_max)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -60414,7 +60414,7 @@ c internal types
IREJ1 = 1
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
+ WRITE(LUN,*)
+ ' SAMPLE_TARGET: NW,XCHG,LRMNT,XJET,IREC,IREJ',
+ NW,XCHG,(KRMNT(j),j=1,NW),(XJET(j),j=1,NW),IREC,IREJ
@@ -60520,7 +60520,7 @@ c pz = (0.5D0*SQS*X2(JJ))**2-PXT(JJ)**2-PYT(JJ)**2
enddo
ENDIF
ENDDO
- IF(NDEBUG.GT.3) CALL PRNT_PRTN_STCK
+ IF(IDEBUG.GT.3) CALL PRNT_PRTN_STCK
IREJ = 0
END
@@ -60876,8 +60876,8 @@ C......................................................................
DOUBLE PRECISION SSIG,SLOPE,ALPHA,ALAM
DOUBLE PRECISION SIG1,SIGEL1,SIGQE1,SIGSD1,SIGQSD1
DOUBLE PRECISION SIGT,SIGEL,SIGQE,SIGSD,SIGQSD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
SAVE
IF(IAT.eq.0.or.IAT.gt.18) THEN
WRITE(LUN,'(//,1X,A)')
@@ -61117,8 +61117,8 @@ C
C---------------------------------------------------------
IMPLICIT NONE
c Commons
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -61153,7 +61153,7 @@ C internals
DOUBLE PRECISION S_RNDM
SAVE
- IF(NDEBUG.gt.0)write(lun,*)
+ IF(IDEBUG.gt.0)write(lun,*)
& ' REMOVE_PI0: Rate,Mode:',xrate,IPAR(50)
C select exchange model
IF(IPAR(50).eq.1)THEN
@@ -61167,7 +61167,7 @@ c replace with pi+ or pi-
LL=LL+IFPI0*(2-INT(MIN((2.D0+XRATE)*S_RNDM(LA),
& 3.D0-EPS10)))
LLIST(I) = LL
- IF(NDEBUG.gt.1)
+ IF(IDEBUG.gt.1)
& WRITE(LUN,*) ' REMOVE_PI0: LA,IFPI0,LNEW:',LA,IFPI0,LL
ENDDO
ENDIF
@@ -61180,8 +61180,8 @@ C=======================================================================
C-----------------------------------------------------------------------
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -61216,7 +61216,7 @@ C--------------------------------------------------------------------
SAVE
DATA ICNT2 /0/
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
& write(lun,*)' SAMPLE_SEA_INDV: called with ',
& '(KRMNT,XMINA,XMINA_SEA,NSEA,XREM0,ALPHA,ASUP,XQMASS,XMAX):',
& KRMNT,XMINA,XMINA_SEA,NSEA,XREM0,ALPHA,ASUP,XQMASS,XMAX
@@ -61233,7 +61233,7 @@ C--------------------------------------------------------------------
ENDIF
IF(ITRY(4).gt.Nsea/2*NREJ(4))THEN
ICNT2 = ICNT2 + 1
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
IF(ICNT2.le.5)THEN
write(lun,*)' SAMPLE_SEA_INDV: rejection!'
write(lun,*)' reached max. no. of trials!', NREJ(4)
@@ -61255,12 +61255,12 @@ c scale for interactions other than first if Nw>1
XX(jj1) = x1
XX(jj2) = x2
XREM = XREM - XX(jj1) - XX(jj2)
- IF(NDEBUG.gt.2)
+ IF(IDEBUG.gt.2)
& WRITE(LUN,*) ' x-frac: JW,X3,X4,XREM',
& J,XX(jj1),XX(jj2),XREM
ENDDO
ITRY(4) = ITRY(4) + 1
- IF(NDEBUG.gt.1) WRITE(LUN,*)
+ IF(IDEBUG.gt.1) WRITE(LUN,*)
& ' SAMPLE_SEA_INDV: ISMPL,XREM0,XREM,XMINA,XMINSEA',
& ITRY(4),XREM0,XREM,XMINA,XMINA_SEA
ENDDO
@@ -61274,8 +61274,8 @@ C=======================================================================
C-----------------------------------------------------------------------
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -61454,7 +61454,7 @@ c put back on mass shell
endif
enddo
ENDIF
- if(ndebug.ge.5) CALL SIB_LIST(6)
+ if(idebug.ge.5) CALL SIB_LIST(6)
END
C=======================================================================
@@ -61472,8 +61472,8 @@ C-------------------------------------------------------------------
DOUBLE PRECISION XCHG,XJET
INTEGER KID,NW,KRMNT,IREJ
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -61537,7 +61537,7 @@ c 2: resample non-diff event (Ns,Nh)..
c 3: resample event (Nw,diff,ndiff)..
IREJ = 1
- IF(ndebug.gt.2)
+ IF(idebug.gt.2)
+ WRITE(LUN,*)
+ ' SAMPLE_BEAM: KID,NW,XCHG,KRMNT,XJET,IREJ',
+ KID,NW,XCHG,KRMNT,XJET,IREJ
@@ -61588,7 +61588,7 @@ c add sea quarks
ICSTDX(2*(J-1)+2,2) = Iref
c add parton index to cache
ENDDO
- IF(NDEBUG.GT.3) CALL PRNT_PRTN_STCK
+ IF(IDEBUG.GT.3) CALL PRNT_PRTN_STCK
IREJ = 0
@@ -61605,8 +61605,8 @@ C-----------------------------------------------------------------------
IMPLICIT NONE
INTEGER IDX,LBAD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -61652,7 +61652,7 @@ C kinematic variables
c set diffraction code of system (1:beam,2:target,3:double)
JDIFF = ABS(IPID)/10
- IF(NDEBUG.gt.1) WRITE(LUN,*)' FRAG_INCHRNT_DIFF: IDX,EE,L0',
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' FRAG_INCHRNT_DIFF: IDX,EE,L0',
& IDX,EE,L0
IPFLAG = -1
@@ -61662,11 +61662,11 @@ c set diffraction code of system (1:beam,2:target,3:double)
c diffractive interaction in center-of-mass system of (sea,rmnt)-nuc
CALL SIB_DIFF(L0,JDIFF,EE,0,LXBAD)
IF(LXBAD.ne.0) THEN
- IF(NDEBUG.gt.1)
+ IF(IDEBUG.gt.1)
& WRITE(LUN,*)' FRAG_INCHRNT_DIFF: fragmentation rejection'
RETURN
ENDIF
- IF(NDEBUG.gt.1)
+ IF(IDEBUG.gt.1)
& WRITE(LUN,*)' FRAG_INCHRNT_DIFF: particles before/after :',
& NOLD,NP
@@ -61732,8 +61732,8 @@ C IAT : mass number of target
INTEGER KB,IAT,KT
COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -61770,7 +61770,7 @@ C--------------------------------------------------------------------
COMMON /SIB_FAC/ FACN
SAVE
- if(Ndebug.gt.1) WRITE(LUN,*)
+ if(Idebug.gt.1) WRITE(LUN,*)
& ' SAMPLE_MINIJETS: (L,NW,NNJET,NNSOF):',
& L,NW,(NNJET(ii),ii=1,nw),(NNSOF(ii),ii=1,nw)
@@ -61792,7 +61792,7 @@ C...hard sea-sea interactions
CALL SAMPLE_HARD (L,X1Jj,X2Jj,PTJET)
X1JET = X1JET + X1Jj
X2JET(JW) = X2JET(JW)+X2Jj
- if(Ndebug.gt.2) THEN
+ if(Idebug.gt.2) THEN
WRITE(LUN,*)
& ' SAMPLE_MINIJETS: hard JJ,JW,X1JET,X2JET(JW):',
& JJ,JW,X1JET,X2JET(JW)
@@ -61800,7 +61800,7 @@ C...hard sea-sea interactions
& ' X1,X2,PT:',X1JJ,X2JJ,PTJET
ENDIF
IF ((X2JET(JW).GT.0.9D0).OR.(X1JET.GT.0.9D0)) then
- if(Ndebug.gt.2) WRITE(LUN,*)
+ if(Idebug.gt.2) WRITE(LUN,*)
& ' SAMPLE_MINIJETS: not enough phase space',
& ' (Ncall,Njet,lbad):',Ncall,Njet,lBAD
return
@@ -61859,7 +61859,7 @@ c set up references
c minijet --> gluon1 --> gluon2 --> minijet
CALL ADD_REF(Irefg1,Irefg2)
CALL ADD_REF(Iref,Irefg1)
- IF(Ndebug.gt.2)THEN
+ IF(Idebug.gt.2)THEN
WRITE(LUN,*)
& ' SAMPLE_MINIJETS: soft JJ,JW,X1JET,X2JET(JW):',
& JJ,JW,X1JET,X2JET(JW)
@@ -61867,7 +61867,7 @@ c minijet --> gluon1 --> gluon2 --> minijet
& ' X1,X2,PT:',X1s,X2s,PTSOF
ENDIF
ELSE
- IF(Ndebug.gt.1) WRITE(LUN,*)
+ IF(Idebug.gt.1) WRITE(LUN,*)
& ' SAMPLE_MINIJETS: not enough phase space',
& ' (Ncall,Nsof,lbad):',Ncall,Njet,lBAD
RETURN
@@ -61903,8 +61903,8 @@ c external types
DOUBLE PRECISION SQS,SIGT,SIGEL,SIGINEL,SLOPE,RHO
INTEGER L0
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -62144,7 +62144,7 @@ C
AL = LOG10(SQS)
J1 = INT((AL-1.D0)*10.D0 + 1)
if((j1.lt.1).or.(j1.gt.NSQS)) then
- if (ndebug .gt. 0)
+ if (idebug .gt. 0)
* write (LUN,'(1x,a,i3,1p,e12.3)')
& ' SIB_SIGMA_EXT: energy out of range ',L,sqs
endif
@@ -62191,8 +62191,8 @@ C-------------------------------------------------------------------
IMPLICIT NONE
C include COMMONs
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -62292,7 +62292,7 @@ C.. initialization
KRMNT = MIN(LRMNT,1)
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
+ WRITE(LUN,*)
+ ' SAMPLE_PROJECTILE: KID,KINT,KRMNT,XCHG,XJET,IREJ',
+ KID,KINT,KRMNT,XCHG,XJET,IREJ
@@ -62309,7 +62309,7 @@ c if remnant is not resolved then on pair less is needed
c (valence pair takes role of one sea pair)
Nsea = 2*(KINT-(1-KRMNT))
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
+ WRITE(LUN,*)
+ ' SAMPLE_PROJECTILE: number of partons to sample ',
+ '(tot,val,sea):',Nval+Nsea,Nval,Nsea
@@ -62325,7 +62325,7 @@ c change proton splitting in case no remnant is present
IF(ITRY(3).gt.NREJ(3)) THEN
ICNT1 = ICNT1 + 1
IF(ICNT1.lt.10)THEN
- if (NDEBUG.gt.0) then
+ if (IDEBUG.gt.0) then
WRITE(LUN,*)' SAMPLE_PROJECTILE: trials exceeded! return..'
WRITE(LUN,*)
+ ' KID,KINT,KRMNT,XCHG,XJET,XVAL,IREJ,NCALL',
@@ -62388,7 +62388,7 @@ C.. check if enough energy left to sample all partons
IF (1.D0-XJET.LT.(Nsea*XMINA_SEA+2.D0*XMINA))THEN
ICNT2 = ICNT2 + 1
IF(ICNT2.le.10)THEN
- IF(NDEBUG.gt.3)THEN
+ IF(IDEBUG.gt.3)THEN
write(lun,*)' SAMPLE_PROJECTILE: rejection!'
write(lun,*)' too little energy to sample all partons!'
write(lun,*)' (NW,Ntot,Nval,Nsea,XMIN,XMIN*N',
@@ -62426,14 +62426,14 @@ c sample from 1/x individually then reject if too large
ALPHA = 1.D0
DO WHILE ( XREM .lt. 2*XMINA )
XREM = 1.D0-XJET
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*) ' N,XREM,XMINA,XMIN2,XMAX2,ALPHA',
& Nsea,XREM,XMINA_SEA,XMINA2,XMAX2,ALPHA
DO j=1,Nsea
jj = 2 + j
IF(KRMNT.eq.0) jj = 4+j
XX(jj) = XM2DIS(XMINA2,XMAX2,ALPHA)
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*) ' J,X,XREM',JJ,XX(JJ),XREM
XREM = XREM - XX(jj)
ENDDO
@@ -62483,7 +62483,7 @@ c scale up for mesons
IFL(j3) = INT(1.5D0+S_RNDM(KID))
ENDIF
IFL(j4) = -IFL(j3)
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*) ' flavor: JW,FLV1,FLV2',J,IFL(j3),IFL(j4)
C... sample sea pt
@@ -62526,13 +62526,13 @@ c avg: PAR(132)
c limit parton virtuality
PTS = MAX(PX(j3)**2+PY(j3)**2,EPS10)
IF((XX(j3)**2+XX(J4)**2)/PTS.lt.8.D0*PAR(122)/S) GOTO 33
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*)' pt: JW,PX,PY,pt',J,Px(j3),Py(j3),sqrt(pts)
ENDDO
C... Prepare the valence partons
100 XVAL=1.D0-XSEAJET
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_PROJECTILE: val. fraction remaining:',
& XVAL
@@ -62570,7 +62570,7 @@ C FOR MESONS, SPLIT ENERGY SYMMETRICALLY.
IF (IABS(KID).LT.13.AND.S_RNDM(0).LE.0.5D0)
& XX(IDXVAL)=XVAL-XX(IDXVAL)
XX(IDXVAL+1)=XVAL-XX(IDXVAL)
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_PROJECTILE: val. sampled (x1,x2):',
& XX(IDXVAL),XX(IDXVAL+1)
c for baryons force diq distribution
@@ -62590,7 +62590,7 @@ C... val. quark transverse momentum
CALL PTDIS_4FLV (10,PX(IDXVAL),PY(IDXVAL))
PX(IDXVAL+1) = -PX(IDXVAL)
PY(IDXVAL+1) = -PY(IDXVAL)
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_PROJECTILE: val. pt (px,py):',
& PX(IDXVAL),PY(IDXVAL)
@@ -62662,14 +62662,14 @@ c reject spin1,isospin singlett
ENDDO
enddo
ENDIF
- IF(ndebug.gt.3)THEN
+ IF(idebug.gt.3)THEN
WRITE(LUN,*)' SAMPLE_PROJECTILE: rmnt PID,NTRY: ',KID1,ITRY(3)
WRITE(LUN,*)' SAMPLE_PROJECTILE: output: I,FLV,PX,PY,X,XSUM'
ENDIF
XSUM = XJET
DO j=IDXVAL,2*(KINT+Krmnt)+2*(1-Krmnt)
XSUM = XSUM + XX(j)
- IF(NDEBUG.gt.3) WRITE(LUN,*) j,IFL(j),PX(J),PY(J),XX(j),XSUM
+ IF(IDEBUG.gt.3) WRITE(LUN,*) j,IFL(j),PX(J),PY(J),XX(j),XSUM
ENDDO
IF(ABS(XSUM-1.D0).gt.EPS3) THEN
WRITE(LUN,*)' SAMPLE_PROJECTILE: parton sum incomplete!',
@@ -62692,8 +62692,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION CBR
INTEGER KDEC,LBARP,IDB
COMMON /S_CSYDEC/ CBR(223+16+12+8), KDEC(1338+6*(16+12+8)),
@@ -62803,8 +62803,8 @@ c external types
DIMENSION SIGDIF(3)
INTEGER L0
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NS_max, NH_max
PARAMETER (NS_max = 20, NH_max = 80)
@@ -62861,7 +62861,7 @@ c internal type declarations
AL = LOG10(SQS)
J1 = INT((AL-1.D0)*10.D0 + 1)
if((j1.lt.1).or.(j1.gt.NSQS)) then
- if(ndebug.gt.0)
+ if(idebug.gt.0)
& write (LUN,'(1x,a,i3,1p,e12.3)')
& ' SIB_SIGMA_HP: energy out of range ',L,sqs
endif
@@ -62910,8 +62910,8 @@ c external types
DIMENSION SIGDIF(3,2)
INTEGER L
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NS_max, NH_max
PARAMETER (NS_max = 20, NH_max = 80)
@@ -62945,7 +62945,7 @@ c internal types
AL = dLOG10(SQS)
J1 = INT((AL - 1.D0)*10.D0 + 1)
if((j1.lt.1).or.(j1.gt.NSQS)) then
- if(ndebug.gt.0)write(LUN,'(1x,a,i3,1p,e12.3)')
+ if(idebug.gt.0)write(LUN,'(1x,a,i3,1p,e12.3)')
& ' SIB_SIGMA_HP2: energy out of range ',L,sqs
endif
if((j1.lt.1).or.(j1.ge.NSQS)) then
@@ -62988,8 +62988,8 @@ Cf2py double precision, intent(in) :: SQS
Cf2py double precision, intent(out) :: SIGprod,SIGbdif
IMPLICIT NONE
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NS_max, NH_max
PARAMETER (NS_max = 20, NH_max = 80)
@@ -63018,7 +63018,7 @@ c internal
AL = LOG10(SQS)
J1 = INT((AL - 1.D0)*10.D0 + 1)
if((j1.lt.1).or.(j1.gt.NSQS)) then
- if (ndebug .gt. 0)
+ if (idebug .gt. 0)
& write (LUN,'(1x,a,i3,1p,e12.3)')
& ' SIB_SIGMA_HAIR: energy out of range ',L,sqs
endif
@@ -63056,8 +63056,8 @@ Cf2py double precision, intent(out) :: SIGprod,SIGbdif,SIGela
INTEGER NS_max, NH_max
PARAMETER (NS_max = 20, NH_max = 80)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION SSIG,PJETC,SSIGN,SSIGNSD,SSIGNEL,ALINT,ASQSMIN,
& ASQSMAX,DASQS
INTEGER NSQS
@@ -63108,7 +63108,7 @@ c read cross section from table
AL = LOG10(SQS)
J1 = INT((AL - 1.D0)*10.D0 + 1)
if((j1.lt.1).or.(j1.gt.NSQS)) then
- if (ndebug .gt. 0)
+ if (idebug .gt. 0)
& write (LUN,'(1x,a,i3,1p,e12.3)')
& ' SIB_SIGMA_HNUC: energy out of range ',L,sqs
endif
@@ -63122,7 +63122,7 @@ c read cross section from table
SIGela = SSIGNEL(J1,L,IK)*(1.D0-T) + SSIGNEL(J1+1,L,IK)*T
ELSEIF(IAT.lt.19)THEN
c calculate cross section
- IF(ndebug.gt.0)THEN
+ IF(idebug.gt.0)THEN
WRITE(LUN,'(1X,A,2I4,F8.2)')
& 'SIB_SIGMA_HNUC: L,IAT,SQS:',L,IAT,SQS
ENDIF
@@ -63140,7 +63140,7 @@ C quasi elastic + elastic singl. diff cross section
SIGbdif = SIGQSD
c elastic cross section
SIGela = SIGel
- if(ndebug.gt.0)THEN
+ if(idebug.gt.0)THEN
WRITE(LUN,'(1X,A,3F8.2)')
& 'SIB_SIGMA_HNUC: SIGprod, SIGbdif, ALAM:',
& SIGprod, SIGbdif, ALAM
@@ -63199,8 +63199,8 @@ C IAT : mass number of target
DOUBLE PRECISION SQS,S,PTmin,XMIN,ZMIN
INTEGER KB,IAT,KT
COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -63218,7 +63218,7 @@ c internal types
DOUBLE PRECISION Z1,Z2,SIG,S_RNDM,Q2,ZSAMPLE
SAVE
- IF(ndebug.gt.2)then
+ IF(idebug.gt.2)then
write(lun,*) ' SAMPLE_HARD: (SQS,S,PTmin,Xmin,Zmin)',
& SQS,S,PTmin,Xmin,Zmin
endif
@@ -63234,7 +63234,7 @@ c internal types
IF(Q2.gt.S*X1*X2) goto 100
PT=dSQRT(Q2*(1.D0-Q2/(S*X1*X2)))
- IF(ndebug.gt.2)then
+ IF(idebug.gt.2)then
write(lun,*) ' SAMPLE_HARD: (X1,X2,PT)',
& X1,X2,PT
endif
@@ -63605,8 +63605,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -63619,21 +63619,21 @@ C-----------------------------------------------------------------------
SAVE
IF(IPAR(8).eq.0) THEN
- if (ndebug .gt. 0 ) WRITE(LUN,*)
+ if (idebug .gt. 0 ) WRITE(LUN,*)
* ' PDF_INI: calcuLating pdf table using Eichten param..'
CALL ZSAMPLE_INI
ELSEIF(IPAR(8).eq.2) THEN
- if (ndebug .gt. 0 ) then
+ if (idebug .gt. 0 ) then
WRITE(LUN,*)' PDF_INI: calculating pdf table using GRV',
* ' param..'
WRITE(LUN,*)' does not work with -fbounds-check !!'
endif
CALL GRV_INI
ELSE
- if (ndebug .gt. 0 ) WRITE(LUN,*)
+ if (idebug .gt. 0 ) WRITE(LUN,*)
* ' PDF_INI: using common table of GRV parametrization..'
ENDIF
- if (ndebug .gt. 0 ) THEN
+ if (idebug .gt. 0 ) THEN
WRITE(LUN,*)APART(1),FFA(1),FFB(1),DX(1),DZ(1)
WRITE(LUN,*)APART(2),FFA(2),FFB(2),DX(2),DZ(2)
ENDIF
@@ -64062,8 +64062,8 @@ C IAT : mass number of target
INTEGER KB,IAT,KT
COMMON /S_RUN/ SQS, S, PTmin, XMIN, ZMIN, KB, KT(NW_max), IAT
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION CCHIK
COMMON /S_CPSPL/ CCHIK(4,99)
@@ -64189,8 +64189,8 @@ C-----------------------------------------------------------------------
c external variables
INTEGER IDX,LBAD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C The final particle output is contained in COMMON /S_PLIST/
C NP : number of final particles
C P(1:NP, 1:5) : 4-momenta + masses of the final particles
@@ -64253,7 +64253,7 @@ c initial parameters
PAR2_def = PAR(2) ! ud/s rate
PAR8_def = PAR(8) ! popcorn rate
PAR24_def = PAR(24) ! c/s rate
- if(ndebug.gt.1)
+ if(idebug.gt.1)
& WRITE(LUN,*) ' EXCTDEC: IDX,IREJ',IDX,LBAD
c read remnant 4momentum from stack
@@ -64266,7 +64266,7 @@ c read remnant 4momentum from stack
write(lun,*) ' EXCTDEC: reference loop broken!',IDX,IR
CALL SIB_REJECT('EXCTDEC ')
endif
- IF(NDEBUG.GT.2)THEN
+ IF(IDEBUG.GT.2)THEN
WRITE(LUN,*) ' EXCTDEC: P0:' , (P0(kk),kk=1,5)
WRITE(LUN,*) ' EXCTDEC: PR1:' , (PR1(kk),kk=1,5)
WRITE(LUN,*) ' EXCTDEC: PR2:' , (PR2(kk),kk=1,5)
@@ -64290,7 +64290,7 @@ c phasespace decay region: EMIN2 < DELTAE < EMIN3
c string decay region: EMIN3 < DELTAE < EMIN4
EMIN4 = PAR(79)
- IF(NDEBUG.gt.2)THEN
+ IF(IDEBUG.gt.2)THEN
WRITE(LUN,*)
& ' EXCTDEC: MASS,IFL1,IFL2,PID',P0(5),IFLR1,IFLR2,KD
WRITE(LUN,*) ' EXCTDEC: DELTAE,EMIN1,EMIN2,EMIN3',
@@ -64314,7 +64314,7 @@ c popcorn rate in remnant
IF(DELTAE.lt.EMIN2)THEN
c beam or resonance region
- IF(NDEBUG.gt.1) then
+ IF(IDEBUG.gt.1) then
if(DELTAE.lt.EMIN1)then
WRITE(LUN,*)' EXCTDEC: fallback to beam..'
else
@@ -64337,7 +64337,7 @@ c beam or resonance region
ELSEIF(DELTAE.lt.EMIN3)THEN
c phasespace decay region
- IF(NDEBUG.gt.1) WRITE(LUN,*)' EXCTDEC: phasespace decay ..'
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' EXCTDEC: phasespace decay ..'
IPFLAG = IPID/iabs(IPID) + ISIGN(1000,IPID)
c set charge exchange probability,
c i.e. prob for p* -> n + pip
@@ -64347,7 +64347,7 @@ c i.e. prob for p* -> n + pip
PAR(8) = PAR8_def
PAR(24) = PAR24_def
IF(IFAIL.eq.1) THEN
- IF(ndebug.gt.0)
+ IF(idebug.gt.0)
& WRITE(LUN,*) ' EXCTDEC: remnant frag. rejection!'
LBAD = 1
RETURN
@@ -64358,7 +64358,7 @@ c i.e. prob for p* -> n + pip
c ELSEIF(DELTAE.lt.EMIN4)THEN
ELSE
C string fragmentation region
- IF(NDEBUG.gt.1) WRITE(LUN,*)' EXCTDEC: string decay ..'
+ IF(IDEBUG.gt.1) WRITE(LUN,*)' EXCTDEC: string decay ..'
N1 = NP+1
IPFLAG = IPFLAG + ISIGN(3000,IPID)
c for meson remnant quark and anti-quark should be treated equally
@@ -64375,7 +64375,7 @@ c turn remnant string around
CALL STRING_FRAG_4FLV
+ (P0(5), IFLR2, IFLR1, 0.D0,0.D0,0.D0,0.D0,IFBAD,1)
IF (IFBAD .EQ. 1)THEN
- IF(ndebug.gt.0)
+ IF(idebug.gt.0)
& WRITE(LUN,*) ' EXCTDEC: remnant frag. rejection!'
LBAD = 1
PAR(2) = PAR2_def
@@ -64409,7 +64409,7 @@ c SID= 2.D0*PTR/P0(5)*COD
COF=COF/ANORF
SIF=SIF/ANORF
ENDIF
- IF(ndebug.gt.3)THEN
+ IF(idebug.gt.3)THEN
write(lun,*)' EXCTDEC: rotation factors (cod,sid,cof,sif):',
& cod,sid,cof,sif
write(lun,*)' EXCTDEC: rotation angles (theta,phi):',
@@ -64424,7 +64424,7 @@ c rotate string final state
enddo
ENDDO
c boost to hadron-hadron center-of-mass
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' EXCTDEC: boost to had-had (gabe,gam):',
& (gabe(j),j=1,4)
DO K=N1,NP
@@ -64535,8 +64535,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -64637,7 +64637,7 @@ c rc4a charm pt
PPT02(J) = PTCHB
ENDDO
- IF(ndebug.gt.2)THEN
+ IF(idebug.gt.2)THEN
WRITE(LUN,*)' PTSETUP_4FLV: (sqs,(u,d),s,diq,pom,cm,cb)',sqs
+ ,ppt02(1),ppt02(3),ppt02(11), ppt02(10),ppt02(4),ppt02(34)
ENDIF
@@ -64682,8 +64682,8 @@ C-----------------------------------------------------------------------
IMPLICIT NONE
c include COMMON blocks
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -64738,7 +64738,7 @@ c local type definitions
GAMMAX = xgam
XMINA = 2.D0*STR_mass/SQS
IF(IPAR(73).eq.1.and.KINT.gt.1) GAMMAX = PAR(119)
- IF(ndebug.gt.3) THEN
+ IF(idebug.gt.3) THEN
WRITE(LUN,*)' IMRG2HAD: called with ',
& '(KRMNT,KINT,NSEA,XGAM,XJET,STR_MASS):',
& KRMNT,KINT,NSEA,XGAM,XJET,STR_MASS
@@ -64753,11 +64753,11 @@ c sample total fraction for sea partons..
Z=(Z1+AC)*(1.D0+R*(((Z2+AC)/(Z1+AC))**NSEA-1.D0))
& **(1.D0/DBLE(NSEA))-AC
XSEA = XMINA*EXP(Z)
- IF(ndebug.gt.3) WRITE(LUN,*) ' total SEA fraction:' , xsea
+ IF(idebug.gt.3) WRITE(LUN,*) ' total SEA fraction:' , xsea
IF ( (1.D0-XSEA)**GAMMAX .LT. S_RNDM(1)) GOTO 50
c maximal fraction remaining for valence..
60 XREM = XSEA - DBLE(Nsea)*XMINA
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& WRITE(LUN,*) ' Xsea,xval,xjet:',
& xsea,1.D0-XSEA-XJET,xjet
@@ -64773,7 +64773,7 @@ c IF(KINT.gt.1.and.j.gt.2*KRMNT) XA=SIGN(ABS(XA)**PAR(116),XA)
XX(jj) = XMINA + XA
c new remainder
XREM = XREM - XA
- IF(ndebug.gt.3) write(lun,*)' x1,j,rem,xa',xX(jj),jj,xrem,xa
+ IF(idebug.gt.3) write(lun,*)' x1,j,rem,xa',xX(jj),jj,xrem,xa
ENDDO
c last parton..
ilast = 2+Nsea
@@ -64785,13 +64785,13 @@ c first interaction takes most energy
IF(KINT.gt.1.and.IPAR(71).eq.1)THEN
JJ = 3
IF(KRMNT.eq.0) JJ = 5
- if(ndebug.gt.4) write(lun,*) ' x1+x2,p*xeq:',
+ if(idebug.gt.4) write(lun,*) ' x1+x2,p*xeq:',
& XX(JJ)+XX(JJ+1),PAR(117)*XSEA/KINT
IF(XX(JJ)+XX(JJ+1).lt.PAR(117)*XSEA/KINT) GOTO 60
ENDIF
XSJ = XSJ + XSEA
- IF(ndebug.gt.3)THEN
+ IF(idebug.gt.3)THEN
write(lun,*)' x1,N,rem',xx(ilast),ilast,xrem
write(lun,*) ' xseajet',xsj
endif
@@ -64853,8 +64853,8 @@ c external types
INTEGER K_beam, NW, Irec, IREJ
c COMMONs
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NW_max
PARAMETER (NW_max = 20)
C--------------------------------------------------------------------
@@ -64956,7 +64956,7 @@ c remember position on parton stack
c set interaction properties
c IF(Irec.ne.1) CALL INI_EVENT(ECM,K_beam,Idm,Irec)
- IF(ndebug.gt.0)then
+ IF(idebug.gt.0)then
IF(Irec.eq.0)THEN
WRITE(LUN,*)
& ' SIB_NDIFF: recursive call with (ecm,kb,kt,np,jdif):',
@@ -65022,7 +65022,7 @@ c NCALL = NCALL + 1
GOTO 2000
ELSEIF(LXBAD.eq.1)THEN
c reject entire event
- if(Ndebug.gt.0)
+ if(Idebug.gt.0)
& WRITE(LUN,*)' SIB_NDIFF: minijet rejection (Ncall):',Ncall
c restore initial state
NP = NP_0
@@ -65046,7 +65046,7 @@ c NCALL = NCALL + 1
GOTO 2000
ELSEIF(LXBAD.eq.1)THEN
c reject entire event
- if(Ndebug.gt.0)
+ if(Idebug.gt.0)
& WRITE(LUN,*)' SIB_NDIFF: rmnt rejection (Ncall,NW):',Ncall,NW
c restore initial state
NP = NP_0
@@ -65074,8 +65074,8 @@ c restore initial state
S = SQS*SQS
RETURN
ENDIF
- IF(NDEBUG.gt.0) THEN
- IF(NDEBUG.gt.1) CALL PRNT_PRTN_STCK
+ IF(IDEBUG.gt.0) THEN
+ IF(IDEBUG.gt.1) CALL PRNT_PRTN_STCK
WRITE(LUN,*) ' SIB_NDIFF: entering fragmentation stage...'
ENDIF
@@ -65149,8 +65149,8 @@ c external type declarations
INTEGER KBEAM,NW,IREC,LBAD
C COMMONs
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -65209,7 +65209,7 @@ c internals
SAVE
DATA PREM /0.D0/ , PREM_NUC /0.D0/
- IF(Ndebug.gt.1)
+ IF(Idebug.gt.1)
& WRITE(LUN,*)' SAMPLE_RMNT: called with (Kbeam,NW,X1JET,',
& 'X2JET,JDIF,Irec):',Kbeam,NW,X1JET,(X2JET(JJ),JJ=1,NW),
& (JDIF(JJ),JJ=1,NW),Irec
@@ -65294,7 +65294,7 @@ c no target side excitation if recursive call (irec=0)!
ELSE
ITGRMNT(K) = 0
ENDIF
- IF(Ndebug.gt.1)
+ IF(Idebug.gt.1)
& WRITE(LUN,'(2X,A,1X,I2,1X,F5.3,1X,I2,1X,I2,1X,I2,1X,I2)')
& 'SAMPLE_RMNT: (JW,PREM,NS,NH,IBMRMNT,LTGRMNT):',
& K,PREM,NNSOF(k),NNJET(k),IBMRMNT,ITGRMNT(k)
@@ -65311,7 +65311,7 @@ c multiple de-excitations not possible..
IBMRMNT=MIN(IBMRMNT,1)
IBMRMNT=IABS(IBMRMNT-1)*Irec
ENDIF
- IF(Ndebug.gt.1)
+ IF(Idebug.gt.1)
& WRITE(LUN,*)
& ' SAMPLE_RMNT: remnant sampling (PREM,NW,LBMRMNT,LTGRMNT): ',
& PREM,NW,IBMRMNT,(ITGRMNT(j),j=1,NW)
@@ -65328,7 +65328,7 @@ c retry without counting
c 22 CONTINUE
IF(ITRY(2).gt.NREJ(2))THEN
LBAD = 2
- IF(ndebug.gt.1)then
+ IF(idebug.gt.1)then
WRITE(LUN,*)' SAMPLE_RMNT: number of trials exceeded'
WRITE(LUN,*)' resample minijets...(IREJ,NW,NCALL)',
& LBAD, NW, NCALL
@@ -65509,8 +65509,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION PI,TWOPI,CMBARN
COMMON /SIB_CST/ PI,TWOPI,CMBARN
DOUBLE PRECISION AM,AM2
@@ -65663,8 +65663,8 @@ C-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DIMENSION SIGDIF(3)
DIMENSION Xpos1(96),Xwgh1(96),Xpos2(96),Xwgh2(96)
DOUBLE PRECISION xil,xiu,tl,tu
@@ -66292,8 +66292,8 @@ C*********************************************************************
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
IMPLICIT INTEGER(I-N)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION P,E
SAVE
@@ -66491,8 +66491,8 @@ C------------------------------------------------
c external types
INTEGER KF,KP1,KP2
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -70332,8 +70332,8 @@ C-----------------------------------------------------------------------
COMMON /CSPA/ ICSPA2(3)
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
C--------------------------------------------------------------------
C SIBYLL utility common blocks containing constants \FR'14
@@ -70352,7 +70352,7 @@ C--------------------------------------------------------------------
DATA AVOG /6.0221367D-04/
DATA ATARGET /14.514D0/ ! effective masss of air
- IF(NDEBUG.gt.0)
+ IF(IDEBUG.gt.0)
& write(lun,*) ' SIGMA_INI: using cross section model no.',
& (ICSPA2(i),i=1,3)
@@ -70402,7 +70402,7 @@ C...K-air
ALINT(J,3) = 1.D0/(AVOG*SSIGA(J,3)/ATARGET)
ENDDO
- if (ndebug .gt. 0 ) THEN
+ if (idebug .gt. 0 ) THEN
WRITE(LUN,'(1X,A)')
& ' SIGMA_INI: NUCLIB interaction lengths [g/cm**2]'
WRITE(LUN,'(1X,A)')
@@ -70532,8 +70532,8 @@ c external types
DIMENSION P1(5),P2(5)
INTEGER IMOD,LBAD
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NIPAR_max,NPAR_max
PARAMETER (NPAR_max=200,NIPAR_max=100)
DOUBLE PRECISION PAR
@@ -70567,7 +70567,7 @@ C--------------------------------------------------------------------
DATA ASLOP /31.10362D0/ ! fit to the slope parameter.
DATA BSLOP /-15.29012D0/
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*) ' TRANSFONSHELL: called with (Ecm,M1,M2,XMAX):',
& ECM,XM1in,XM2in,XMAX
@@ -70625,14 +70625,14 @@ c no pt
PTSWTCH = 0.D0
SLOPE = 1.D0
ENDIF
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& WRITE(LUN,*) ' TRANSFONSHELL: (SLOP0,SLOPE,log(M**2)):',
& SLOP0,SLOPE,ALX
T = -DLOG(MAX(EPS10,S_RNDM(0)))/SLOPE
PTS = T*X1*PTSWTCH
PZB2 = S*0.5D0*0.5D0*X1*X1-XMB2-PTS
PZT2 = S*0.5D0*0.5D0*X2*X2-XMT2-PTS
- IF(NDEBUG.gt.3)
+ IF(IDEBUG.gt.3)
& WRITE(LUN,*) ' TRANSFONSHELL: (PTS,PZB2,PZT2):',PTS,PZB2,PZT2
c IF (ABS(PZB2)-PZT2.GT.EPS10) GOTO 60
IF (PZB2.lt.0.D0.or.PZT2.LT.0.D0) GOTO 60
@@ -70652,7 +70652,7 @@ c IF (ABS(PZB2)-PZT2.GT.EPS10) GOTO 60
P1(ii) = -P2(ii)
enddo
P1(5) = XMB
- IF(NDEBUG.gt.3) THEN
+ IF(IDEBUG.gt.3) THEN
WRITE(LUN,*) ' TRANSFONSHELL: (P1):',(p1(ii),ii=1,5)
WRITE(LUN,*) ' TRANSFONSHELL: (P2):',(p2(ii),ii=1,5)
ENDIF
@@ -70673,8 +70673,8 @@ Cf2py double precision, intent(out) :: X1,X2,PT
IMPLICIT NONE
c include COMMONs
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
INTEGER NW_max
PARAMETER (NW_max = 20)
C--------------------------------------------------------------------
@@ -70717,14 +70717,14 @@ c internal types
DOUBLE PRECISION XMINA,XM2DIS,XR,SLOPE,S_RNDM,XRNDM
SAVE
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_SEA: alpha,asup,qmass,xmax',
& ALPHA,ASUP,XMASS,XMAX
c min. momentum fraction for massive quarks
c i.e. sample from 1/(x+x_min)
XMINA = 2.D0*XMASS/SQS
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_SEA: xmina:',XMINA
c exponent of large x suppression: (1-x)**b, b=0 or b>1
IF(ABS(ASUP).lt.EPS3)THEN
@@ -70739,7 +70739,7 @@ c quark
100 X1 = XM2DIS(XMINA,XMAX,ALPHA) ! ~(1/x)**alpha
XR = LOG(1.D0-X1)-LOG(1.D0-XMINA)
XRNDM = S_RNDM(1)
- IF(ndebug.gt.4)
+ IF(idebug.gt.4)
& write(lun,*) ' X1,XR,SLOPE*XR:',X1,XR,SLOPE*XR
if(SLOPE*XR.le.LOG(max(XRNDM,eps10))) goto 100
@@ -70747,7 +70747,7 @@ c anti-quark
200 X2 = XM2DIS(XMINA,XMAX,ALPHA) ! ~(1/x)**alpha
XR = log(1.D0-X2)-log(1.D0-XMINA)
XRNDM = S_RNDM(2)
- IF(ndebug.gt.4)
+ IF(idebug.gt.4)
& write(lun,*) ' X2,XR,SLOPE*XR,XRNDM:',
& X2,XR,SLOPE*XR,XRNDM
if(SLOPE*XR.le.log(max(XRNDM,eps10))) goto 200
@@ -70760,7 +70760,7 @@ c anti-quark
c sample pt
c not yet implemented... to avoid problem with virtual partons
pt = 0.D0
- IF(ndebug.gt.3)
+ IF(idebug.gt.3)
& write(lun,*) ' SAMPLE_SEA: X1,X2,PT:',X1,X2,PT
END
@@ -73149,8 +73149,8 @@ c-----------------------------------------------------------------------
implicit double precision (a-h,o-z)
#include "conex.h"
#include "conex.incnex"
- INTEGER NCALL, NDEBUG, LUN
- COMMON /S_DEBUG/ NCALL, NDEBUG, LUN
+ INTEGER NCALL, IDEBUG, LUN
+ COMMON /S_DEBUG/ NCALL, IDEBUG, LUN
DOUBLE PRECISION CBR
INTEGER KDEC,LBARP,IDB
COMMON /S_CSYDEC/ CBR(223+16+12+8), KDEC(1338+6*(16+12+8)),
@@ -73168,10 +73168,10 @@ c-----------------------------------------------------------------------
#if !__CXCORSIKA__ && !__CORSIKA8__
- Ndebug=0
+ Idebug=0
#endif
#ifdef __CXDEBUG__
- if(isx.ge.8)Ndebug=isx-7
+ if(isx.ge.8)Idebug=isx-7
#endif
#if !__CXCORSIKA__ && !__CORSIKA8__