Changeset 3260 for branches

Show
Ignore:
Timestamp:
01/26/12 21:06:44 (4 months ago)
Author:
erkn
Message:

Sources: updated maxwellian parameters for ILL, HZB and FRM2

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • branches/mcxtrace-1.0/nlib/sources/Source_gen.comp

    r3256 r3260  
    6363* 
    6464* <b>Some sources parameters:</b> 
    65 * PSI cold source T1= 296.16, I1=8.5E11, T2=40.68, I2=5.2E11 
    66 * ILL VCS cold source T1=216.8,I1=1.24e+13,T2=33.9,I2=1.02e+13 
    67 *        (H1)         T3=16.7 ,I3=3.0423e+12 
    68 * ILL HCS cold source T1=413.5,I1=10.22e12,T2=145.8,I2=3.44e13 
    69 *        (H5)         T3=40.1 ,I3=2.78e13 
    70 * ILL Thermal tube    T1=683.7,I1=0.5874e+13,T2=257.7,I2=2.5099e+13 
    71 *        (H12)        T3=16.7 ,I3=1.0343e+12 
    72 * ILL Hot source      T1=1695,I1=1.74e13,T2=708,I2=3.9e12 
     65* PSI cold source     T1=296.2,I1=8.5E11, T2=40.68,I2=5.2E11 
     66* ILL VCS cold source T1=362.7,I1=7.95e12,T2=168.9,I2=1.23e+13, 
     67*        (H1)         T3=34.4 ,I3=1.12e13,yheight=0.22,xwidth=0.14, 
     68* ILL HCS cold source T1=311.9,I1=1.52e13,T2=131.7,I2=1.30e13, 
     69*        (H5)         T3=35.8 ,I3=1.27e13,radius=0.105 
     70* ILL Thermal tube    T1=229.6,I1=5.32e13,T2=1102, I2=4.37e12, 
     71*        (H12)        T3=437.1,I3=3.84e13,radius=0.1 
     72* ILL Hot source      T1=1695, I1=1.74e13,T2=708,  I2=3.9e12 
     73* HZB cold source     T1=43.7 ,I1=1.4e12, T2=137.2,I2=2.08e12,radius=.155 
     74* HZB bi-spectral     T1=43.7, I1=1.4e12, T2=137.2,I2=2.08e12,T3=293.0,I3=1.77e12 
     75* HZB thermal tube    T1=293.0,I1=2.64e12 
     76* FRM2 cold           T1=35.0, I1=9.38e12,T2=547.5,I2=2.23e12,T3=195.4,I3=1.26e13 
     77* FRM2 thermal        T1=285.6,I1=3.06e13,T2=3000, I2=1.68e12,T3=429.9,I3=6.77e12 
    7378* 
    7479* %VALIDATION 
    … …  
    104109* focus_ah: (deg) maximal (uniform) vert. height divergence 
    105110* T1:       [K]   Temperature of the Maxwellian source, 0=none 
    106 * I1:       [1/(cm**2*st*AA)] Source flux per solid angle, area and Angstrom 
     111* I1:       [1/(cm**2*sr*AA)] Source flux per solid angle, area and Angstrom 
    107112*                 if I1=0, the source emits 1 in 4*PI whole space. 
    108113* T2:       [K]   Second Maxwellian source Temperature, 0=none 
    109 * I2:       [1/(cm**2*st*AA)] Second Maxwellian Source flux 
     114* I2:       [1/(cm**2*sr*AA)] Second Maxwellian Source flux 
    110115* T3:       [K]   Third Maxwellian source Temperature, 0=none 
    111 * I3:       [1/(cm**2*st*AA)] Third Maxwellian Source flux 
     116* I3:       [1/(cm**2*sr*AA)] Third Maxwellian Source flux 
    112117* flux_file:[str] Name of a two columns [lambda flux] text file that contains 
    113118*                 the wavelength distribution of the flux in [1/(s*cm**2*st*AA)]. 
    … …  
    385390    dlambda = (Lmax-Lmin)/2; 
    386391  } 
    387   if ((lambda0 < dlambda) || (dlambda < 0)) 
     392  if (lambda0 <= 0 || (lambda0 < dlambda) || (dlambda < 0)) 
    388393  { fprintf(stderr,"Source_gen: %s: Error: Wavelength range %.3f +/- %.3f AA calculated \n", 
    389394      NAME_CURRENT_COMP, lambda0, dlambda);