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