Adding column density per zone

[Morisset]

Following a comment by gary on the Cloudy group, one could do this easily in PyNeb:

"You can save the density of the species in each zone and use the SAVE DR
command to get the thickness of each zone. If you then write a little
script to multiply those two numbers, you get the column density per zone."