31 IMPLICIT DOUBLE PRECISION (
a-
h,
o-
z)
49 2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
50 3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
51 4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
52 5/mn7int/
x(mni) ,xt(mni) ,dirin(mni)
53 6/mn7fx2/
xs(mni) ,xts(mni) ,dirins(mni)
54 7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
55 8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
56 9/mn7fx1/ ipfix(mni) ,npfix
57 a/mn7var/ vhmat(mnihl)
58 b/mn7vat/ vthmat(mnihl)
59 c/mn7sim/
p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
61 parameter(maxdbg=10, maxstk=10, maxcwd=20,
maxp=30, maxcpt=101)
64 d/mn7npr/ maxint ,npar ,maxext ,nu
65 e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
66 e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
67 f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
68 g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
69 h/mn7min/ amin ,
up ,edm ,fval3 ,epsi ,apsi ,dcovar
70 i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
72 k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
73 l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
74 m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
75 n/mn7cpt/ chpt(maxcpt)
76 o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
77 CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
78 + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
79 LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
80 CHARACTER copt*1, corg*(*), cmes*(*), ctyp*7
82 CHARACTER origin(maxmes,2)*10, warmes(maxmes,2)*60
83 common/mn7wrc/origin, warmes
84 common/mn7wri/nfcwar(maxmes,2),icirc(2)
87 IF (corg(1:3).EQ.
'SHO' .AND. cmes(1:3).EQ.
'SHO')
GO TO 200
89 IF (copt .EQ.
'W')
THEN 92 WRITE (isyswr,
'(A,A/A,A)')
' MINUIT WARNING IN ',corg,
93 +
' ============== ',cmes
99 WRITE (isyswr,
'(A,A/A,A)')
' MINUIT DEBUG FOR ',corg,
100 +
' ============== ',cmes
105 IF (nwrmes(ityp) .EQ. 0) icirc(ityp) = 0
106 nwrmes(ityp) = nwrmes(ityp) + 1
107 icirc(ityp) = icirc(ityp) + 1
108 IF (icirc(ityp) .GT. maxmes) icirc(ityp) = 1
110 origin(ic,ityp) = corg
111 warmes(ic,ityp) = cmes
112 nfcwar(ic,ityp) = nfcn
117 IF (copt .EQ.
'W')
THEN 124 IF (nwrmes(ityp) .GT. 0)
THEN 125 englsh =
' WAS SUPPRESSED. ' 126 IF (nwrmes(ityp) .GT. 1) englsh =
'S WERE SUPPRESSED.' 127 WRITE (isyswr,
'(/1X,I5,A,A,A,A/)') nwrmes(ityp),
128 +
' MINUIT ',ctyp,
' MESSAGE', englsh
131 IF (nm .GT. maxmes)
THEN 132 WRITE (isyswr,
'(A,I2,A)')
' ONLY THE MOST RECENT ',
133 + maxmes,
' WILL BE LISTED BELOW.' 137 WRITE (isyswr,
'(A)')
' CALLS ORIGIN MESSAGE' 140 IF (ic .GT. maxmes) ic = 1
141 WRITE (isyswr,
'(1X,I6,1X,A,1X,A)')
142 + nfcwar(ic,ityp),origin(ic,ityp),warmes(ic,ityp)
145 WRITE (isyswr,
'(1H )')
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data h g *is for param c g data up(2, 1)/7.0d0/
dE dx *! Nuc Int sampling table e
real(4), dimension(:), allocatable, save h
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
dE dx *! Nuc Int sampling table d
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
dE dx *! Nuc Int sampling table g
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
dE dx *! Nuc Int sampling table f
dE dx *! Nuc Int sampling table c