1 | !     ==================================================================!
   2 | !     File used for generating an initial solution in PRISSMA           !
   3 | !     Values are given at the i-th node "_n(i)"                         !
   4 | !     Auxiliary real type available variables: r, x_ref, y_ref, z_ref   !
   5 | !                                              R_ref, L_ref, T_ref      !
   6 | !     ==================================================================!
   7 |       R_ref = 0.3
   8 |       L_ref = 1.2
   9 |       r = SQRT(z(i)*z(i)+y(i)*y(i))
  10 | 
  11 |       T_n(i)     = 0.97*800.+1.17*1200.*(1.-r/R_ref)*(x(i)/L_ref)
  12 | !     0,97 and 1,17 coeff to obtain the good profile on the mesh
  13 |  
  14 |       P_n(i)     = 1.
  15 |       XH2O_n(i)  = 0.05*(1.-2. *(x(i)/L_ref-0.5)**2)*(2.-r/R_ref)
  16 |       XCO2_n(i)  = 0.04*(1.-3. *(x(i)/L_ref-0.5)**2)*(2.5-r/R_ref)
  17 |       XCO_n(i)   = 0.
  18 |       XO2_n(i)   = 0.2*(1.-XH2O_n(i)-XCO2_n(i)-XCO_n(i))
  19 |       XN2_n(i)   = 0.8*(1.-XH2O_n(i)-XCO2_n(i)-XCO_n(i))
  20 |       XSOOT_n(i) = 0.   
  21 | 
  22 |       kscattering_n(i) = 0.
  23 |       kabsorption_n(i) = 0.1
  24 |