1 |
{ |
2 |
The following was generated by the command: |
3 |
|
4 |
tabfunc -i rho1 rho2 rho3 rho4 rho5 rho6 rho7 rho8 < radfactor.dat |
5 |
|
6 |
Used by scripts illumcal, reflcal. |
7 |
} |
8 |
interp_arr`(x,f)=if(x-1,if(f(0)-x,interp_arr2`(floor(x),x,f),f(f(0))),f(1)); |
9 |
interp_arr2`(i,x,f)=(i+1-x)*f(i)+(x-i)*f(i+1); |
10 |
rho1(x)=interp_arr`((x-380)/10+1,rho1`p2810); |
11 |
rho1`p2810(i)=select(i,0.219,0.252,0.256,0.252,0.244, |
12 |
0.237,0.23,0.225,0.22,0.216,0.214,0.216,0.223,0.226,0.225, |
13 |
0.227,0.236,0.253,0.272,0.298,0.341,0.39,0.424,0.442,0.45, |
14 |
0.451,0.451,0.45,0.451,0.453,0.455,0.458,0.462,0.464,0.466, |
15 |
0.466,0.467,0.467,0.467,0.467,0.467); |
16 |
rho2(x)=interp_arr`((x-380)/10+1,rho2`p2810); |
17 |
rho2`p2810(i)=select(i,0.07,0.089,0.111,0.118,0.121, |
18 |
0.122,0.123,0.127,0.131,0.138,0.15,0.174,0.207,0.242,0.26, |
19 |
0.267,0.272,0.282,0.299,0.322,0.335,0.341,0.342,0.342,0.341, |
20 |
0.339,0.338,0.336,0.334,0.332,0.331,0.329,0.328,0.326,0.324, |
21 |
0.324,0.322,0.32,0.316,0.315,0.314); |
22 |
rho3(x)=interp_arr`((x-380)/10+1,rho3`p2810); |
23 |
rho3`p2810(i)=select(i,0.065,0.07,0.073,0.074,0.074, |
24 |
0.073,0.073,0.074,0.077,0.085,0.109,0.148,0.198,0.241,0.278, |
25 |
0.339,0.392,0.4,0.38,0.349,0.315,0.285,0.264,0.252,0.241,0.229, |
26 |
0.22,0.216,0.219,0.23,0.251,0.288,0.34,0.39,0.431,0.46,0.481, |
27 |
0.493,0.5,0.505,0.516); |
28 |
rho4(x)=interp_arr`((x-380)/10+1,rho4`p2810); |
29 |
rho4`p2810(i)=select(i,0.074,0.093,0.116,0.124,0.128, |
30 |
0.135,0.144,0.161,0.186,0.229,0.281,0.332,0.37,0.39,0.395, |
31 |
0.385,0.367,0.341,0.312,0.28,0.247,0.214,0.185,0.169,0.16, |
32 |
0.154,0.151,0.148,0.148,0.151,0.158,0.165,0.17,0.17,0.166, |
33 |
0.164,0.168,0.177,0.185,0.192,0.197); |
34 |
rho5(x)=interp_arr`((x-380)/10+1,rho5`p2810); |
35 |
rho5`p2810(i)=select(i,0.295,0.31,0.313,0.319,0.326, |
36 |
0.334,0.346,0.36,0.381,0.403,0.415,0.419,0.413,0.403,0.389, |
37 |
0.372,0.353,0.331,0.308,0.284,0.26,0.232,0.12,0.194,0.185,0.18, |
38 |
0.176,0.175,0.175,0.18,0.186,0.192,0.199,0.199,0.196,0.195, |
39 |
0.197,0.203,0.208,0.215,0.219); |
40 |
rho6(x)=interp_arr`((x-380)/10+1,rho6`p2810); |
41 |
rho6`p2810(i)=select(i,0.151,0.265,0.41,0.492,0.517, |
42 |
0.531,0.544,0.556,0.554,0.541,0.519,0.488,0.45,0.414,0.377, |
43 |
0.341,0.309,0.279,0.253,0.234,0.225,0.221,0.22,0.22,0.223, |
44 |
0.233,0.244,0.258,0.268,0.278,0.283,0.291,0.302,0.325,0.351, |
45 |
0.376,0.401,0.425,0.447,0.469,0.485); |
46 |
rho7(x)=interp_arr`((x-380)/10+1,rho7`p2810); |
47 |
rho7`p2810(i)=select(i,0.378,0.524,0.551,0.559,0.561, |
48 |
0.556,0.544,0.522,0.488,0.448,0.408,0.363,0.324,0.301,0.283, |
49 |
0.265,0.257,0.259,0.26,0.256,0.254,0.27,0.302,0.344,0.377,0.4, |
50 |
0.42,0.438,0.452,0.462,0.468,0.473,0.483,0.496,0.511,0.525, |
51 |
0.539,0.553,0.565,0.575,0.581); |
52 |
rho8(x)=interp_arr`((x-380)/10+1,rho8`p2810); |
53 |
rho8`p2810(i)=select(i,0.104,0.17,0.319,0.462,0.49, |
54 |
0.482,0.462,0.439,0.413,0.382,0.352,0.325,0.299,0.283,0.27, |
55 |
0.256,0.25,0.254,0.264,0.272,0.278,0.295,0.348,0.434,0.528, |
56 |
0.604,0.648,0.676,0.693,0.705,0.712,0.717,0.721,0.719,0.725, |
57 |
0.729,0.73,0.73,0.73,0.73,0.73); |