Source code for grid_plot_elem_angles
#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""Plot a nice histogram of element angles to plot_element_angles.jpg. The
script expects an elem.dat file in the present working directory.
Examples
--------
>>> ls -1
elec.dat
elem.dat
>>> grid_plot_elem_angles
>>> ls -1
elec.dat
elem.dat
plot_element_angles.jpg
"""
import os
import crtomo.grid as CRGrid
[docs]
def main():
if not os.path.isfile('elem.dat'):
raise Exception('elem.dat not found!')
grid = CRGrid.crt_grid()
grid.load_elem_file('elem.dat')
fig, ax = grid.analyze_internal_angles(return_plot=True)
fig.savefig('plot_element_angles.jpg', dpi=300)