This commit was generated by cvs2svn to compensate for changes in r7641,
[u/mrichter/AliRoot.git] / DPMJET / flukapro / (ISOTOP)
1 *$ CREATE ISOTOP.ADD
2 *COPY ISOTOP
3 *
4 *=== isotop ===========================================================*
5 *
6 *----------------------------------------------------------------------*
7 *                                                                      *
8 *                                                                      *
9 *     Created on 23 september 1990 by    Alfredo Ferrari & Paola Sala  *
10 *                                                   Infn - Milan       *
11 *                                                                      *
12 *     Last change on 17-dec-96     by    Alfredo Ferrari               *
13 *                                                                      *
14 *     description of the common block(s) and variable(s)               *
15 *                                                                      *
16 *                                                                      *
17 *        Isondx = Initial and final indexes in the abuiso and isomnm   *
18 *                 for a given atomic number                            *
19 *        Isomnm = Mass numbers of the isotopes                         *
20 *        Abuiso = Natural abundances of the isotopes                   *
21 *        Astlin(1,iz) = "average" A of the stability line versus the   *
22 *                       atomic number Z                                *
23 *        Astlin(2,iz) = Dispersion of A of the stability line versus   *
24 *                       the atomic number Z                            *
25 *        Zstlin(1,ia) = "average" Z of the stability line versus the   *
26 *                       mass number A                                  *
27 *        Zstlin(2,ia) = Dispersion of Z of the stability line versus   *
28 *                       the mass number A                              *
29 *          Amssst(iz) = Atomic weight of iz_th element according to    *
30 *                       the tabulated mass tables                      *
31 *        Waps (ia,iz) = Atomic excess mass (MeV), for mass number Ia   *
32 *                       and atomic number given by Z = Z_start+Iz-1,   *
33 *                       where Z_start is stored into Inwaps (ia)       *
34 *       T12nuc(ia,iz) = Half-life (s) for mass number Ia and Z like    *
35 *                       for Waps                                       *
36 *       Brdecy(ia,iz) = Branching ratio for the first decay channel,   *
37 *                       for mass number Ia and Z like for Waps         *
38 *       Jspnuc(ia,iz) = Spin (hbar/2 units) for mass number Ia and Z   *
39 *                       like for Waps                                  *
40 *       Jptnuc(ia,iz) = Parity for mass number Ia and Z like for Waps  *
41 *     Idcnuc(k,ia,iz) = k_th (k=1,2) decay channel for mass number Ia  *
42 *                       and Z like for Waps                            *
43 *        Wapism (ism) = Atomic excess mass (MeV) of the ism_th isomer  *
44 *        T12ism (ism) = Half-life (s) of the ism_th isomer             *
45 *        Wapism (ism) = Atomic excess mass of the ism_th isomer        *
46 *        Bdcism (ism) = Branching ratio for the first decay channel    *
47 *                       of the ism_th isomer                           *
48 *        Jspism (ism) = Spin (hbar/2 units) of the ism_th isomer       *
49 *        Jptism (ism) = Parity (hbar/2 units) of the ism_th isomer     *
50 *       Idcism(k,ism) = k_th (k=1,2) decay channel for the ism_th      *
51 *                       isomer                                         *
52 *        Inwism (ia)  = (Cumulative) number of isomers with mass       *
53 *                       number =< Ia                                   *
54 *         IdcydA (j)  = Delta-A for the j_th type decay channel        *
55 *         IdcydZ (j)  = Delta-Z for the j_th type decay channel        *
56 *                                                                      *
57 *          Decay channel index:                                        *
58 *                                                                      *
59 *                    -1:        Unknown                                *
60 *                     0:        Stable                                 *
61 *                     1: IT     Internal Transition                    *
62 *                     2: A      Alpha                                  *
63 *                     3: N      Neutron emission                       *
64 *                     4: P      Proton  emission                       *
65 *                     5: 2P     Proton  + Proton emission              *
66 *                     6: 12C    12-C emission                          *
67 *                     7: 14C    14-C emission                          *
68 *                     8: B+     Beta+                                  *
69 *                     9: EC     Electron Capture                       *
70 *                    10: ECP    Electron Capture + Proton              *
71 *                    11: ECA    Electron Capture + Alpha               *
72 *                    12: EC2P   Electron Capture + Proton + Proton     *
73 *                    13: EC2A   Electron Capture + Alpha  + Alpha      *
74 *                    14: ECF    Electron Capture + Fission             *
75 *                    15: B-     Beta-                                  *
76 *                    16: B-N    Beta- + Neutron                        *
77 *                    17: B-A    Beta- + Alpha                          *
78 *                    18: B-2N   Beta- + Neutron + Neutron              *
79 *                    19: B-NA   Beta- + Neutron + Alpha                *
80 *                    20: B-2A   Beta- + Alpha   + Alpha                *
81 *                    21: B-B-   Beta- + Beta-                          *
82 *                    22: SF     Spontaneous Fission                    *
83 *                                                                      *
84 *----------------------------------------------------------------------*
85 *
86 *     PARAMETER ( NAMSMX = 250 )
87 *     PARAMETER ( NZGVAX =  10 )
88 *
89       PARAMETER ( NAMSMX = 270 )
90       PARAMETER ( NZGVAX =  16 )
91 *  Maximum number of isomers:
92       PARAMETER ( NISMMX = 591 )
93 *  Maximum number of types of decay channels:
94       PARAMETER ( NDCYMX =  22 )
95       COMMON / ISOTOP / WAPS   (NAMSMX,NZGVAX),  T12NUC (NAMSMX,NZGVAX),
96      &                  BRDECY (NAMSMX,NZGVAX),  WAPISM (NISMMX),
97      &                  T12ISM (NISMMX), BRDISM (NISMMX),
98      &                  ABUISO (NSTBIS), ASTLIN (2,100), ZSTLIN (2,260),
99      &                  AMSSST (100)  , ISOMNM (NSTBIS), ISONDX (2,100),
100      &                  JSPNUC (NAMSMX,NZGVAX), JPTNUC (NAMSMX,NZGVAX),
101      &                  IDCNUC (2,NAMSMX,NZGVAX), INWAPS (NAMSMX),
102      &                  JSPISM (NISMMX), JPTISM (NISMMX),
103      &                  IDCISM (2,NISMMX), IZWISM (NISMMX),
104      &                  INWISM (0:NAMSMX), IDCYDZ (NDCYMX),
105      &                  IDCYDA (NDCYMX)
106       COMMON / CHISTP / CHDECY (NDCYMX)
107       CHARACTER CHDECY*5
108