In this example, we run 20 models (using 8 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.
You can see step by step the workflow of this small script using ipython notebook. You can see it directly on here: https://github.com/Morisset/pyCloudy/blob/master/pyCloudy/docs/Using_pyCloudy_2.ipynb