*$ CREATE FLUOXR.ADD *COPY FLUOXR * *=== fluoxr ===========================================================* * *----------------------------------------------------------------------* * * * Copyright (C) 1990-2012 by Alfredo Ferrari * * All Rights Reserved. * * * * * * Include file Fluoxr: created on 4-October-1990 by A. Ferrari * * INFN - Milan * * * * Last change on 17-Feb-2012 by Alfredo Ferrari * * * * Description of the variables (not complete): * * * * Brrlxl(ix,il,iz) = Relative branching ratios of the * * ix-th Lil line (il=1-3, LI,LII,LIII) * * Brrkxl(ix,iz) = Relative branching ratios of the * * ix-th K line * * Flyldk(iz) = Fluorescence yield for the K-shell * * of the iz-th element * * Flyldl(il,iz) = Fluorescence yield for the il sub- * * shell of the L-shell (LI,LII,LII) * * of the iz-th element * * Zelshl(ks,iz) = Number of electrons in the ks_th sub-* * shell of the element iz * * Zceshl(ks,iz) = Number of conduction electrons (if * * any) in the ks_th sub-shell of the * * element iz * * Imnshl(ks) = Main shell (K,L,M,N,O,P,Q) of the * * ks_th sub-shell * * id=1 --> K * * id=2 - 4 --> L * * id=5 - 9 --> M * * id=10-16 --> N * * id=17-23 --> O * * id=24-28 --> P * * id=29-31 --> Q * * Ishell(id) = Main shell (K,L,M,N,O,P,Q) of the * * id-th level * * id=1 --> K * * id=2 - 4 --> L * * id=5 - 9 --> M * * id=10-16 --> N * * id=17-21 --> O * * id=22-27 --> P * * id=28-29 --> Q * * Ktoxtr(ix) = end level index for the ix-th X * * transitions from the K shell * * Ltoxtr(ix,il) = end level index for the ix-th X * * transitions from Lil subshell * * Nshell(is) = Number of sub-shell for the * * is-th (main) shell * * (is=1-7 --> K,L,M,N,O,P,Q) * * Nflndx(1,iz) = Last (sub)edge index for the iz ele- * * ment in the global numbering * * (Nlfndx(1,iz)-Nflndx(1,iz-1) is the * * total number of (sub)edges of the iz * * element for which a tabulated fit * * exists). These indeces are used in * * the Zphoto and Sgphpa arrays * * Nflndx(2,iz) = Last (main)edge (shell) index for the* * iz element in the global numbering * * (Nlfndx(2,iz)-Nflndx(2,iz-1) is the * * total number of tabulated main shells* * of the iz element). These indeces are* * used in the Azphot, Bzphot and Eiedge* * arrays (there are 455 main shells for* * all iz elements) * * Nflndx(3,iz) = Last (energy) level index for the * * iz element in the global numbering * * (Nlfndx(3,iz) - Nflndx(3,iz-1) is * * the total number of tabulated energy * * levels of the iz element). The total * * number of recorded energy levels is * * 1612 and includes also levels below * * 1 keV. * * The indeces between Nflndx(3,iz-1)+1 * * and Nlfndx(3,iz) point to the corres-* * ponding elements of the Ziedge array * * Ktbedg(ie) = Starting blank common location (0 * * address, 0th element is included) for* * edge-by-edge photoelectric tabulati- * * ons * * Ielshl(ks,iz) = is_th (edge) sub-shell corresponding * * to the ks_th subshell of the element * * iz * * Kelshl(is,iz) = ks_th subshell corresponding to the * * is_th (edge) sub-shell of the element* * iz * * Cheshl(ks,iz) = ks_th sub-shell nomenclature * * Lnwpex = logical flag for using the NeW edge- * * by-edge PhotoElectric X-section tabu-* * lations. If true, Nflndx(1,iz) and * * Nflndx(3,iz) coincides * * * * ks_th sub-shell: it corresponds to the sequence * * K (1S1/2) * * L1 (2S1/2) * * L2 (2P1/2) * * L3 (2P3/2) * * M1 (3S1/2) * * M2 (3P1/2) * * M3 (3P3/2) * * M4 (3D3/2) * * M5 (3D5/2) * * N1 (4S1/2) * * N2 (4P1/2) * * N3 (4P3/2) * * N4 (4D3/2) * * N5 (4D5/2) * * N6 (4F5/2) * * N7 (4F7/2) * * O1 (5S1/2) * * O2 (5P1/2) * * O3 (5P3/2) * * O4 (5D3/2) * * O5 (5D5/2) * * O6 (5F5/2) * * O7 (5F7/2) * * P1 (6S1/2) * * P2 (6P1/2) * * P3 (6P3/2) * * P4 (6D3/2) * * P5 (6D5/2) * * Q1 (7S1/2) * * Q2 (7P1/2) * * Q3 (7P3/2) * * is_th sub-shell: it corresponds to the same sequence * * skipping empty sub-shells, it is the * * sequence of photoelectric edges (and * * energy levels) * * * *----------------------------------------------------------------------* * PARAMETER ( EHYDR0 = 1.3605697953179906 D-05 ) PARAMETER ( SIGTHO = 6.6524615280813428 D-01 ) * PARAMETER ( TWOPIP = TWOTWO * PIPIPI ) * This is 4 sqrt(2) PARAMETER ( TWO3O2 = 5.656854249492380 D+00 ) * PARAMETER ( PHCONS = TWO3O2 * SIGTHO * FSCTO4 ) PARAMETER ( PHCONS = 1.0671336148164197D-08 ) PARAMETER ( PHCON2 = 64.D+00 * SIGTHO / FSCTO3 ) * PARAMETER ( PHCON2 = XXXXXXXXXXXXXXX ) * This is 2 x pi x exp (-4) PARAMETER ( PHEDG0 = 0.1150805531573017 D+00 ) * This is 3 / ( 8 x sqrt (2) ) PARAMETER ( SAUCST = 0.2651650429449553 D+00 ) * Parameters for new photoelectric tabulations: PARAMETER ( EMXED0 = 1.D+11 ) * PARAMETER ( NTBED0 = 15 ) * PARAMETER ( NTBED1 = 55 ) * PARAMETER ( NTBD12 = 2000 ) * PARAMETER ( Z0ED12 = 11.7 D+00 ) * PARAMETER ( S2ED12 = ONEONE ) * PARAMETER ( NTBD20 = 400 ) * PARAMETER ( Z0ED20 = 19.9 D+00 ) * PARAMETER ( S2ED20 = TWOTWO ) * PARAMETER ( NTBD27 = 350 ) * PARAMETER ( Z0ED27 = 27.0 D+00 ) * PARAMETER ( S2ED27 = SIXSIX ) PARAMETER ( NTBED0 = 15 ) PARAMETER ( NTBED1 = 58 ) PARAMETER ( NTBD12 = 2000 ) PARAMETER ( Z0ED12 = 11.35 D+00 ) PARAMETER ( S2ED12 = TWOTWO ) PARAMETER ( NTBD20 = 600 ) PARAMETER ( Z0ED20 = 19.7 D+00 ) PARAMETER ( S2ED20 = THRTWO ) PARAMETER ( NTBD27 = 450 ) PARAMETER ( Z0ED27 = 27.35 D+00 ) PARAMETER ( S2ED27 = TENTEN ) * PARAMETER ( EDGPOA = THRTHR ) PARAMETER ( EDGPOD = 0.0 D+00 ) * t = tabulation abscissa * x = E/E_edge * * 1/a * 1: t = ln[x] * a * E = E_edge exp[t ] * * * 2: t = ln[1+ln(x)] * * E = E_edge exp[exp(t)-1] * * 1/a * 3: t = ln[1+ln(x) ] * a * E = E_edge exp{[exp(t)-1] } * PARAMETER ( IEDGTB = 3 ) * LOGICAL LNWPEX COMMON / FLUOXR / ZPHOTO (920) , AZPHOT (455) , BZPHOT (455), & ZIEDGE (1612) , ZIEDFT (1612) , EIEDGE (455), & BRRLXL (4,3,100), BRRKXL (8,100), FLYLDK (100), & FLYLDL (3,100) , SGPHPA (1612) , HALLCR (100), & ASYPOW (1612) , AXMLMX (1612) , EHYDRO , & ZELSHL (MXESHL,100), ZCESHL (MXESHL,100), & EMXEDG, & KTOXTR (8) , LTOXTR (4,3) , & NFLNDX (3,0:100), NFLNFT (0:100), NSHELL (0:7), & ISHELL (0:29) , KTBEDG (1612) , NTBEDG(1612), & IDMIN (100) , IMNSHL (MXESHL), & IELSHL (MXESHL,100), KELSHL (MXESHL,100), & LNWPEX COMMON / CHFLUO / CHKXLN (8), CHLXLN (4,3), CHESHL (MXESHL), & CHMNSL (7) CHARACTER*8 CHKXLN, CHLXLN CHARACTER*10 CHESHL CHARACTER*1 CHMNSL SAVE / FLUOXR / SAVE / CHFLUO /