MTEX 4.0.6 (show documentation) Import pole figure data Import EBSD data Import ODF data Uninstall MTEX >> Multi_phases_seismic_anisotropy_harmonic_Hill_MTEX4 ebsd = EBSD (show methods, plot) Phase Orientations Mineral Color Symmetry Crystal reference frame 0 58485 (24%) notIndexed 1 152345 (62%) Forsterite dark green mmm 2 26058 (11%) Enstatite light red mmm 3 9064 (3.7%) Diopside dark blue 12/m1 X||a*, Y||b*, Z||c Properties: bands, bc, bs, error, mad, x, y Scan unit : um Crystal Symmetry (CS) : names of variables Variable name = Forsterite_CS Variable name = Enstatite_CS Variable name = Diopside_CS Use mineral_CS to defined your CS in m-file odf_olivine = ODF (show methods, plot) crystal symmetry : Forsterite (mmm) specimen symmetry: 1 Portion specified by Fourier coefficients: degree: 28 weight: 1 pfs_olivine_uvw = Miller (show methods, plot) size: 1 x 3 mineral: Forsterite (mmm) u 1 0 0 v 0 1 0 w 0 0 1 ebsd_corrected = EBSD (show methods, plot) Phase Orientations Mineral Color Symmetry Crystal reference frame 0 58485 (24%) notIndexed 1 152345 (62%) Forsterite dark green mmm 2 26058 (11%) Enstatite light red mmm 3 9064 (3.7%) Diopside dark blue 12/m1 X||a*, Y||b*, Z||c Properties: bands, bc, bs, error, mad, x, y Scan unit : um Area fractions for all phases in the EBSD object # Phase Points Area fraction 1 notIndexed 58485 0.2378 2 Forsterite 152345 0.6194 3 Enstatite 26058 0.1059 4 Diopside 9064 0.0369 Total area fraction = 1.0000 Percentage indexed points in map = 76.2 Area fractions based on minerals with EBSD indexed points only # Phase Points Area fraction 2 Forsterite 152345 0.8126 3 Enstatite 26058 0.1390 4 Diopside 9064 0.0483 Total area fraction = 1.0000 ********************************* Revised Area fraction calculation ********************************* You may use the volume fractions determined by the script from your *.ctf EBSD file. Alternatively you can enter area fractions determined by another method or remove some phases (e.g. optical microscope, electron probe etc.) MAKE SURE YOU DO NOT HAVE THE Si CALIBRATION STANDARD IN YOUR LIST OF INDEXED PHASES - Silicon Points > 0 IF SO REMOVE IT BY ANSWERING NO below and set Si area fraction to 0.0 Use the volume fractions given by the script from your EBSD file(Y/N) ?Y Area fractions based on minerals with EBSD indexed points only # Phase Points Area fraction 2 Forsterite 152345 0.8126 3 Enstatite 26058 0.1390 4 Diopside 9064 0.0483 Total area fraction = 1.0000 Area fractions : names of variables Variable name = af_Forsterite Area Fraction = 0.8126 Variable name = af_Enstatite Area Fraction = 0.1390 Variable name = af_Diopside Area Fraction = 0.0483 The af_mineral_name above will be used for area fractions in this m-file CS_Tensor_olivine = crystalSymmetry (show methods, plot) mineral : Olivine symmetry: mmm a, b, c : 4.8, 10, 6 C_olivine = elastic stiffness tensor (show methods, plot) unit : GPa density: 3.355 rank : 4 (3 x 3 x 3 x 3) mineral: Olivine (mmm) tensor in Voigt matrix representation: 320.5 68.2 71.6 0 0 0 68.2 196.5 76.8 0 0 0 71.6 76.8 233.5 0 0 0 0 0 0 64 0 0 0 0 0 0 77 0 0 0 0 0 0 78.7 Elastic tensor (GPa): Single crystal Olivine 320.50 68.15 71.60 0.00 0.00 0.00 68.15 196.50 76.80 0.00 0.00 0.00 71.60 76.80 233.50 0.00 0.00 0.00 0.00 0.00 0.00 64.00 0.00 0.00 0.00 0.00 0.00 0.00 77.00 0.00 0.00 0.00 0.00 0.00 0.00 78.70 Density (g/cm3) = 3.3550 Anisotropy Parameters P-wave anisotropy (%) = 24.34 S-wave maximum anisotropy (%) = 17.96 S-wave minimum anisotropy (%) = 0.09 dVs=Vs1-Vs2 maximum (km/s) = 0.90 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 13.04 S2-wave anisotropy (%) = 10.10 Vp/Vs1 anisotropy (%) = 25.59 Vp/Vs2 anisotropy (%) = 16.47 Plot has successfully terminated Elastic tensor (GPa): Single crystal Enstatite 236.90 79.60 63.20 0.00 0.00 0.00 79.60 180.50 56.80 0.00 0.00 0.00 63.20 56.80 230.40 0.00 0.00 0.00 0.00 0.00 0.00 84.30 0.00 0.00 0.00 0.00 0.00 0.00 79.40 0.00 0.00 0.00 0.00 0.00 0.00 80.10 Density (g/cm3) = 3.3060 Anisotropy Parameters P-wave anisotropy (%) = 13.60 S-wave maximum anisotropy (%) = 13.18 S-wave minimum anisotropy (%) = 0.01 dVs=Vs1-Vs2 maximum (km/s) = 0.62 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 3.70 S2-wave anisotropy (%) = 13.85 Vp/Vs1 anisotropy (%) = 16.53 Vp/Vs2 anisotropy (%) = 22.97 Plot has successfully terminated Elastic tensor (GPa): Single crystal Diopside 228.10 78.80 70.20 0.00 7.90 0.00 78.80 181.10 61.10 0.00 5.90 0.00 70.20 61.10 245.40 0.00 39.70 0.00 0.00 0.00 0.00 78.90 0.00 6.40 7.90 5.90 39.70 0.00 68.20 0.00 0.00 0.00 0.00 6.40 0.00 78.10 Density (g/cm3) = 3.2860 Anisotropy Parameters P-wave anisotropy (%) = 25.51 S-wave maximum anisotropy (%) = 22.54 S-wave minimum anisotropy (%) = 0.09 dVs=Vs1-Vs2 maximum (km/s) = 1.06 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 17.32 S2-wave anisotropy (%) = 20.68 Vp/Vs1 anisotropy (%) = 29.80 Vp/Vs2 anisotropy (%) = 36.70 Plot has successfully terminated ref_vec_1 = vector3d (show methods, plot) size: 1 x 1 x y z 1 0 0 ref_vec_2 = vector3d (show methods, plot) size: 1 x 1 x y z 0 1 0 ref_vec_3 = vector3d (show methods, plot) size: 1 x 1 x y z 0 0 -1 100% Olivine : Voigt, Hill and Reuss averages C_Voigt_Ol = elastic stiffness tensor (show methods, plot) unit : GPa density: 3.355 rank : 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 204.61 77.05 76.83 0.28 6.59 -2.99 77.05 247.35 81.57 5.85 0.83 -2.99 76.83 81.57 260.72 5.85 6.59 -0.88 0.28 5.85 5.85 86.23 -1.06 2.88 6.59 0.83 6.59 -1.06 77.92 2.79 -2.99 -2.99 -0.88 2.88 2.79 74.46 C_Reuss_Ol = elastic stiffness tensor (show methods, plot) unit : GPa density: 3.355 rank : 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 204.29 76.25 75.68 -0.15 4.79 -2.39 76.25 240.25 79.8 5.05 0.25 -2.7 75.68 79.8 251.59 5.28 5.71 -0.61 -0.15 5.05 5.28 83.19 -0.89 2.52 4.79 0.25 5.71 -0.89 75.32 2.18 -2.39 -2.7 -0.61 2.52 2.18 72.48 C_Hill_Ol = elastic stiffness tensor (show methods, plot) unit : GPa density: 3.355 rank : 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 204.45 76.65 76.26 0.07 5.69 -2.69 76.65 243.8 80.68 5.45 0.54 -2.85 76.26 80.68 256.15 5.57 6.15 -0.74 0.07 5.45 5.57 84.71 -0.97 2.7 5.69 0.54 6.15 -0.97 76.62 2.48 -2.69 -2.85 -0.74 2.7 2.48 73.47 Elastic tensor (GPa): 100%% Olivine : Hill average 204.45 76.65 76.26 0.07 5.69 -2.69 76.65 243.80 80.68 5.45 0.54 -2.85 76.26 80.68 256.15 5.57 6.15 -0.74 0.07 5.45 5.57 84.71 -0.97 2.70 5.69 0.54 6.15 -0.97 76.62 2.48 -2.69 -2.85 -0.74 2.70 2.48 73.47 Density (g/cm3) = 3.3550 Anisotropy Parameters P-wave anisotropy (%) = 13.66 S-wave maximum anisotropy (%) = 8.80 S-wave minimum anisotropy (%) = 0.06 dVs=Vs1-Vs2 maximum (km/s) = 0.43 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 5.61 S2-wave anisotropy (%) = 3.99 Vp/Vs1 anisotropy (%) = 8.60 Vp/Vs2 anisotropy (%) = 9.96 Plot has successfully terminated 100% Enstatite : Voigt, Hill and Reuss averages C_Voigt_Opx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 219.2 64.95 64.72 -1.17 0.16 -0.42 64.95 219.96 62.05 0.05 0.55 -0.42 64.72 62.05 224.4 0.05 0.16 -0.39 -1.17 0.05 0.05 80.06 -0.02 -0.19 0.16 0.55 0.16 -0.02 78.54 0.61 -0.42 -0.42 -0.39 -0.19 0.61 77.31 C_Reuss_Opx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 215.1 64 63.98 -1.25 0.17 -0.49 64 216.33 61.13 0.13 0.59 -0.49 63.98 61.13 222.01 0.14 0.16 -0.44 -1.25 0.13 0.14 79.01 -0.03 -0.21 0.17 0.59 0.16 -0.03 77.24 0.68 -0.49 -0.49 -0.44 -0.21 0.68 75.84 C_Hill_Opx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 217.15 64.48 64.35 -1.21 0.17 -0.46 64.48 218.14 61.59 0.09 0.57 -0.46 64.35 61.59 223.21 0.1 0.16 -0.41 -1.21 0.09 0.1 79.54 -0.03 -0.2 0.17 0.57 0.16 -0.03 77.89 0.65 -0.46 -0.46 -0.41 -0.2 0.65 76.58 Elastic tensor (GPa): 100%% Enstatite : Hill average 217.15 64.48 64.35 -1.21 0.17 -0.46 64.48 218.14 61.59 0.09 0.57 -0.46 64.35 61.59 223.21 0.10 0.16 -0.41 -1.21 0.09 0.10 79.54 -0.03 -0.20 0.17 0.57 0.16 -0.03 77.89 0.65 -0.46 -0.46 -0.41 -0.20 0.65 76.58 Density (g/cm3) = 3.3060 Anisotropy Parameters P-wave anisotropy (%) = 1.51 S-wave maximum anisotropy (%) = 2.08 S-wave minimum anisotropy (%) = 0.01 dVs=Vs1-Vs2 maximum (km/s) = 0.10 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 0.89 S2-wave anisotropy (%) = 1.21 Vp/Vs1 anisotropy (%) = 1.52 Vp/Vs2 anisotropy (%) = 1.22 Warning: print2array generating a 30.6M pixel image. This could be slow and might also cause memory problems. > In print2array at 62 In export_fig at 334 In saveFigure at 54 In Func_Seismic_multiplot_gfs_MTEX4 at 400 In MTEX4_multi_phases_seismic_anisotropy_harmonic_Hill_MTEX4 at 722 Plot has successfully terminated 100% Diopside : Voigt, Hill and Reuss averages K_BW4 = DirichletKernel (show methods, plot) bandwidth: 4 halfwidth: 29° diopside_ODF_Dirichlet = ODF (show methods, plot) crystal symmetry : Diopside (12/m1, X||a*, Y||b*, Z||c) specimen symmetry: 1 Portion specified by Fourier coefficients: degree: 10 weight: 1 C_Voigt_Cpx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 215.46 70.45 69.37 -0.45 0.95 0.8 70.45 226.39 69.1 0.02 0.38 0.8 69.37 69.1 215.11 0.02 0.95 0.4 -0.45 0.02 0.02 75.83 0.2 0.28 0.95 0.38 0.95 0.2 72.96 0.24 0.8 0.8 0.4 0.28 0.24 75.24 C_Reuss_Cpx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 206.78 67.2 66.31 -0.71 0.86 1.16 67.2 212.7 65.19 -0.45 0.45 1.18 66.31 65.19 204.55 -0.44 0.85 0.6 -0.71 -0.45 -0.44 71.69 0.27 0.21 0.86 0.45 0.85 0.27 69.68 0.14 1.16 1.18 0.6 0.21 0.14 71.26 C_Hill_Cpx = tensor (show methods, plot) rank: 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 211.12 68.82 67.84 -0.58 0.9 0.98 68.82 219.54 67.14 -0.22 0.41 0.99 67.84 67.14 209.83 -0.21 0.9 0.5 -0.58 -0.22 -0.21 73.76 0.24 0.25 0.9 0.41 0.9 0.24 71.32 0.19 0.98 0.99 0.5 0.25 0.19 73.25 Elastic tensor (GPa): 100%% Diopside : Hill average 211.12 68.82 67.84 -0.58 0.90 0.98 68.82 219.54 67.14 -0.22 0.41 0.99 67.84 67.14 209.83 -0.21 0.90 0.50 -0.58 -0.22 -0.21 73.76 0.24 0.25 0.90 0.41 0.90 0.24 71.32 0.19 0.98 0.99 0.50 0.25 0.19 73.25 Density (g/cm3) = 3.2860 Anisotropy Parameters P-wave anisotropy (%) = 2.71 S-wave maximum anisotropy (%) = 1.79 S-wave minimum anisotropy (%) = 0.01 dVs=Vs1-Vs2 maximum (km/s) = 0.08 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 0.49 S2-wave anisotropy (%) = 1.33 Vp/Vs1 anisotropy (%) = 2.55 Vp/Vs2 anisotropy (%) = 1.50 Warning: print2array generating a 30.6M pixel image. This could be slow and might also cause memory problems. > In print2array at 62 In export_fig at 334 In saveFigure at 54 In Func_Seismic_multiplot_gfs_MTEX4 at 400 In MTEX4_multi_phases_seismic_anisotropy_harmonic_Hill_MTEX4 at 745 Plot has successfully terminated Aggregate density (g/cm3) = 3.3449 Aggregate (Olivine,Enstatite,Diopside) : Hill average C_Hill_Ag = elastic stiffness tensor (show methods, plot) unit : GPa density: 3.355 rank : 4 (3 x 3 x 3 x 3) tensor in Voigt matrix representation: 206.54 74.58 74.19 -0.14 4.69 -2.2 74.58 239.06 77.38 4.43 0.54 -2.33 74.19 77.38 249.33 4.53 5.06 -0.64 -0.14 4.43 4.53 83.46 -0.78 2.18 4.69 0.54 5.06 -0.78 76.54 2.12 -2.2 -2.33 -0.64 2.18 2.12 73.89 Elastic tensor (GPa): Aggregate (Olivine,Enstatite,Diopside) : Hill average 206.54 74.58 74.19 -0.14 4.69 -2.20 74.58 239.06 77.38 4.43 0.54 -2.33 74.19 77.38 249.33 4.53 5.06 -0.64 -0.14 4.43 4.53 83.46 -0.78 2.18 4.69 0.54 5.06 -0.78 76.54 2.12 -2.20 -2.33 -0.64 2.18 2.12 73.89 Density (g/cm3) = 3.3449 Anisotropy Parameters P-wave anisotropy (%) = 11.38 S-wave maximum anisotropy (%) = 7.46 S-wave minimum anisotropy (%) = 0.04 dVs=Vs1-Vs2 maximum (km/s) = 0.36 dVs=Vs1-Vs2 minimum (km/s) = 0.00 S1-wave anisotropy (%) = 4.61 S2-wave anisotropy (%) = 3.32 Vp/Vs1 anisotropy (%) = 7.09 Vp/Vs2 anisotropy (%) = 8.23 Plot has successfully terminated Script has successfully terminated >>