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