Simulate the volumetric cooling of a surface or a mesh.

Input parameters:

1. Sh → Surface or mesh to analyse.
2. Un → Optional voxel accuracy. Smaller values drastically increase computation time.
3. Gr → Normalised temperature gradient.
4. Fn → Optional function for post-processing subsurface temperature gradient values.
5. Sc → Optional planar sections through the shape.

Output parameters:

1. Sf → Cooling gradient on shape surface.
2. Sc → Sections with normalised temperature distributions.