23 IMPLICIT DOUBLE PRECISION (
a-
h,
o-
z)
55 2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
56 3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
57 4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
58 5/mn7int/
x(mni) ,xt(mni) ,dirin(mni)
59 6/mn7fx2/
xs(mni) ,xts(mni) ,dirins(mni)
60 7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
61 8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
62 9/mn7fx1/ ipfix(mni) ,npfix
63 a/mn7var/ vhmat(mnihl)
64 b/mn7vat/ vthmat(mnihl)
65 c/mn7sim/
p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
67 parameter(maxdbg=10, maxstk=10, maxcwd=20,
maxp=30, maxcpt=101)
70 d/mn7npr/ maxint ,npar ,maxext ,nu
71 e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
72 e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
73 f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
74 g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
75 h/mn7min/ amin ,
up ,edm ,fval3 ,epsi ,apsi ,dcovar
76 i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
78 k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
79 l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
80 m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
81 n/mn7cpt/ chpt(maxcpt)
82 o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
83 CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
84 + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
85 LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
90 IF (iuext .EQ. 0)
GO TO 100
91 IF (iuext .LT. 0)
THEN 94 IF (iint .GT. npar)
GO TO 100
100 IF (iext .EQ. 0)
GO TO 100
101 IF (iext .GT. nu)
GO TO 100
107 IF (nvl .LT. 0)
GO TO 100
110 IF (iint .GT. 0) err = werr(iint)
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