profile.inc [SRC] [CPP] [JOB] [SCAN] TOOLS / INCLUDE

```   1 | !     ==================================================================!
2 | !     The abscisa relative to the initial point is "x", you can use     !
3 | !     this value to generate the desired profile. Pressure is in [bar]  !
4 | !     ==================================================================!
5 |
6 | !     Profiles where calculated with cold wall at the left. This line
7 | !     puts the profiles in the "right" direction.
8 |
9 |       x = 2.0 - x
10 |
11 | !     Profiles calculated using a polynomial regression:
12 |
13 |       tempe(i) = 314.74 - 964.64 * x + 12418 * x**2 - 63841 * x**3 +    &
14 |      &           1.6589e+05 * x**4 - 2.4259e+05 * x**5 +                &
15 |      &           2.0819e+05 * x**6 - 1.0292e+05 * x**7 +                &
16 |      &           26964 * x**8 - 2889.2 * x**9
17 |
18 |       press(i) = 1.
19 |       XH2O(i)  = 0.001367 - 0.089525 * x + 1.1534 * x**2 - 5.9345 * x**3&
20 |      &         + 15.436 * x**4 - 22.606 * x**5 + 19.441 * x**6 -        &
21 |      &           9.6362 * x**7 + 2.5319 * x**8 - 0.27215 * x**9
22 |
23 |       XCO2(i)  = 0.00075944 - 0.049735 * x + 0.64075 * x**2 -           &
24 |      &           3.2969 * x**3 + 8.5755 * x**4 - 12.559 * x**5 +        &
25 |      &           10.8 * x**6 - 5.3533 * x** 7 + 1.4066 * x**8 -         &
26 |      &           0.15119 * x**9
27 |
28 |       XCO(i)   = 0.
29 |
30 |       XO2(i)   = 0.20927 - 0.082682 * x + 0.71791 * x**2 - 2.4765 * x**3&
31 |      &         + 4.1997 * x**4 - 3.71 * x**5 + 1.5718 * x**6 -          &
32 |      &           0.2514 * x**7
33 |
34 |       XN2(i)   = 0.79261 + 0.0040393 * x - 0.033259 * x**2 +            &
35 |      &           0.098691 * x**3 - 0.11162 * x**4 + 0.031811 * x**5
36 |
37 |       XSOOT(i) = 0.
38 |
39 |       if ( XH2O(i).lt.0. ) then
40 |         XH2O(i) = 0.
41 |       endif
42 |       if ( XCO2(i).lt.0. ) then
43 |         XCO2(i) = 0.
44 |       endif
45 |       if ( XO2 (i).lt.0. ) then
46 |         XO2 (i) = 0.
47 |       endif
48 |       if ( XN2 (i).lt.0. ) then
49 |         XN2 (i) = 0.
50 |       endif
51 |       if ( tempe(i).lt.300. ) then
52 |         tempe(i) = 300.
53 |       endif
```