SUNCAT-Center/catmap

how to mapped the production rate of one product?

zc-belief opened this issue · 1 comments

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.

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)`