COSMOS v7.655  COSMOSv7655
(AirShowerMC)
mninpu.f
Go to the documentation of this file.
1 *
2 * $Id: mninpu.F,v 1.1.1.1 1996/03/07 14:31:30 mclareni Exp $
3 *
4 * $Log: mninpu.F,v $
5 * Revision 1.1.1.1 1996/03/07 14:31:30 mclareni
6 * Minuit
7 *
8 *
9  SUBROUTINE mninpu(IUNIT,IERR)
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 the user to SET INPUT to IUNIT,
25 CC an alternative to MNSTIN where the user can specify just
26 CC a logical unit number and he is not interrogated about
27 CC open files and rewinding, all that is the responsibility
28 CC of the user and cannot be fixed interactively.
29 *
30 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
31 *
32 * $Log: d506cm.inc,v $
33 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
34 * Minuit
35 *
36 *
37 *
38 *
39 * d506cm.inc
40 *
41  parameter(mne=100 , mni=50)
42  parameter(mnihl=mni*(mni+1)/2)
43  CHARACTER*10 CPNAM
44  COMMON
45  1/mn7nam/ cpnam(mne)
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)
57 C
58  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
59  parameter(zero=0.0, one=1.0, half=0.5)
60  COMMON
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)
68  j/mn7arg/ word7(maxp)
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
77 C
78  ierr = 0
79 C IUNIT = 0, revert to previous input file
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'
84  ELSE
85  isysrd = istkrd(nstkrd)
86  nstkrd = nstkrd - 1
87  ENDIF
88 C
89 C new input file
90  ELSE
91  IF (nstkrd .GE. maxstk) THEN
92  WRITE (isyswr, '(A)') ' INPUT FILE STACK SIZE EXCEEDED.'
93  GO TO 800
94  ENDIF
95  nstkrd = nstkrd + 1
96  istkrd(nstkrd) = isysrd
97  isysrd = iunit
98  ENDIF
99 C
100  RETURN
101  800 ierr = 1
102  RETURN
103  END
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
z
nodes z
Definition: cstdatmos0.fNew2.h:93
up
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/
e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
i
nodes i
Definition: cstdatmos0.fNew2.h:93
xs
! 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
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
o
struct ob o[NpMax]
Definition: Zprivate.h:34
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
a
nodes a
Definition: cstdatmos0.fNew2.h:93
b
dE dx *! Nuc Int sampling table b
Definition: cblkMuInt.h:130
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
h
dE dx *! Nuc Int sampling table h
Definition: cblkMuInt.h:130
g
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
maxp
!onst int maxp
Definition: Zprivate.h:3
n
integer n
Definition: Zcinippxc.h:1
mninpu
subroutine mninpu(IUNIT, IERR)
Definition: mninpu.f:10
softenpik::half
integer, parameter half
Definition: csoftenPiK.f:108
x
! 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
f
dE dx *! Nuc Int sampling table f
Definition: cblkMuInt.h:130
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130