*$ CREATE EMFBRM.ADD *COPY EMFBRM * *=== Emfbrm ===========================================================* * *----------------------------------------------------------------------* * * * Copyright (C) 1991-2011 by Alfredo Ferrari & Paola Sala * * All Rights Reserved. * * * * * * This is the common for the new EMF bremsstrahlung: the name been * * changed from Brempr (2005) to remark this * * * * * * Last change on 17-apr-09 by Alfredo Ferrari * * * * Mxxmde = maximum number of media in Emf * * * * Variables for the new bremsstrahlung: * * * * Thbrem = total electron/positron threshold energy for * * bremsstrahlung production (or minimum energy * * of emitted photons), MeV * * Althbr = natural logarithm of Thbrem - me * * Ebrm0 = minimum energy for which bremsstrahlung data * * are tabulated * * Albrm0 = natural logarithm of Ebrm0 * * Ebrmrt = ratio between subsequent tabulated energies * * after the initial 0:nktl-3 values * * Abrmrt = natural logarithm of Ebrmrt * * Aktild = k/T values corresponding at the tabulated en- * * ergies ( jth energy = (Thbrem-me) / Aktild (j) )* * Alktld = natural logarithm of Aktild * * Indsum = array used for addressing data * * Jndsum = array used for addressing data * * Ind0br = zero index address of brm. branching ratios * * and of brm. power fit exponents * * Ind1br = last index address of brm. branching ratios * * and of brm. power fit exponents * * Jnd0br = zero index address of brm. tabulated ds/dk * * Jnd1br = last index address of brm. tabulated ds/dk * * Lnwbrm = Logical flag for the new bremss. * * * *----------------------------------------------------------------------* * * PARAMETER ( NKTL = 17 ) * Until Fluka2010.4 included: * Changed on 17-Apr-2009 in order to make thinner the bremss. runtime * tabulations because of LPM ** PARAMETER ( NBRBIN = 75 ) ** PARAMETER ( BREXRT = ONEFOU ) * PARAMETER ( NBRBIN = 150 ) * PARAMETER ( BREXRT = ONEONE / EIGEIG ) PARAMETER ( NBRBIN = 160 ) PARAMETER ( BREXRT = ONEONE / EIGEIG ) * LOGICAL LNWBRM COMMON / EMFBRM / THBREM (MXXMDE), ALTHBR (MXXMDE), & EBRM0 (MXXMDE), ALBRM0 (MXXMDE), & EBRMRT (MXXMDE), ABRMRT (MXXMDE), & AKTILD (0:NBRBIN), ALKTLD (0:NBRBIN), & INDSUM (0:NBRBIN+1), IND0BR (2,MXXMDE), & IND1BR (2,MXXMDE), LND0BR (2,MXXMDE), & LND1BR (2,MXXMDE), JND0BR (2,MXXMDE), & JND1BR (2,MXXMDE), LNWBRM (MXXMDE) DIMENSION JNDSUM (0:NBRBIN) EQUIVALENCE ( INDSUM (1), JNDSUM (0) ) SAVE / EMFBRM /