how to mapped the production rate of one product?
zc-belief opened this issue · 1 comments
zc-belief commented
Because there are two products in my reaction network, if i just want to map one product, please guide me how to set it? Thanks a lot.
gusmaogabriels commented
Here follows a very simple example where you plot only the production rates of 'CH3OH_g'
and 'CH4_g'
, which are species that should exist in you microkinetic model. The CatMAP model object (model
) instantiation was suppressed in this snippet.
from catmap import analyze
from matplotlib import pyplot as plt
labels = ['CH3OH_g','CH4_g']
vm = analyze.VectorMap(model)
vm.include_empty = False
vm.include_labels = labels;
vm.plot_variable = 'production_rate' #tell the model which output to plot
vm.log_scale = True #rates should be plotted on a log-scale
vm.min = 1e-50#minimum rate to plot
vm.max = 1e8 #maximum rate to plot
vm.threshold = 1e-50 #anything below this is considered to be 0
vm.subplots_adjust_kwargs = {'left':0.2,'right':0.95,'bottom':0.15}
fig = plt.figure(figsize=[6.,6./2],dpi=200)
ax_list=[fig.add_subplot('12{}'.format(_)) for _ in range(2)]
ax_list = [ax_list[0],ax_list[1]]
vm.plot(ax_list=ax_list)`