How to extract stress within a certain bounding geometry in an Abaqus ODB using scripting?

Question

I wish to be able to extract stress (S11) and its corresponding element's z-coordinates only in a certain bounding cylinder. My ODB files are relatively big and it would be quite helpful if I could limit the number of elements I have to compute. I need to use scripting and avoid viewport manipulation as I've a large number of ODBs to process.

I have attached a schematic of an example ODB.

Please click here for a schematic of an example ODB.