Figure 5.23:

Dimensionless CO_2 production rate versus dimensionless gas-phase CO and O_2 concentrations.

Code for Figure 5.23

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# Converted from CO3d.m
import numpy as np

xisosamples = 30; yisosamples = 30
x = np.logspace(-2, 2, xisosamples)
y = np.logspace(-2, 3, yisosamples)
xx, yy = np.meshgrid(x, y)
zz = xx * np.sqrt(yy) / (1 + xx + np.sqrt(yy))**2

with open('CO3d.dat', 'w') as f:
    for j in range(yisosamples):
        for i in range(xisosamples):
            f.write(f'{x[i]:f} {y[j]:f} {zz[j, i]:f}\n')
        if j < yisosamples - 1:
            f.write('\n')