In this example, we run 20 models (using 3 processors) changing the density and the O/H for each model. We read the results and plot [OIII]5007/4363 vs. [OII] 3726/29, using symbol colors and sizes related to O/H and the density resp.
The plot will look like this:
You can see step by step the workflow of this small script using ipython notebook. You can download the Using_pyCloudy_2.ipynb and run it in you own notebook, or see it directly on nbviewer here: https://github.com/Morisset/pyCloudy/blob/master/pyCloudy/docs/Using_pyCloudy_2.ipynb