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_1.ipynb and run it in you own notebook, or see it directly on nbviewer here: http://nbviewer.ipython.org/url/sites.google.com/site/pycloudy/examples/example-2/Using_pyCloudy_2.ipynb