gbcd_YAG.txt msym 3 # 32 ax(1) ax(2) ax(3) ang 1.0 0.0 0.0 5 1.0 0.0 0.0 15 1.0 0.0 0.0 25 1.0 0.0 0.0 35 1.0 0.0 0.0 45 1.0 1.0 0.0 10 1.0 1.0 0.0 20 1.0 1.0 0.0 30 1.0 1.0 0.0 40 1.0 1.0 0.0 50 1.0 1.0 0.0 60 1.0 1.0 1.0 10 1.0 1.0 1.0 20 1.0 1.0 1.0 30 1.0 1.0 1.0 40 1.0 1.0 1.0 50 1.0 1.0 1.0 60 3.0 2.0 1.0 10 3.0 2.0 1.0 20 3.0 2.0 1.0 30 3.0 2.0 1.0 40 3.0 2.0 1.0 50 3.0 2.0 1.0 60 1.0 0.0 0.0 36.9 1.0 1.0 1.0 38.2 1.0 1.0 0.0 38.9 1.0 1.0 0.0 50.5 1.0 0.0 0.0 22.6 1.0 1.0 1.0 27.8 2.0 1.0 0.0 48.2 1.0 1.0 0.0 31.6 2.0 1.0 0.0 35.4 1.0 1.0 1.0 60 1.0 1.0 0.0 31.58 1.0 1.0 0.0 38.94 1.0 1.0 1.0 60 1.0 0.02 0.0 90 1.0 1.732 0.0 90 -1.0 1.732 0.0 90 -1.0 0.02 0.0 90 -1.0 -1.732 0.0 90 1.0 -1.732 0.0 90 0.00 -1.0 0.0 90 -1.732 -1.0 0.0 90 -1.732 1.0 0.0 90 0.00 1.0 0.0 90 1.732 1.0 0.0 90 1.732 -1.0 0.0 90 0.02 -1.0 0.0 90 -1.732 -1.0 0.0 90 -1.732 1.0 0.0 90 0.02 1.0 0.0 90 1.732 1.0 0.0 90 1.732 -1.0 0.0 90 0.0 0.0 1.0 30.0 0.0 0.0 1.0 60.0 0.0 0.0 1.0 90.0 0.0 0.0 1.0 120.0 0.0 0.0 1.0 150.0 0.0 0.0 1.0 180.0 1.732 1.0 0.0 60.0 0.0 0.0 1.0 30.0 1.0 0.0 0.0 90.0 The first line is the name of the input file. This file should be in this directory, and it should have the following properties: Each line contains a floating point number. The units of the number are in MRD. The first value corresponds to the GBCD at phi_1 = 0, cap_phi= 0, phi_2 = 0, theta = 0, Phi = 0. Phi is then incremented, followed by theta, etc... phi_1, phi_2, and Phi are incremented in 10 degree intervals between zero and 90. cap_phi and theta are discretized into nine equal intervals of cos(theta). Files produced by calc_gbcd are in this format. If Cubic, msym = 3 If Hexagonal, msym = 2 If Tetragonal, msym = 1 If Trigonal, msym = 4 The axis values (ax(i)) are in vectors (x,y,z), in the cartesian (not the crystal) reference frame and the angle is in degrees. In the plot. produced, the [100] direction is in the plane, pointing to the right. The [010] direction is in the plane, pointing up. For each plot the program writes the max MRD and the MRD for the misorientation (by averaging the values at all plane orientations) to the terminal window. The output of this program are two sets of files: 1. result_5d_n.dat where n is the integer order from the axis angle list. The plots are in the format for 3D field. 2. 5d_orig_degreesn.dat where n is the integer order from the axis angle list. The plots are in the format for GMT. If GMT is installed, you can use the script "Draw_stereograms" to make plots. Here are some lists a axis angle pair that are convenient for copy and paste. 1.0 1.0 0.0 20 1.0 1.0 0.0 30 1.0 1.0 0.0 31.58 1.0 1.0 0.0 38.94 1.0 1.0 0.0 40 1.0 1.0 0.0 50 0.0 0.0 1.0 43.17 1.0 1.0 1.0 60 1.0 0.0 0.0 5.0 1.0 0.0 0.0 10.0 1.0 0.0 0.0 15.0 1.0 0.0 0.0 25.0 1.0 0.0 0.0 35.0 1.0 0.0 0.0 45.0 0.02 -1.0 0.0 30 -1.732 -1.0 0.0 30 -1.732 1.0 0.0 30 0.02 1.0 0.0 30 1.732 1.0 0.0 30 1.732 -1.0 0.0 30 sigma 19c 1.0 0.02 0.0 86.98 low angle 0001 series 0.0 0.0 1.0 5.0 0.0 0.0 1.0 7.0 0.0 0.0 1.0 9.0 0.0 0.0 1.0 11.0 0.0 0.0 1.0 13.0 0.0 0.0 1.0 15.0 high index peak in 60-65 range 11.5 1.0 3.72 61.0 Peak in 40-45 range 5.5 1.0 5.60 42.5 sigma 19b 5.0 0.02 1.867 65.1 sigma 19a 0.0 0.0 1.0 13.17 sigma 13b 1.0 0.02 0.0 57.4 1.0 1.732 0.0 57.4 -1.0 1.732 0.0 57.4 -1.0 0.02 0.0 57.4 -1.0 -1.732 0.0 57.4 1.0 -1.732 0.0 57.4 1.0 0.02 0.0 84.78 1.0 1.732 0.0 84.78 -1.0 1.732 0.0 84.78 -1.0 0.02 0.0 84.78 -1.0 -1.732 0.0 84.78 1.0 -1.732 0.0 84.78 0.02 -1.0 0.0 86.98 -1.732 -1.0 0.0 86.98 -1.732 1.0 0.0 86.98 0.02 1.0 0.0 86.98 1.732 1.0 0.0 86.98 1.732 -1.0 0.0 86.98 0.0 0.0 1.0 5 0.0 0.0 1.0 15 0.0 0.0 1.0 25 0.0 0.0 1.0 35 0.0 0.0 1.0 45 0.0 0.0 1.0 55 This is the [10-10] set 0.02 -1.0 0.0 60 -1.732 -1.0 0.0 60 -1.732 1.0 0.0 60 0.02 1.0 0.0 60 1.732 1.0 0.0 60 1.732 -1.0 0.0 60 This is the [2-1-10] set 1.0 0.02 0.0 30 1.0 1.732 0.0 30 -1.0 1.732 0.0 30 -1.0 0.02 0.0 30 -1.0 -1.732 0.0 30 1.0 -1.732 0.0 30 1.0 1.0 0.0 31.58 1.0 1.0 0.0 38.94 0.0 -1.0 0.0 20 1.0 0.02 0.0 20 -1.0 0.02 0.0 20 0.02 -1.0 0.0 20 1.0 1.0 1.0 10.0 1.0 1.0 1.0 20.0 1.0 1.0 1.0 30.0 1.0 1.0 1.0 40.0 1.0 1.0 1.0 50.0 1.0 1.0 1.0 60.0 1.0 1.0 0.0 10.0 1.0 1.0 0.0 20.0 1.0 1.0 0.0 30.0 1.0 1.0 0.0 40.0 1.0 1.0 0.0 50.0 1.0 1.0 0.0 60.0 1.0 0.0 0.0 5.0 1.0 0.0 0.0 15.0 1.0 0.0 0.0 25.0 1.0 0.0 0.0 35.0 1.0 0.0 0.0 45.0 1.0 0.0 0.0 36.86 (sigma 5) 1.0 1.0 1.0 38.21 (sigma 7) 1.0 1.0 0.0 38.94 (sigma 9) 1.0 1.0 0.0 31.59 (sigma 27a) 2.0 1.0 0.0 35.43 (sigma 27b) 1.0 1.0 1.0 60.0 (sigma 3) 1.0 0.0 0.0 0.0 (low angle) 1.0 1.0 1.0 5.0 1.0 1.0 1.0 10.0 1.0 1.0 1.0 20.0 1.0 1.0 1.0 30.0 1.0 1.0 1.0 40.0 1.0 1.0 1.0 50.0 1.0 1.0 1.0 60.0 1.0 1.0 0.0 10.0 1.0 1.0 0.0 20.0 1.0 1.0 0.0 30.0 1.0 1.0 0.0 40.0 1.0 1.0 0.0 50.0 1.0 1.0 0.0 60.0 1.0 0.0 0.0 10.0 1.0 0.0 0.0 20.0 1.0 0.0 0.0 30.0 1.0 0.0 0.0 40.0 1.0 0.0 0.0 50.0 [maximum # of plots = 9]