Hi Eamonn,
The current version of pyCloudy is axy-simmetric, it only considers morphologies where a dependency in theta (equatorial angle) is defined, no phy (longitude) variations.
I cannot work on this right now, but the idea would be to add an effective Delaunay triangulation method. A draft version of this is already implemented, in pyCloudy.3d.model_3d._get_interp_tri, but I never finished it...
In any case, you will need to run quite a lot of models to reach a decent spatial resolution. I'm actually doing this in an axy-simmetric model and use 180 1D-models to feed the 3D model.
We can try to work on this in September, contact me at that time.
Christophe