COSMOS v7.655  COSMOSv7655
(AirShowerMC)
mncomd.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine mncomd (FCN, CRDBIN, ICONDN, FUTIL)
 

Function/Subroutine Documentation

◆ mncomd()

subroutine mncomd ( external  FCN,
character*(*)  CRDBIN,
  ICONDN,
external  FUTIL 
)

Definition at line 10 of file mncomd.f.

References a, b, c, d, e, f, false, g, h, softenpik::half, i, j, m, maxp, mncrck(), mnexcm(), n, o, p, parameter(), true, up(), x, xs, and z.

Referenced by mnread().

10 *
11 * $Id: d506dp.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
12 *
13 * $Log: d506dp.inc,v $
14 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
15 * Minuit
16 *
17 *
18 *
19 *
20 * d506dp.inc
21 *
22 C ************ DOUBLE PRECISION VERSION *************
23  IMPLICIT DOUBLE PRECISION (a-h,o-z)
24 CC Called by user. 'Reads' a command string and executes.
25 CC Equivalent to MNEXCM except that the command is given as a
26 CC character string.
27 CC
28 CC ICONDN = 0: command executed normally
29 CC 1: command is blank, ignored
30 CC 2: command line unreadable, ignored
31 CC 3: unknown command, ignored
32 CC 4: abnormal termination (e.g., MIGRAD not converged)
33 CC 5: command is a request to read PARAMETER definitions
34 CC 6: 'SET INPUT' command
35 CC 7: 'SET TITLE' command
36 CC 8: 'SET COVAR' command
37 CC 9: reserved
38 CC 10: END command
39 CC 11: EXIT or STOP command
40 CC 12: RETURN command
41 CC
42 *
43 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
44 *
45 * $Log: d506cm.inc,v $
46 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
47 * Minuit
48 *
49 *
50 *
51 *
52 * d506cm.inc
53 *
54  parameter(mne=100 , mni=50)
55  parameter(mnihl=mni*(mni+1)/2)
56  CHARACTER*10 cpnam
57  COMMON
58  1/mn7nam/ cpnam(mne)
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)
70 C
71  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
72  parameter(zero=0.0, one=1.0, half=0.5)
73  COMMON
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)
81  j/mn7arg/ word7(maxp)
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
90  dimension plist(maxp)
91  CHARACTER comand*(maxcwd)
92  CHARACTER clower*26, cupper*26
93  LOGICAL leader
94 C
95  EXTERNAL fcn,futil
96  CHARACTER*(*) crdbin
97  CHARACTER*100 crdbuf
98  DATA clower/'abcdefghijklmnopqrstuvwxyz'/
99  DATA cupper/'ABCDEFGHIJKLMNOPQRSTUVWXYZ'/
100 C
101  lenbuf = len(crdbin)
102  crdbuf = crdbin
103  icondn = 0
104 C record not case-sensitive, get upper case, strip leading blanks
105  leader = .true.
106  ipos = 1
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
111  GO TO 110
112  ENDIF
113  leader = .false.
114  DO 108 ic= 1, 26
115  IF (crdbuf(i:i) .EQ. clower(ic:ic)) crdbuf(i:i)=cupper(ic:ic)
116  108 CONTINUE
117  110 CONTINUE
118  111 CONTINUE
119 C blank or null command
120  IF (ipos .GT. lenbuf) THEN
121  WRITE (isyswr,'(A)') ' BLANK COMMAND IGNORED.'
122  icondn = 1
123  GO TO 900
124  ENDIF
125 C . . preemptive commands
126 C if command is 'PARAMETER'
127  IF (crdbuf(ipos:ipos+2) .EQ. 'PAR') THEN
128  icondn = 5
129  lphead = .true.
130  GO TO 900
131  ENDIF
132 C if command is 'SET INPUT'
133  IF (crdbuf(ipos:ipos+6) .EQ. 'SET INP') THEN
134  icondn = 6
135  lphead = .true.
136  GO TO 900
137  ENDIF
138 C if command is 'SET TITLE'
139  IF (crdbuf(ipos:ipos+6) .EQ. 'SET TIT') THEN
140  icondn = 7
141  lphead = .true.
142  GO TO 900
143  ENDIF
144 C if command is 'SET COVARIANCE'
145  IF (crdbuf(ipos:ipos+6) .EQ. 'SET COV') THEN
146  icondn = 8
147  lphead = .true.
148  GO TO 900
149  ENDIF
150 C crack the command . . . . . . . . . . . . . . . .
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'
155  icondn = 2
156  GO TO 900
157  ENDIF
158 C
159  CALL mnexcm(fcn,comand(1:lnc),plist,llist,ierr,futil)
160  icondn = ierr
161  900 RETURN
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
nodes z
subroutine mncrck(CRDBUF, MAXCWD, COMAND, LNC, MXP, PLIST, LLIST, IERR, ISYSWR)
Definition: mncrck.f:11
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
Definition: cblkMuInt.h:130
subroutine mnexcm(FCN, COMAND, PLIST, LLIST, IERFLG, FUTIL)
Definition: mnexcm.f:25
nodes i
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
Definition: cblkElemag.h:7
real(4), dimension(:), allocatable, save h
Definition: cNRLAtmos.f:28
! 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
Definition: Zglobalc.h:18
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
struct ob o[NpMax]
Definition: Zprivate.h:34
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
real(4), save a
Definition: cNRLAtmos.f:20
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
Definition: cblkElemag.h:7
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
real(4), save b
Definition: cNRLAtmos.f:21
!onst int maxp
Definition: Zprivate.h:3
integer n
Definition: Zcinippxc.h:1
integer, parameter half
Definition: csoftenPiK.f:108
! 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
Definition: Zptcl.h:21
dE dx *! Nuc Int sampling table f
Definition: cblkMuInt.h:130
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function:
Here is the caller graph for this function: