1 | !     ==================================================================!
   2 | !     File used for generating an initial solution in PRISSMA           !
   3 | !     Values are given at the i-th node "_n(i)"                         !
   4 | !     Auziliary real type available variables: r, z_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)     = 800.+1200.*(1.-r/R_ref)*(x(i)/L_ref)
  12 | 
  13 |       P_n(i)     = 1.
  14 |       XH2O_n(i)  = 0.05*(1.-2. *(z(i)/L_ref-0.5)**2)*(2.-r/R_ref)
  15 |       XCO2_n(i)  = 0.04*(1.-3. *(z(i)/L_ref-0.5)**2)*(2.5-r/R_ref)
  16 |       XCO_n(i)   = 0.
  17 |       XO2_n(i)   = 0.2*(1.-XH2O_n(i)-XCO2_n(i)-XCO_n(i))
  18 |       XN2_n(i)   = 0.8*(1.-XH2O_n(i)-XCO2_n(i)-XCO_n(i))
  19 |       XSOOT_n(i) = 1.e-7
  20 | 
  21 |       kscattering_n(i) = 0.
  22 |       kabsorption_n(i) = 0.1
  23 |