23 IMPLICIT DOUBLE PRECISION (
a-
h,
o-
z)
59 2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
60 3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
61 4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
62 5/mn7int/
x(mni) ,xt(mni) ,dirin(mni)
63 6/mn7fx2/
xs(mni) ,xts(mni) ,dirins(mni)
64 7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
65 8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
66 9/mn7fx1/ ipfix(mni) ,npfix
67 a/mn7var/ vhmat(mnihl)
68 b/mn7vat/ vthmat(mnihl)
69 c/mn7sim/
p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
71 parameter(maxdbg=10, maxstk=10, maxcwd=20,
maxp=30, maxcpt=101)
74 d/mn7npr/ maxint ,npar ,maxext ,nu
75 e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
76 e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
77 f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
78 g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
79 h/mn7min/ amin ,
up ,edm ,fval3 ,epsi ,apsi ,dcovar
80 i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
82 k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
83 l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
84 m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
85 n/mn7cpt/ chpt(maxcpt)
86 o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
87 CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
88 + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
89 LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
91 CHARACTER comand*(maxcwd)
92 CHARACTER clower*26, cupper*26
98 DATA clower/
'abcdefghijklmnopqrstuvwxyz'/
99 DATA cupper/
'ABCDEFGHIJKLMNOPQRSTUVWXYZ'/
107 DO 110
i= 1, min(maxcwd,lenbuf)
108 IF (crdbuf(
i:
i) .EQ.
'''')
GO TO 111
109 IF (crdbuf(
i:
i) .EQ.
' ')
THEN 110 IF (leader) ipos = ipos + 1
115 IF (crdbuf(
i:
i) .EQ. clower(ic:ic)) crdbuf(
i:
i)=cupper(ic:ic)
120 IF (ipos .GT. lenbuf)
THEN 121 WRITE (isyswr,
'(A)')
' BLANK COMMAND IGNORED.' 127 IF (crdbuf(ipos:ipos+2) .EQ.
'PAR')
THEN 133 IF (crdbuf(ipos:ipos+6) .EQ.
'SET INP')
THEN 139 IF (crdbuf(ipos:ipos+6) .EQ.
'SET TIT')
THEN 145 IF (crdbuf(ipos:ipos+6) .EQ.
'SET COV')
THEN 151 CALL mncrck(crdbuf(ipos:lenbuf),maxcwd,comand,lnc,
152 +
maxp, plist, llist, ierr,isyswr)
153 IF (ierr .GT. 0)
THEN 154 WRITE (isyswr,
'(A)')
' COMMAND CANNOT BE INTERPRETED' 159 CALL mnexcm(fcn,comand(1:lnc),plist,llist,ierr,futil)
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
subroutine mncrck(CRDBUF, MAXCWD, COMAND, LNC, MXP, PLIST, LLIST, IERR, ISYSWR)
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
subroutine mnexcm(FCN, COMAND, PLIST, LLIST, IERFLG, FUTIL)
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
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
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon ! knockon is considered Obsolete *PhotoProd false
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