Producing Figure 3b from OzBench et al 2008

This script produces Figure 3b from the OzBench et al, 2008 benchmark using a glucifer store saved in OzBenchEtAl-2008.ipynb.

References

  1. OzBench, M.; Regenauer-Lieb, K.; Stegman, D. R.; Morra, G.; Farrington, R.; Hale, A.; May, D. A.; Freeman, J.; Bourgouin, L.; Mühlhaus, H. & Moresi, L. A model comparison study of large-scale mantle-lithosphere dynamics driven by subduction. Physics of the Earth and Planetary Interiors, 2008, 171, 224-234. OzBench, 2008
In [1]:
import os
outputPath = os.path.join(os.path.abspath("."),"output/")
if not os.path.exists(outputPath):
    os.makedirs(outputPath)
In [2]:
import glucifer
saved = glucifer.Viewer(outputPath+'ozBench', quality=3)

Create image for each saved timestep

In [3]:
#Apply the camera script
saved.commands(['translate y 1', 'rotate y 60', 'rotate x 30', 'zoom 0.1'])

#Write an image for each time step
for step in saved.timesteps():
    saved.timestep(step)
    saved["title"] = "Timestep ##"
    saved.image(outputPath + "ozBenchParticle_"+ str(step).zfill(4))
In [4]:
#Display the final image
saved.show()
In [ ]: