*$ CREATE ISOTOP.ADD *COPY ISOTOP * *=== isotop ===========================================================* * *----------------------------------------------------------------------* * * * Copyright (C) 1990-2011 by Alfredo Ferrari & Paola Sala * * All Rights Reserved. * * * * * * Created on 23 September 1990 by Alfredo Ferrari & Paola Sala * * Infn - Milan * * * * Last change on 31-Jan-11 by Alfredo Ferrari * * * * description of the common block(s) and variable(s) * * * * * * Isondx(k,iz) = Initial (k=1) and final (k=2) indexes in the * * abuiso and isomnm arrays for a given atomic * * number (iz) * * Isomnm(i) = Mass number of the i_th stable isotope * * Isoznm(i) = Atomic number of the i_th stable isotope * * Isqndx(k,iz) = Initial (k=1) and final (k=2) indexes in the * * isqmnm array for a given atomic number (iz) * * Lzmetl(iz) = Logical flag identifying metallic elements * * Izmetl(iz) = Number of conduction band electrons for a * * given atomic number (iz) * * Isqmnm(i) = Mass number of the i_th "quasi" stable * * isotope * * Isqznm(i) = Atomic number of the i_th "quasi" stable * * isotope * * Abuiso(i) = Natural abundances of the i_th stable isotope * * Amnsti(i) = Nuclear mass (GeV) of the i_th (quasi) stable * * isotope * * Ammsti(i) = Atomic mass (GeV) of the i_th (quasi) stable * * isotope * * Astlin(1,iz) = "average" A of the stability line versus the * * atomic number Z * * Astlin(2,iz) = Dispersion of A of the stability line versus * * the atomic number Z * * Zstlin(1,ia) = "average" Z of the stability line versus the * * mass number A * * Zstlin(2,ia) = Dispersion of Z of the stability line versus * * the mass number A * * Amssst(iz) = Atomic weight of iz_th element according to * * the tabulated mass tables * * Waps (ia,iz) = Atomic excess mass (MeV), for mass number Ia * * and atomic number given by Z = Z_start+Iz-1, * * where Z_start is stored into Inwaps (ia) * * Tfexms(ia,iz) = Thomas-Fermi (by Myers & Swiatecki) atomic * * excess mass (MeV), for mass number Ia and * * atomic number given by Z = Z_start+Iz-1, where * * Z_start is stored into Inwaps (ia) * * Shllms(ia,iz) = shell correction term (MeV) by Myers & Swiat- * * ecki for mass number Ia and atomic number given* * by Z = Z_start+Iz-1, where Z_start is stored * * into Inwaps (ia) * * Shllnm(ia,iz) = shell correction term (MeV) by Nix & Moller * * for mass number Ia and atomic number given * * by Z = Z_start+Iz-1, where Z_start is stored * * into Inwaps (ia) * * Dfenms(ia,iz) = deformation energy (MeV) by Myers & Swiatecki * * for mass number Ia and atomic number given * * by Z = Z_start+Iz-1, where Z_start is stored * * into Inwaps (ia) * * T12nuc(ia,iz) = Half-life (s) for mass number Ia and Z like * * for Waps * * Thrnci(k,ia,iz) = k_th (k=0,1,2->photon/neutron/proton) inelastic* * reaction minimum threshold (GeV) for mass num- * * ber Ia and Z like for Waps (!!! Note GeV !!!) * * Brdecy(k,ia,iz) = Branching ratio for the k_th decay channel, * * for mass number Ia and Z like for Waps * * Inwaps (ia) = minimum Z for which infos are recorded for a * * given A * * Inwexp(2,ia) = minimum and maximum Z for which experimental * * infos are available for a given A * * Inwadm(2,ia) = minimum and maximum Z for which FRDM mass * * infos are available for a given A * * Jspnuc(ia,iz) = Spin (hbar/2 units) for mass number Ia and Z * * like for Waps * * Jptnuc(ia,iz) = Parity for mass number Ia and Z like for Waps * * Idcnuc(k,ia,iz) = k_th (k=1,...,4) decay channel for mass number * * Ia and Z like for Waps * * Ithnci(k,ia,iz) = k_th (k=0,1,2->photon/neutron/proton) inelastic* * reaction threshold index (0=photon,1=neutron, * * 2=proton,3=deuteron,4=triton,5=3-He,6=4-He) for* * mass number Ia and Z like for Waps * * Nexlev(ia,iz) = number of excited levels for mass number Ia * * and Z like for Waps (including the ground * * state) * * Nmxclv(ia,iz) = number of excited levels up to which the levels* * are complete for mass number Ia and Z like for * * Waps (including the ground state) * * Nmxslv(ia,iz) = number of excited levels up to which the spin * * parity assignments are unique for mass number * * Ia and Z like for Waps (including the ground * * state) * * Kdxnuc(ia,iz) = starting location in blank common (i*4, 0 * * address) for the gamma deexcitation data for * * mass number Ia and Z like for Waps * * The data consist of Nexlev(ia,iz) energy levels* * each jth one with j-1 branching ratios to the * * underlying j-1 levels. The partial number of * * points is Sum^Nexlev_1[j] if Nexlev < Nexmxb, * * Sum^Nexmxb_1[j]+(Nexlev(ia,iz)-Nexmxb)*(Nexmxb+1)* * otherwise. * * They are followed by Nexlev(ia,iz) level width * * (GeV) values. * * They are followed by Nexlev(ia,iz) integer * * flags * * flag = dcy1+10^2*dcy2+10^4*dcy3+10^6*dcy4 * * each jth one with j-1 further flags * * flag = lev_fin+10000*int(br_ic*10^5) * * for the transitions to the underlying j-1 * * levels. * * The partial number of points for this part is * * Sum^Nexlev_1[j] if Nexlev < Nexmxb, * * Sum^Nexmxb_1[j]+(Nexlev(ia,iz)-Nexmxb)*(Nexmxb+1)* * otherwise. * * They are followed by Nexlev(ia,iz) integers, * * each one giving: * * spin1 + spin2 * 100 + (parity + 1) * 10000 * * where spin1-spin2 is the spin range of that * * level in hbar/2 units (for well known levels * * spin1=spin2 of course) and parity is the parity* * (-1,1 or 0 if unknown) * * Therefore the total amount of data is: * * 2 Sum^Nexlev_1[j] + 2 Nexlev(ia,iz) * * (Nexlev(ia,iz) =< Nexmxb) * * or * * 2 Sum^Nexmxb_1[j] * * + 2 (Nexlev(ia,iz) - Nexmxb) * (Nexmxb + 1) * * + 2 Nexlev(ia,iz) * * (Nexlev(ia,iz) > Nexmxb) * * (note that everything is stored as single pre- * * cision data) * * Ngmlns(ia,iz) = number of recorded gamma lines for mass number * * Ia and Z like for Waps * * Kgmlns(ia,iz) = starting location in blank common (i*4, 0 * * address) for the gamma lines data for mass * * number Ia and Z like for Waps. * * The data consist of Ngmlns(ia,iz) data sets * * each one with the line branching, and the line * * energy (GeV). Negative branchings mean 511 keV * * gammas from positron annihilation. * * (note that everything is stored as single pre- * * cision data) * * Ncelns(ia,iz) = number of recorded CE lines for mass number * * Ia and Z like for Waps * * Kcelns(ia,iz) = starting location in blank common (i*4, 0 * * address) for the CE lines data for mass * * number Ia and Z like for Waps. * * The data consist of Ncelns(ia,iz) data sets * * each one with the line branching, and the line * * energy (GeV). * * (note that everything is stored as single pre- * * cision data) * * Nallns(ia,iz) = number of recorded alpha lines for mass number * * Ia and Z like for Waps * * Kallns(ia,iz) = starting location in blank common (i*4, 0 * * address) for the alpha lines data for mass * * number Ia and Z like for Waps. * * The data consist of Nallns(ia,iz) data sets * * each one with the line branching, and the line * * energy (GeV). * * (note that everything is stored as single pre- * * cision data) * * Nbtspc(ia,iz) = number of recorded beta spectra for mass number* * Ia and Z like for Waps * * Kbtspc(ia,iz) = starting location in blank common (i*4, 0 * * address) for the beta spectra data for mass * * number Ia and Z like for Waps. * * The data consist of Nbtspc(ia,iz) data sets * * each one with the spectrum branching, the spec-* * trum average energy (GeV), and the spectrum end* * -point energy (GeV). The average energy sign * * is giving the electron sign. * * (note that everything is stored as single pre- * * cision data) * * Ngxbzr(ia,iz) = number of energy intervals for (g,x) cross * * section Bezier fits for mass number Ia and * * Z like for Waps * * Kgxbzr(ia,iz) = starting location in blank common (i*4, 0 * * address) for the (g,x) cross section Bezier fit* * data for mass number Ia and Z like for Waps * * The data consist of Ngxbzr(ia,iz) energy limits* * each one with 4 Bezier coefficients and 1 fun- * * ction maximum (the energy limit is the last of * * the 6 data). The energy limits are actually * * Ngxbzr(ia,iz)+1, the first one being the star- * * ting threshold. * * Therefore the total amount of data is: * * 6 x Ngxbzr(ia,iz) + 1 * * (note that everything is stored as single pre- * * cision data) * * Kqmdin(ia,iz) = starting location in blank common (i*4, 0 * * address) for the QMD config. initialization * * data for mass number Ia and Z like for Waps * * Izwism (ism) = Z of the ism_th isomer * * Wapism (ism) = Atomic excess mass (MeV) of the ism_th isomer * * T12ism (ism) = Half-life (s) of the ism_th isomer * * Wapism (ism) = Atomic excess mass of the ism_th isomer * * Bdcism (k,ism) = Branching ratio for the k_th decay channel * * of the ism_th isomer * * Jspism (ism) = Spin (hbar/2 units) of the ism_th isomer * * Jptism (ism) = Parity (hbar/2 units) of the ism_th isomer * * Idcism(k,ism) = k_th (k=1,...,4) decay channel for the ism_th * * isomer * * Inwism (ia) = (Cumulative) number of isomers with mass * * number =< Ia * * Ngmism (ism) = number of recorded gamma lines for the ism_th * * isomer * * Kgmism (ism) = starting location in blank common (i*4, 0 * * address) for the gamma lines data for the * * ism_th isomer. * * The data consist of Ngmism(ism) data sets * * each one with the line branching, and the line * * energy (GeV). Negative branchings mean 511 keV * * gammas from positron annihilation. * * (note that everything is stored as single pre- * * cision data) * * Nceism (ism) = number of recorded CE lines for the ism_th * * isomer * * Kceism (ism) = starting location in blank common (i*4, 0 * * address) for the CE lines data for the * * ism_th isomer. * * The data consist of Nceism(ism) data sets * * each one with the line branching, and the line * * energy (GeV). * * (note that everything is stored as single pre- * * cision data) * * Nalism (ism) = number of recorded alpha lines for the ism_th * * isomer * * Kalism (ism) = starting location in blank common (i*4, 0 * * address) for the alpha lines data for the * * ism_th isomer. * * The data consist of Nalism(ism) data sets * * each one with the line branching, and the line * * energy (GeV). * * (note that everything is stored as single pre- * * cision data) * * Nbtism (ism) = number of recorded beta spectra for the ism_th * * isomer * * Kbtism (ism) = starting location in blank common (i*4, 0 * * address) for the beta spectra data for the * * ism_th isomer. * * The data consist of Nbtism(ism) data sets * * each one with the spectrum branching, the spec-* * trum end-point energy (GeV), and the spectrum * * dose * * (note that everything is stored as single pre- * * cision data) * * Is2zwa (is) = Z index in the Waps array for the is_th stable * * or "quasi" stable isotope (1