diff --git a/src/2DPhaseDiagram_PlotScript.py b/src/2DPhaseDiagram_PlotScript.py
new file mode 100644
index 0000000000000000000000000000000000000000..0079b6995cc92bac5e6db3f10a004c745bc42e58
--- /dev/null
+++ b/src/2DPhaseDiagram_PlotScript.py
@@ -0,0 +1,910 @@
+## trace generated using paraview version 5.7.0
+#
+# To ensure correct image size when batch processing, please search
+# for and uncomment the line `# renderView*.ViewSize = [*,*]`
+
+#### import the simple module from the paraview
+from paraview.simple import *
+#### disable automatic camera reset on 'Show'
+paraview.simple._DisableFirstRenderCameraReset()
+
+# create a new 'XML Structured Grid Reader'
+phaseDiagram2DGamma_infinity55_4000spvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram2DGamma_infinity-5+5_4000sp.vts'])
+phaseDiagram2DGamma_infinity55_4000spvts.CellArrayStatus = []
+phaseDiagram2DGamma_infinity55_4000spvts.PointArrayStatus = ['Type', 'angles', 'curvature']
+
+# get active view
+renderView1 = GetActiveViewOrCreate('RenderView')
+# uncomment following to set a specific view size
+# renderView1.ViewSize = [1380, 547]
+
+# show data in view
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay = Show(phaseDiagram2DGamma_infinity55_4000spvts, renderView1)
+
+# get color transfer function/color map for 'Type'
+typeLUT = GetColorTransferFunction('Type')
+typeLUT.AutomaticRescaleRangeMode = "Grow and update on 'Apply'"
+typeLUT.InterpretValuesAsCategories = 0
+typeLUT.AnnotationsInitialized = 0
+typeLUT.ShowCategoricalColorsinDataRangeOnly = 0
+typeLUT.RescaleOnVisibilityChange = 0
+typeLUT.EnableOpacityMapping = 0
+typeLUT.RGBPoints = [1.0, 0.001462, 0.000466, 0.013866, 1.007844, 0.002258, 0.001295, 0.018331, 1.015686, 0.003279, 0.002305, 0.023708, 1.02353, 0.004512, 0.00349, 0.029965, 1.0313720000000002, 0.00595, 0.004843, 0.03713, 1.0392160000000001, 0.007588, 0.006356, 0.044973, 1.0470579999999998, 0.009426, 0.008022, 0.052844, 1.0549019999999998, 0.011465, 0.009828, 0.06075, 1.0627460000000002, 0.013708, 0.011771, 0.068667, 1.0705879999999999, 0.016156, 0.01384, 0.076603, 1.0784319999999998, 0.018815, 0.016026, 0.084584, 1.086274, 0.021692, 0.01832, 0.09261, 1.094118, 0.024792, 0.020715, 0.100676, 1.10196, 0.028123, 0.023201, 0.108787, 1.109804, 0.031696, 0.025765, 0.116965, 1.117648, 0.03552, 0.028397, 0.125209, 1.12549, 0.039608, 0.03109, 0.133515, 1.133334, 0.04383, 0.03383, 0.141886, 1.1411760000000002, 0.048062, 0.036607, 0.150327, 1.1490200000000002, 0.05232, 0.039407, 0.158841, 1.1568619999999998, 0.056615, 0.04216, 0.167446, 1.1647059999999998, 0.060949, 0.044794, 0.176129, 1.1725500000000002, 0.06533, 0.047318, 0.184892, 1.1803919999999999, 0.069764, 0.049726, 0.193735, 1.1882359999999998, 0.074257, 0.052017, 0.20266, 1.196078, 0.078815, 0.054184, 0.211667, 1.203922, 0.083446, 0.056225, 0.220755, 1.211764, 0.088155, 0.058133, 0.229922, 1.219608, 0.092949, 0.059904, 0.239164, 1.2274500000000002, 0.097833, 0.061531, 0.248477, 1.2352940000000001, 0.102815, 0.06301, 0.257854, 1.243138, 0.107899, 0.064335, 0.267289, 1.2509800000000002, 0.113094, 0.065492, 0.276784, 1.2588240000000002, 0.118405, 0.066479, 0.286321, 1.2666659999999998, 0.123833, 0.067295, 0.295879, 1.2745099999999998, 0.12938, 0.067935, 0.305443, 1.282352, 0.135053, 0.068391, 0.315, 1.290196, 0.140858, 0.068654, 0.324538, 1.2980399999999999, 0.146785, 0.068738, 0.334011, 1.305882, 0.152839, 0.068637, 0.343404, 1.313726, 0.159018, 0.068354, 0.352688, 1.321568, 0.165308, 0.067911, 0.361816, 1.329412, 0.171713, 0.067305, 0.370771, 1.3372540000000002, 0.178212, 0.066576, 0.379497, 1.3450980000000001, 0.184801, 0.065732, 0.387973, 1.352942, 0.19146, 0.064818, 0.396152, 1.3607839999999998, 0.198177, 0.063862, 0.404009, 1.3686280000000002, 0.204935, 0.062907, 0.411514, 1.3764699999999999, 0.211718, 0.061992, 0.418647, 1.3843139999999998, 0.218512, 0.061158, 0.425392, 1.392156, 0.225302, 0.060445, 0.431742, 1.4, 0.232077, 0.059889, 0.437695, 1.4078439999999999, 0.238826, 0.059517, 0.443256, 1.415686, 0.245543, 0.059352, 0.448436, 1.42353, 0.25222, 0.059415, 0.453248, 1.431372, 0.258857, 0.059706, 0.45771, 1.439216, 0.265447, 0.060237, 0.46184, 1.4470580000000002, 0.271994, 0.060994, 0.46566, 1.4549020000000001, 0.278493, 0.061978, 0.46919, 1.462746, 0.284951, 0.063168, 0.472451, 1.4705880000000002, 0.291366, 0.064553, 0.475462, 1.4784320000000002, 0.29774, 0.066117, 0.478243, 1.4862739999999999, 0.304081, 0.067835, 0.480812, 1.4941179999999998, 0.310382, 0.069702, 0.483186, 1.50196, 0.316654, 0.07169, 0.48538, 1.509804, 0.322899, 0.073782, 0.487408, 1.5176479999999999, 0.329114, 0.075972, 0.489287, 1.52549, 0.335308, 0.078236, 0.491024, 1.533334, 0.341482, 0.080564, 0.492631, 1.541176, 0.347636, 0.082946, 0.494121, 1.54902, 0.353773, 0.085373, 0.495501, 1.5568619999999997, 0.359898, 0.087831, 0.496778, 1.5647060000000002, 0.366012, 0.090314, 0.49796, 1.5725500000000001, 0.372116, 0.092816, 0.499053, 1.5803919999999998, 0.378211, 0.095332, 0.500067, 1.5882360000000002, 0.384299, 0.097855, 0.501002, 1.5960779999999999, 0.390384, 0.100379, 0.501864, 1.6039219999999998, 0.396467, 0.102902, 0.502658, 1.611764, 0.402548, 0.10542, 0.503386, 1.619608, 0.408629, 0.10793, 0.504052, 1.62745, 0.414709, 0.110431, 0.504662, 1.635294, 0.420791, 0.11292, 0.505215, 1.643138, 0.426877, 0.115395, 0.505714, 1.6509800000000001, 0.432967, 0.117855, 0.50616, 1.658824, 0.439062, 0.120298, 0.506555, 1.6666660000000002, 0.445163, 0.122724, 0.506901, 1.6745100000000002, 0.451271, 0.125132, 0.507198, 1.6823519999999998, 0.457386, 0.127522, 0.507448, 1.6901960000000003, 0.463508, 0.129893, 0.507652, 1.6980399999999998, 0.46964, 0.132245, 0.507809, 1.705882, 0.47578, 0.134577, 0.507921, 1.7137259999999999, 0.481929, 0.136891, 0.507989, 1.721568, 0.488088, 0.139186, 0.508011, 1.729412, 0.494258, 0.141462, 0.507988, 1.737254, 0.500438, 0.143719, 0.50792, 1.745098, 0.506629, 0.145958, 0.507806, 1.752942, 0.512831, 0.148179, 0.507648, 1.7607840000000001, 0.519045, 0.150383, 0.507443, 1.768628, 0.52527, 0.152569, 0.507192, 1.7764699999999998, 0.531507, 0.154739, 0.506895, 1.7843140000000002, 0.537755, 0.156894, 0.506551, 1.7921559999999999, 0.544015, 0.159033, 0.506159, 1.7999999999999998, 0.550287, 0.161158, 0.505719, 1.8078440000000002, 0.556571, 0.163269, 0.50523, 1.815686, 0.562866, 0.165368, 0.504692, 1.8235299999999999, 0.569172, 0.167454, 0.504105, 1.831372, 0.57549, 0.16953, 0.503466, 1.839216, 0.581819, 0.171596, 0.502777, 1.847058, 0.588158, 0.173652, 0.502035, 1.854902, 0.594508, 0.175701, 0.501241, 1.862746, 0.600868, 0.177743, 0.500394, 1.8705880000000001, 0.607238, 0.179779, 0.499492, 1.878432, 0.613617, 0.181811, 0.498536, 1.8862740000000002, 0.620005, 0.18384, 0.497524, 1.8941179999999997, 0.626401, 0.185867, 0.496456, 1.9019599999999999, 0.632805, 0.187893, 0.495332, 1.9098039999999998, 0.639216, 0.189921, 0.49415, 1.9176479999999998, 0.645633, 0.191952, 0.49291, 1.92549, 0.652056, 0.193986, 0.491611, 1.9333339999999999, 0.658483, 0.196027, 0.490253, 1.941176, 0.664915, 0.198075, 0.488836, 1.94902, 0.671349, 0.200133, 0.487358, 1.956862, 0.677786, 0.202203, 0.485819, 1.964706, 0.684224, 0.204286, 0.484219, 1.97255, 0.690661, 0.206384, 0.482558, 1.9803920000000002, 0.697098, 0.208501, 0.480835, 1.9882360000000001, 0.703532, 0.210638, 0.479049, 1.9960779999999998, 0.709962, 0.212797, 0.477201, 2.003922, 0.716387, 0.214982, 0.47529, 2.0117640000000003, 0.722805, 0.217194, 0.473316, 2.019608, 0.729216, 0.219437, 0.471279, 2.02745, 0.735616, 0.221713, 0.46918, 2.035294, 0.742004, 0.224025, 0.467018, 2.043138, 0.748378, 0.226377, 0.464794, 2.05098, 0.754737, 0.228772, 0.462509, 2.058824, 0.761077, 0.231214, 0.460162, 2.0666659999999997, 0.767398, 0.233705, 0.457755, 2.07451, 0.773695, 0.236249, 0.455289, 2.082352, 0.779968, 0.238851, 0.452765, 2.0901959999999997, 0.786212, 0.241514, 0.450184, 2.09804, 0.792427, 0.244242, 0.447543, 2.1058820000000003, 0.798608, 0.24704, 0.444848, 2.1137259999999998, 0.804752, 0.249911, 0.442102, 2.121568, 0.810855, 0.252861, 0.439305, 2.1294120000000003, 0.816914, 0.255895, 0.436461, 2.137254, 0.822926, 0.259016, 0.433573, 2.145098, 0.828886, 0.262229, 0.430644, 2.152942, 0.834791, 0.26554, 0.427671, 2.160784, 0.840636, 0.268953, 0.424666, 2.168628, 0.846416, 0.272473, 0.421631, 2.17647, 0.852126, 0.276106, 0.418573, 2.184314, 0.857763, 0.279857, 0.415496, 2.1921559999999998, 0.86332, 0.283729, 0.412403, 2.2, 0.868793, 0.287728, 0.409303, 2.2078439999999997, 0.874176, 0.291859, 0.406205, 2.215686, 0.879464, 0.296125, 0.403118, 2.2235300000000002, 0.884651, 0.30053, 0.400047, 2.231372, 0.889731, 0.305079, 0.397002, 2.239216, 0.8947, 0.309773, 0.393995, 2.247058, 0.899552, 0.314616, 0.391037, 2.254902, 0.904281, 0.31961, 0.388137, 2.262746, 0.908884, 0.324755, 0.385308, 2.270588, 0.913354, 0.330052, 0.382563, 2.278432, 0.917689, 0.3355, 0.379915, 2.2862739999999997, 0.921884, 0.341098, 0.377376, 2.294118, 0.925937, 0.346844, 0.374959, 2.3019600000000002, 0.929845, 0.352734, 0.372677, 2.3098039999999997, 0.933606, 0.358764, 0.370541, 2.317648, 0.937221, 0.364929, 0.368567, 2.3254900000000003, 0.940687, 0.371224, 0.366762, 2.333334, 0.944006, 0.377643, 0.365136, 2.341176, 0.94718, 0.384178, 0.363701, 2.3490200000000003, 0.95021, 0.39082, 0.362468, 2.356862, 0.953099, 0.397563, 0.361438, 2.364706, 0.955849, 0.4044, 0.360619, 2.37255, 0.958464, 0.411324, 0.360014, 2.380392, 0.960949, 0.418323, 0.35963, 2.388236, 0.96331, 0.42539, 0.359469, 2.396078, 0.965549, 0.432519, 0.359529, 2.4039219999999997, 0.967671, 0.439703, 0.35981, 2.411764, 0.96968, 0.446936, 0.360311, 2.419608, 0.971582, 0.45421, 0.36103, 2.4274500000000003, 0.973381, 0.46152, 0.361965, 2.435294, 0.975082, 0.468861, 0.363111, 2.4431380000000003, 0.97669, 0.476226, 0.364466, 2.45098, 0.97821, 0.483612, 0.366025, 2.458824, 0.979645, 0.491014, 0.367783, 2.466666, 0.981, 0.498428, 0.369734, 2.47451, 0.982279, 0.505851, 0.371874, 2.4823519999999997, 0.983485, 0.51328, 0.374198, 2.490196, 0.984622, 0.520713, 0.376698, 2.49804, 0.985693, 0.528148, 0.379371, 2.5058819999999997, 0.9867, 0.535582, 0.38221, 2.513726, 0.987646, 0.543015, 0.38521, 2.5215680000000003, 0.988533, 0.550446, 0.388365, 2.5294119999999998, 0.989363, 0.557873, 0.391671, 2.537254, 0.990138, 0.565296, 0.395122, 2.5450980000000003, 0.990871, 0.572706, 0.398714, 2.552942, 0.991558, 0.580107, 0.402441, 2.560784, 0.992196, 0.587502, 0.406299, 2.568628, 0.992785, 0.594891, 0.410283, 2.57647, 0.993326, 0.602275, 0.41439, 2.584314, 0.993834, 0.609644, 0.418613, 2.592156, 0.994309, 0.616999, 0.42295, 2.6, 0.994738, 0.62435, 0.427397, 2.607844, 0.995122, 0.631696, 0.431951, 2.615686, 0.99548, 0.639027, 0.436607, 2.6235299999999997, 0.99581, 0.646344, 0.441361, 2.631372, 0.996096, 0.653659, 0.446213, 2.6392160000000002, 0.996341, 0.660969, 0.45116, 2.647058, 0.99658, 0.668256, 0.456192, 2.654902, 0.996775, 0.675541, 0.461314, 2.6627460000000003, 0.996925, 0.682828, 0.466526, 2.670588, 0.997077, 0.690088, 0.471811, 2.678432, 0.997186, 0.697349, 0.477182, 2.686274, 0.997254, 0.704611, 0.482635, 2.694118, 0.997325, 0.711848, 0.488154, 2.7019599999999997, 0.997351, 0.719089, 0.493755, 2.709804, 0.997351, 0.726324, 0.499428, 2.717648, 0.997341, 0.733545, 0.505167, 2.7254899999999997, 0.997285, 0.740772, 0.510983, 2.733334, 0.997228, 0.747981, 0.516859, 2.7411760000000003, 0.997138, 0.75519, 0.522806, 2.74902, 0.997019, 0.762398, 0.528821, 2.756862, 0.996898, 0.769591, 0.534892, 2.7647060000000003, 0.996727, 0.776795, 0.541039, 2.77255, 0.996571, 0.783977, 0.547233, 2.780392, 0.996369, 0.791167, 0.553499, 2.788236, 0.996162, 0.798348, 0.55982, 2.796078, 0.995932, 0.805527, 0.566202, 2.803922, 0.99568, 0.812706, 0.572645, 2.811764, 0.995424, 0.819875, 0.57914, 2.8196079999999997, 0.995131, 0.827052, 0.585701, 2.82745, 0.994851, 0.834213, 0.592307, 2.835294, 0.994524, 0.841387, 0.598983, 2.8431379999999997, 0.994222, 0.84854, 0.605696, 2.85098, 0.993866, 0.855711, 0.612482, 2.8588240000000003, 0.993545, 0.862859, 0.619299, 2.866666, 0.99317, 0.870024, 0.626189, 2.87451, 0.992831, 0.877168, 0.633109, 2.882352, 0.99244, 0.88433, 0.640099, 2.890196, 0.992089, 0.89147, 0.647116, 2.89804, 0.991688, 0.898627, 0.654202, 2.905882, 0.991332, 0.905763, 0.661309, 2.913726, 0.99093, 0.912915, 0.668481, 2.9215679999999997, 0.99057, 0.920049, 0.675675, 2.929412, 0.990175, 0.927196, 0.682926, 2.9372540000000003, 0.989815, 0.934329, 0.690198, 2.9450979999999998, 0.989434, 0.94147, 0.697519, 2.952942, 0.989077, 0.948604, 0.704863, 2.9607840000000003, 0.988717, 0.955742, 0.712242, 2.968628, 0.988367, 0.962878, 0.719649, 2.97647, 0.988033, 0.970012, 0.727077, 2.984314, 0.987691, 0.977154, 0.734536, 2.992156, 0.987387, 0.984288, 0.742002, 3.0, 0.987053, 0.991438, 0.749504]
+typeLUT.UseLogScale = 0
+typeLUT.ColorSpace = 'RGB'
+typeLUT.UseBelowRangeColor = 0
+typeLUT.BelowRangeColor = [0.0, 0.0, 0.0]
+typeLUT.UseAboveRangeColor = 0
+typeLUT.AboveRangeColor = [0.5, 0.5, 0.5]
+typeLUT.NanColor = [0.0, 1.0, 0.0]
+typeLUT.NanOpacity = 1.0
+typeLUT.Discretize = 1
+typeLUT.NumberOfTableValues = 256
+typeLUT.ScalarRangeInitialized = 1.0
+typeLUT.HSVWrap = 0
+typeLUT.VectorComponent = 0
+typeLUT.VectorMode = 'Magnitude'
+typeLUT.AllowDuplicateScalars = 1
+typeLUT.Annotations = []
+typeLUT.ActiveAnnotatedValues = []
+typeLUT.IndexedColors = []
+typeLUT.IndexedOpacities = []
+
+# get opacity transfer function/opacity map for 'Type'
+typePWF = GetOpacityTransferFunction('Type')
+typePWF.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.4276729533762702, 0.625, 0.5, 0.0, 1.7232704162597656, 0.6691176295280457, 0.5, 0.0, 2.3332666599993335, 0.8676470518112183, 0.5, 0.0, 2.471698045730591, 0.6911764740943909, 0.5, 0.0, 3.0, 0.8014705777168274, 0.5, 0.0]
+typePWF.AllowDuplicateScalars = 1
+typePWF.UseLogScale = 0
+typePWF.ScalarRangeInitialized = 1
+
+# trace defaults for the display properties.
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Representation = 'Surface'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.AmbientColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ColorArrayName = ['POINTS', 'Type']
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DiffuseColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.LookupTable = typeLUT
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.MapScalars = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.MultiComponentsMapping = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.InterpolateScalarsBeforeMapping = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Opacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PointSize = 2.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.LineWidth = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.RenderLinesAsTubes = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.RenderPointsAsSpheres = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Interpolation = 'Gouraud'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Specular = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SpecularColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SpecularPower = 100.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Luminosity = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Ambient = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Diffuse = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.EdgeColor = [0.0, 0.0, 0.5]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.BackfaceRepresentation = 'Follow Frontface'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.BackfaceAmbientColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.BackfaceDiffuseColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.BackfaceOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Position = [0.0, 0.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Scale = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Orientation = [0.0, 0.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Origin = [0.0, 0.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Pickable = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Texture = None
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Triangulate = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseShaderReplacements = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ShaderReplacements = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.NonlinearSubdivisionLevel = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseDataPartitions = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayUseScaleArray = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayScaleArray = 'Type'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayScaleFunction = 'PiecewiseFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayMaterial = 'None'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Orient = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OrientationMode = 'Direction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectOrientationVectors = 'None'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Scaling = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleMode = 'No Data Scaling Off'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleFactor = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectScaleArray = 'Type'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType = 'Arrow'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseGlyphTable = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphTableIndexArray = 'Type'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseCompositeGlyphTable = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseGlyphCullingAndLOD = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.LODValues = []
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ColorByLODIndex = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GaussianRadius = 0.05
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ShaderPreset = 'Sphere'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.CustomTriangleScale = 3
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.CustomShader = """ // This custom shader code define a gaussian blur
+ // Please take a look into vtkSMPointGaussianRepresentation.cxx
+ // for other custom shader examples
+ //VTK::Color::Impl
+ float dist2 = dot(offsetVCVSOutput.xy,offsetVCVSOutput.xy);
+ float gaussian = exp(-0.5*dist2);
+ opacity = opacity*gaussian;
+"""
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.Emissive = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleByArray = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SetScaleArray = ['POINTS', 'Type']
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleArrayComponent = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.UseScaleFunction = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleTransferFunction = 'PiecewiseFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityByArray = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityArray = ['POINTS', 'Type']
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityArrayComponent = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityTransferFunction = 'PiecewiseFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid = 'GridAxesRepresentation'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelColor = [0.0, 1.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelFontSize = 18
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelJustification = 'Left'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionCellLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelColor = [1.0, 1.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelFontSize = 18
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelJustification = 'Left'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectionPointLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes = 'PolarAxesRepresentation'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScalarOpacityFunction = typePWF
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScalarOpacityUnitDistance = 0.03988178566826777
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SelectMapper = 'Projected tetra'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SamplingDimensions = [128, 128, 128]
+
+# init the 'PiecewiseFunction' selected for 'OSPRayScaleFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayScaleFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.213836476688135, 0.625, 0.5, 0.0, 1.3616352081298828, 0.6691176295280457, 0.5, 0.0, 1.6666333299996667, 0.8676470518112183, 0.5, 0.0, 1.7358490228652954, 0.6911764740943909, 0.5, 0.0, 2.0, 0.8014705777168274, 0.5, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OSPRayScaleFunction.UseLogScale = 0
+
+# init the 'Arrow' selected for 'GlyphType'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.TipResolution = 6
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.TipRadius = 0.1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.TipLength = 0.35
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.ShaftResolution = 6
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.ShaftRadius = 0.03
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.GlyphType.Invert = 0
+
+# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleTransferFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.4276729533762702, 0.625, 0.5, 0.0, 1.7232704162597656, 0.6691176295280457, 0.5, 0.0, 2.3332666599993335, 0.8676470518112183, 0.5, 0.0, 2.471698045730591, 0.6911764740943909, 0.5, 0.0, 3.0, 0.8014705777168274, 0.5, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.ScaleTransferFunction.UseLogScale = 0
+
+# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityTransferFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.4276729533762702, 0.625, 0.5, 0.0, 1.7232704162597656, 0.6691176295280457, 0.5, 0.0, 2.3332666599993335, 0.8676470518112183, 0.5, 0.0, 2.471698045730591, 0.6911764740943909, 0.5, 0.0, 3.0, 0.8014705777168274, 0.5, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.OpacityTransferFunction.UseLogScale = 0
+
+# init the 'GridAxesRepresentation' selected for 'DataAxesGrid'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitle = 'X Axis'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitle = 'Y Axis'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitle = 'Z Axis'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XTitleOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YTitleOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZTitleOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.FacesToRender = 63
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.CullBackface = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.CullFrontface = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.GridColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ShowGrid = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ShowEdges = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ShowTicks = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.LabelUniqueEdgesOnly = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.AxesToLabel = 63
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XAxisNotation = 'Mixed'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XAxisPrecision = 2
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XAxisUseCustomLabels = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.XAxisLabels = []
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YAxisNotation = 'Mixed'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YAxisPrecision = 2
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YAxisUseCustomLabels = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.YAxisLabels = []
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZAxisNotation = 'Mixed'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZAxisPrecision = 2
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZAxisUseCustomLabels = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.ZAxisLabels = []
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.UseCustomBounds = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.DataAxesGrid.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
+
+# init the 'PolarAxesRepresentation' selected for 'PolarAxes'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Visibility = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Translation = [0.0, 0.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Scale = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Orientation = [0.0, 0.0, 0.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.EnableCustomBounds = [0, 0, 0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.EnableCustomRange = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.CustomRange = [0.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialAxesVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.DrawRadialGridlines = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarArcsVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.DrawPolarArcsGridlines = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.NumberOfRadialAxes = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.AutoSubdividePolarAxis = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.NumberOfPolarAxis = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.MinimumRadius = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.MinimumAngle = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.MaximumAngle = 90.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialAxesOriginToPolarAxis = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Ratio = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarArcsColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryPolarArcsColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitle = 'Radial Distance'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleLocation = 'Bottom'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarLabelVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarLabelFormat = '%-#6.3g'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarLabelExponentLocation = 'Labels'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialLabelVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialLabelFormat = '%-#3.1f'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialLabelLocation = 'Bottom'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.RadialUnitsVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ScreenSize = 10.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTitleFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisLabelFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTextFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextColor = [1.0, 1.0, 1.0]
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextOpacity = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextFontFamily = 'Arial'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextFontFile = ''
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextBold = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextItalic = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextShadow = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SecondaryRadialAxesTextFontSize = 12
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.EnableDistanceLOD = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.DistanceLODThreshold = 0.7
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.EnableViewAngleLOD = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ViewAngleLODThreshold = 0.7
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.SmallestVisiblePolarAngle = 0.5
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarTicksVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcTicksOriginToPolarAxis = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.TickLocation = 'Both'
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.AxisTickVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.AxisMinorTickVisibility = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcTickVisibility = 1
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcMinorTickVisibility = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.DeltaAngleMajor = 10.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.DeltaAngleMinor = 5.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisMajorTickSize = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTickRatioSize = 0.3
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisMajorTickThickness = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.PolarAxisTickRatioThickness = 0.5
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisMajorTickSize = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTickRatioSize = 0.3
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisMajorTickThickness = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.LastRadialAxisTickRatioThickness = 0.5
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcMajorTickSize = 0.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcTickRatioSize = 0.3
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcMajorTickThickness = 1.0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.ArcTickRatioThickness = 0.5
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.Use2DMode = 0
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.PolarAxes.UseLogAxis = 0
+
+# reset view to fit data
+renderView1.ResetCamera()
+
+#changing interaction mode based on data extents
+renderView1.InteractionMode = '2D'
+renderView1.CameraPosition = [0.0, 10010.0, 0.5]
+renderView1.CameraFocalPoint = [0.0, 10.0, 0.5]
+renderView1.CameraViewUp = [1.0, 0.0, 0.0]
+
+# show color bar/color legend
+phaseDiagram2DGamma_infinity55_4000spvtsDisplay.SetScalarBarVisibility(renderView1, True)
+
+# update the view to ensure updated data information
+renderView1.Update()
+
+# create a new 'Transform'
+transform1 = Transform(Input=phaseDiagram2DGamma_infinity55_4000spvts)
+transform1.Transform = 'Transform'
+transform1.TransformAllInputVectors = 1
+
+# init the 'Transform' selected for 'Transform'
+transform1.Transform.Translate = [0.0, 0.0, 0.0]
+transform1.Transform.Rotate = [0.0, 0.0, 0.0]
+transform1.Transform.Scale = [1.0, 1.0, 1.0]
+
+# Properties modified on transform1.Transform
+transform1.Transform.Scale = [0.1, 1.0, 1.0]
+
+# show data in view
+transform1Display = Show(transform1, renderView1)
+
+# trace defaults for the display properties.
+transform1Display.Representation = 'Surface'
+transform1Display.AmbientColor = [1.0, 1.0, 1.0]
+transform1Display.ColorArrayName = ['POINTS', 'Type']
+transform1Display.DiffuseColor = [1.0, 1.0, 1.0]
+transform1Display.LookupTable = typeLUT
+transform1Display.MapScalars = 1
+transform1Display.MultiComponentsMapping = 0
+transform1Display.InterpolateScalarsBeforeMapping = 1
+transform1Display.Opacity = 1.0
+transform1Display.PointSize = 2.0
+transform1Display.LineWidth = 1.0
+transform1Display.RenderLinesAsTubes = 0
+transform1Display.RenderPointsAsSpheres = 0
+transform1Display.Interpolation = 'Gouraud'
+transform1Display.Specular = 0.0
+transform1Display.SpecularColor = [1.0, 1.0, 1.0]
+transform1Display.SpecularPower = 100.0
+transform1Display.Luminosity = 0.0
+transform1Display.Ambient = 0.0
+transform1Display.Diffuse = 1.0
+transform1Display.EdgeColor = [0.0, 0.0, 0.5]
+transform1Display.BackfaceRepresentation = 'Follow Frontface'
+transform1Display.BackfaceAmbientColor = [1.0, 1.0, 1.0]
+transform1Display.BackfaceDiffuseColor = [1.0, 1.0, 1.0]
+transform1Display.BackfaceOpacity = 1.0
+transform1Display.Position = [0.0, 0.0, 0.0]
+transform1Display.Scale = [1.0, 1.0, 1.0]
+transform1Display.Orientation = [0.0, 0.0, 0.0]
+transform1Display.Origin = [0.0, 0.0, 0.0]
+transform1Display.Pickable = 1
+transform1Display.Texture = None
+transform1Display.Triangulate = 0
+transform1Display.UseShaderReplacements = 0
+transform1Display.ShaderReplacements = ''
+transform1Display.NonlinearSubdivisionLevel = 1
+transform1Display.UseDataPartitions = 0
+transform1Display.OSPRayUseScaleArray = 0
+transform1Display.OSPRayScaleArray = 'Type'
+transform1Display.OSPRayScaleFunction = 'PiecewiseFunction'
+transform1Display.OSPRayMaterial = 'None'
+transform1Display.Orient = 0
+transform1Display.OrientationMode = 'Direction'
+transform1Display.SelectOrientationVectors = 'None'
+transform1Display.Scaling = 0
+transform1Display.ScaleMode = 'No Data Scaling Off'
+transform1Display.ScaleFactor = 0.1
+transform1Display.SelectScaleArray = 'Type'
+transform1Display.GlyphType = 'Arrow'
+transform1Display.UseGlyphTable = 0
+transform1Display.GlyphTableIndexArray = 'Type'
+transform1Display.UseCompositeGlyphTable = 0
+transform1Display.UseGlyphCullingAndLOD = 0
+transform1Display.LODValues = []
+transform1Display.ColorByLODIndex = 0
+transform1Display.GaussianRadius = 0.005
+transform1Display.ShaderPreset = 'Sphere'
+transform1Display.CustomTriangleScale = 3
+transform1Display.CustomShader = """ // This custom shader code define a gaussian blur
+ // Please take a look into vtkSMPointGaussianRepresentation.cxx
+ // for other custom shader examples
+ //VTK::Color::Impl
+ float dist2 = dot(offsetVCVSOutput.xy,offsetVCVSOutput.xy);
+ float gaussian = exp(-0.5*dist2);
+ opacity = opacity*gaussian;
+"""
+transform1Display.Emissive = 0
+transform1Display.ScaleByArray = 0
+transform1Display.SetScaleArray = ['POINTS', 'Type']
+transform1Display.ScaleArrayComponent = ''
+transform1Display.UseScaleFunction = 1
+transform1Display.ScaleTransferFunction = 'PiecewiseFunction'
+transform1Display.OpacityByArray = 0
+transform1Display.OpacityArray = ['POINTS', 'Type']
+transform1Display.OpacityArrayComponent = ''
+transform1Display.OpacityTransferFunction = 'PiecewiseFunction'
+transform1Display.DataAxesGrid = 'GridAxesRepresentation'
+transform1Display.SelectionCellLabelBold = 0
+transform1Display.SelectionCellLabelColor = [0.0, 1.0, 0.0]
+transform1Display.SelectionCellLabelFontFamily = 'Arial'
+transform1Display.SelectionCellLabelFontFile = ''
+transform1Display.SelectionCellLabelFontSize = 18
+transform1Display.SelectionCellLabelItalic = 0
+transform1Display.SelectionCellLabelJustification = 'Left'
+transform1Display.SelectionCellLabelOpacity = 1.0
+transform1Display.SelectionCellLabelShadow = 0
+transform1Display.SelectionPointLabelBold = 0
+transform1Display.SelectionPointLabelColor = [1.0, 1.0, 0.0]
+transform1Display.SelectionPointLabelFontFamily = 'Arial'
+transform1Display.SelectionPointLabelFontFile = ''
+transform1Display.SelectionPointLabelFontSize = 18
+transform1Display.SelectionPointLabelItalic = 0
+transform1Display.SelectionPointLabelJustification = 'Left'
+transform1Display.SelectionPointLabelOpacity = 1.0
+transform1Display.SelectionPointLabelShadow = 0
+transform1Display.PolarAxes = 'PolarAxesRepresentation'
+transform1Display.ScalarOpacityFunction = typePWF
+transform1Display.ScalarOpacityUnitDistance = 0.005557396853320044
+transform1Display.SelectMapper = 'Projected tetra'
+transform1Display.SamplingDimensions = [128, 128, 128]
+
+# init the 'PiecewiseFunction' selected for 'OSPRayScaleFunction'
+transform1Display.OSPRayScaleFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.213836476688135, 0.625, 0.5, 0.0, 1.3616352081298828, 0.6691176295280457, 0.5, 0.0, 1.6666333299996667, 0.8676470518112183, 0.5, 0.0, 1.7358490228652954, 0.6911764740943909, 0.5, 0.0, 2.0, 0.8014705777168274, 0.5, 0.0]
+transform1Display.OSPRayScaleFunction.UseLogScale = 0
+
+# init the 'Arrow' selected for 'GlyphType'
+transform1Display.GlyphType.TipResolution = 6
+transform1Display.GlyphType.TipRadius = 0.1
+transform1Display.GlyphType.TipLength = 0.35
+transform1Display.GlyphType.ShaftResolution = 6
+transform1Display.GlyphType.ShaftRadius = 0.03
+transform1Display.GlyphType.Invert = 0
+
+# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
+transform1Display.ScaleTransferFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.4276729533762702, 0.625, 0.5, 0.0, 1.7232704162597656, 0.6691176295280457, 0.5, 0.0, 2.3332666599993335, 0.8676470518112183, 0.5, 0.0, 2.471698045730591, 0.6911764740943909, 0.5, 0.0, 3.0, 0.8014705777168274, 0.5, 0.0]
+transform1Display.ScaleTransferFunction.UseLogScale = 0
+
+# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
+transform1Display.OpacityTransferFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.4276729533762702, 0.625, 0.5, 0.0, 1.7232704162597656, 0.6691176295280457, 0.5, 0.0, 2.3332666599993335, 0.8676470518112183, 0.5, 0.0, 2.471698045730591, 0.6911764740943909, 0.5, 0.0, 3.0, 0.8014705777168274, 0.5, 0.0]
+transform1Display.OpacityTransferFunction.UseLogScale = 0
+
+# init the 'GridAxesRepresentation' selected for 'DataAxesGrid'
+transform1Display.DataAxesGrid.XTitle = 'X Axis'
+transform1Display.DataAxesGrid.YTitle = 'Y Axis'
+transform1Display.DataAxesGrid.ZTitle = 'Z Axis'
+transform1Display.DataAxesGrid.XTitleColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.XTitleFontFamily = 'Arial'
+transform1Display.DataAxesGrid.XTitleFontFile = ''
+transform1Display.DataAxesGrid.XTitleBold = 0
+transform1Display.DataAxesGrid.XTitleItalic = 0
+transform1Display.DataAxesGrid.XTitleFontSize = 12
+transform1Display.DataAxesGrid.XTitleShadow = 0
+transform1Display.DataAxesGrid.XTitleOpacity = 1.0
+transform1Display.DataAxesGrid.YTitleColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.YTitleFontFamily = 'Arial'
+transform1Display.DataAxesGrid.YTitleFontFile = ''
+transform1Display.DataAxesGrid.YTitleBold = 0
+transform1Display.DataAxesGrid.YTitleItalic = 0
+transform1Display.DataAxesGrid.YTitleFontSize = 12
+transform1Display.DataAxesGrid.YTitleShadow = 0
+transform1Display.DataAxesGrid.YTitleOpacity = 1.0
+transform1Display.DataAxesGrid.ZTitleColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.ZTitleFontFamily = 'Arial'
+transform1Display.DataAxesGrid.ZTitleFontFile = ''
+transform1Display.DataAxesGrid.ZTitleBold = 0
+transform1Display.DataAxesGrid.ZTitleItalic = 0
+transform1Display.DataAxesGrid.ZTitleFontSize = 12
+transform1Display.DataAxesGrid.ZTitleShadow = 0
+transform1Display.DataAxesGrid.ZTitleOpacity = 1.0
+transform1Display.DataAxesGrid.FacesToRender = 63
+transform1Display.DataAxesGrid.CullBackface = 0
+transform1Display.DataAxesGrid.CullFrontface = 1
+transform1Display.DataAxesGrid.GridColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.ShowGrid = 0
+transform1Display.DataAxesGrid.ShowEdges = 1
+transform1Display.DataAxesGrid.ShowTicks = 1
+transform1Display.DataAxesGrid.LabelUniqueEdgesOnly = 1
+transform1Display.DataAxesGrid.AxesToLabel = 63
+transform1Display.DataAxesGrid.XLabelColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.XLabelFontFamily = 'Arial'
+transform1Display.DataAxesGrid.XLabelFontFile = ''
+transform1Display.DataAxesGrid.XLabelBold = 0
+transform1Display.DataAxesGrid.XLabelItalic = 0
+transform1Display.DataAxesGrid.XLabelFontSize = 12
+transform1Display.DataAxesGrid.XLabelShadow = 0
+transform1Display.DataAxesGrid.XLabelOpacity = 1.0
+transform1Display.DataAxesGrid.YLabelColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.YLabelFontFamily = 'Arial'
+transform1Display.DataAxesGrid.YLabelFontFile = ''
+transform1Display.DataAxesGrid.YLabelBold = 0
+transform1Display.DataAxesGrid.YLabelItalic = 0
+transform1Display.DataAxesGrid.YLabelFontSize = 12
+transform1Display.DataAxesGrid.YLabelShadow = 0
+transform1Display.DataAxesGrid.YLabelOpacity = 1.0
+transform1Display.DataAxesGrid.ZLabelColor = [1.0, 1.0, 1.0]
+transform1Display.DataAxesGrid.ZLabelFontFamily = 'Arial'
+transform1Display.DataAxesGrid.ZLabelFontFile = ''
+transform1Display.DataAxesGrid.ZLabelBold = 0
+transform1Display.DataAxesGrid.ZLabelItalic = 0
+transform1Display.DataAxesGrid.ZLabelFontSize = 12
+transform1Display.DataAxesGrid.ZLabelShadow = 0
+transform1Display.DataAxesGrid.ZLabelOpacity = 1.0
+transform1Display.DataAxesGrid.XAxisNotation = 'Mixed'
+transform1Display.DataAxesGrid.XAxisPrecision = 2
+transform1Display.DataAxesGrid.XAxisUseCustomLabels = 0
+transform1Display.DataAxesGrid.XAxisLabels = []
+transform1Display.DataAxesGrid.YAxisNotation = 'Mixed'
+transform1Display.DataAxesGrid.YAxisPrecision = 2
+transform1Display.DataAxesGrid.YAxisUseCustomLabels = 0
+transform1Display.DataAxesGrid.YAxisLabels = []
+transform1Display.DataAxesGrid.ZAxisNotation = 'Mixed'
+transform1Display.DataAxesGrid.ZAxisPrecision = 2
+transform1Display.DataAxesGrid.ZAxisUseCustomLabels = 0
+transform1Display.DataAxesGrid.ZAxisLabels = []
+transform1Display.DataAxesGrid.UseCustomBounds = 0
+transform1Display.DataAxesGrid.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
+
+# init the 'PolarAxesRepresentation' selected for 'PolarAxes'
+transform1Display.PolarAxes.Visibility = 0
+transform1Display.PolarAxes.Translation = [0.0, 0.0, 0.0]
+transform1Display.PolarAxes.Scale = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.Orientation = [0.0, 0.0, 0.0]
+transform1Display.PolarAxes.EnableCustomBounds = [0, 0, 0]
+transform1Display.PolarAxes.CustomBounds = [0.0, 1.0, 0.0, 1.0, 0.0, 1.0]
+transform1Display.PolarAxes.EnableCustomRange = 0
+transform1Display.PolarAxes.CustomRange = [0.0, 1.0]
+transform1Display.PolarAxes.PolarAxisVisibility = 1
+transform1Display.PolarAxes.RadialAxesVisibility = 1
+transform1Display.PolarAxes.DrawRadialGridlines = 1
+transform1Display.PolarAxes.PolarArcsVisibility = 1
+transform1Display.PolarAxes.DrawPolarArcsGridlines = 1
+transform1Display.PolarAxes.NumberOfRadialAxes = 0
+transform1Display.PolarAxes.AutoSubdividePolarAxis = 1
+transform1Display.PolarAxes.NumberOfPolarAxis = 0
+transform1Display.PolarAxes.MinimumRadius = 0.0
+transform1Display.PolarAxes.MinimumAngle = 0.0
+transform1Display.PolarAxes.MaximumAngle = 90.0
+transform1Display.PolarAxes.RadialAxesOriginToPolarAxis = 1
+transform1Display.PolarAxes.Ratio = 1.0
+transform1Display.PolarAxes.PolarAxisColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.PolarArcsColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.LastRadialAxisColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.SecondaryPolarArcsColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.SecondaryRadialAxesColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.PolarAxisTitleVisibility = 1
+transform1Display.PolarAxes.PolarAxisTitle = 'Radial Distance'
+transform1Display.PolarAxes.PolarAxisTitleLocation = 'Bottom'
+transform1Display.PolarAxes.PolarLabelVisibility = 1
+transform1Display.PolarAxes.PolarLabelFormat = '%-#6.3g'
+transform1Display.PolarAxes.PolarLabelExponentLocation = 'Labels'
+transform1Display.PolarAxes.RadialLabelVisibility = 1
+transform1Display.PolarAxes.RadialLabelFormat = '%-#3.1f'
+transform1Display.PolarAxes.RadialLabelLocation = 'Bottom'
+transform1Display.PolarAxes.RadialUnitsVisibility = 1
+transform1Display.PolarAxes.ScreenSize = 10.0
+transform1Display.PolarAxes.PolarAxisTitleColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.PolarAxisTitleOpacity = 1.0
+transform1Display.PolarAxes.PolarAxisTitleFontFamily = 'Arial'
+transform1Display.PolarAxes.PolarAxisTitleFontFile = ''
+transform1Display.PolarAxes.PolarAxisTitleBold = 0
+transform1Display.PolarAxes.PolarAxisTitleItalic = 0
+transform1Display.PolarAxes.PolarAxisTitleShadow = 0
+transform1Display.PolarAxes.PolarAxisTitleFontSize = 12
+transform1Display.PolarAxes.PolarAxisLabelColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.PolarAxisLabelOpacity = 1.0
+transform1Display.PolarAxes.PolarAxisLabelFontFamily = 'Arial'
+transform1Display.PolarAxes.PolarAxisLabelFontFile = ''
+transform1Display.PolarAxes.PolarAxisLabelBold = 0
+transform1Display.PolarAxes.PolarAxisLabelItalic = 0
+transform1Display.PolarAxes.PolarAxisLabelShadow = 0
+transform1Display.PolarAxes.PolarAxisLabelFontSize = 12
+transform1Display.PolarAxes.LastRadialAxisTextColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.LastRadialAxisTextOpacity = 1.0
+transform1Display.PolarAxes.LastRadialAxisTextFontFamily = 'Arial'
+transform1Display.PolarAxes.LastRadialAxisTextFontFile = ''
+transform1Display.PolarAxes.LastRadialAxisTextBold = 0
+transform1Display.PolarAxes.LastRadialAxisTextItalic = 0
+transform1Display.PolarAxes.LastRadialAxisTextShadow = 0
+transform1Display.PolarAxes.LastRadialAxisTextFontSize = 12
+transform1Display.PolarAxes.SecondaryRadialAxesTextColor = [1.0, 1.0, 1.0]
+transform1Display.PolarAxes.SecondaryRadialAxesTextOpacity = 1.0
+transform1Display.PolarAxes.SecondaryRadialAxesTextFontFamily = 'Arial'
+transform1Display.PolarAxes.SecondaryRadialAxesTextFontFile = ''
+transform1Display.PolarAxes.SecondaryRadialAxesTextBold = 0
+transform1Display.PolarAxes.SecondaryRadialAxesTextItalic = 0
+transform1Display.PolarAxes.SecondaryRadialAxesTextShadow = 0
+transform1Display.PolarAxes.SecondaryRadialAxesTextFontSize = 12
+transform1Display.PolarAxes.EnableDistanceLOD = 1
+transform1Display.PolarAxes.DistanceLODThreshold = 0.7
+transform1Display.PolarAxes.EnableViewAngleLOD = 1
+transform1Display.PolarAxes.ViewAngleLODThreshold = 0.7
+transform1Display.PolarAxes.SmallestVisiblePolarAngle = 0.5
+transform1Display.PolarAxes.PolarTicksVisibility = 1
+transform1Display.PolarAxes.ArcTicksOriginToPolarAxis = 1
+transform1Display.PolarAxes.TickLocation = 'Both'
+transform1Display.PolarAxes.AxisTickVisibility = 1
+transform1Display.PolarAxes.AxisMinorTickVisibility = 0
+transform1Display.PolarAxes.ArcTickVisibility = 1
+transform1Display.PolarAxes.ArcMinorTickVisibility = 0
+transform1Display.PolarAxes.DeltaAngleMajor = 10.0
+transform1Display.PolarAxes.DeltaAngleMinor = 5.0
+transform1Display.PolarAxes.PolarAxisMajorTickSize = 0.0
+transform1Display.PolarAxes.PolarAxisTickRatioSize = 0.3
+transform1Display.PolarAxes.PolarAxisMajorTickThickness = 1.0
+transform1Display.PolarAxes.PolarAxisTickRatioThickness = 0.5
+transform1Display.PolarAxes.LastRadialAxisMajorTickSize = 0.0
+transform1Display.PolarAxes.LastRadialAxisTickRatioSize = 0.3
+transform1Display.PolarAxes.LastRadialAxisMajorTickThickness = 1.0
+transform1Display.PolarAxes.LastRadialAxisTickRatioThickness = 0.5
+transform1Display.PolarAxes.ArcMajorTickSize = 0.0
+transform1Display.PolarAxes.ArcTickRatioSize = 0.3
+transform1Display.PolarAxes.ArcMajorTickThickness = 1.0
+transform1Display.PolarAxes.ArcTickRatioThickness = 0.5
+transform1Display.PolarAxes.Use2DMode = 0
+transform1Display.PolarAxes.UseLogAxis = 0
+
+# hide data in view
+Hide(phaseDiagram2DGamma_infinity55_4000spvts, renderView1)
+
+# show color bar/color legend
+transform1Display.SetScalarBarVisibility(renderView1, True)
+
+# update the view to ensure updated data information
+renderView1.Update()
+
+# toggle 3D widget visibility (only when running from the GUI)
+Hide3DWidgets(proxy=transform1.Transform)
+
+# reset view to fit data
+renderView1.ResetCamera()
+
+# change representation type
+transform1Display.SetRepresentationType('Point Gaussian')
+
+# Properties modified on transform1Display
+transform1Display.ShaderPreset = 'Plain circle'
+
+# set scalar coloring
+ColorBy(transform1Display, ('POINTS', 'angles'))
+
+# Hide the scalar bar for this color map if no visible data is colored by it.
+HideScalarBarIfNotNeeded(typeLUT, renderView1)
+
+# rescale color and/or opacity maps used to include current data range
+transform1Display.RescaleTransferFunctionToDataRange(True, False)
+
+# show color bar/color legend
+transform1Display.SetScalarBarVisibility(renderView1, True)
+
+# get color transfer function/color map for 'angles'
+anglesLUT = GetColorTransferFunction('angles')
+anglesLUT.AutomaticRescaleRangeMode = "Grow and update on 'Apply'"
+anglesLUT.InterpretValuesAsCategories = 0
+anglesLUT.AnnotationsInitialized = 0
+anglesLUT.ShowCategoricalColorsinDataRangeOnly = 0
+anglesLUT.RescaleOnVisibilityChange = 0
+anglesLUT.EnableOpacityMapping = 0
+anglesLUT.RGBPoints = [0.0, 0.001462, 0.000466, 0.013866, 0.006160663193689554, 0.002258, 0.001295, 0.018331, 0.012319755591052443, 0.003279, 0.002305, 0.023708, 0.018480418784742, 0.004512, 0.00349, 0.029965, 0.024639511182104886, 0.00595, 0.004843, 0.03713, 0.03080017437579444, 0.007588, 0.006356, 0.044973, 0.03695926677315698, 0.009426, 0.008022, 0.052844, 0.043119929966846535, 0.011465, 0.009828, 0.06075, 0.04928059316053644, 0.013708, 0.011771, 0.068667, 0.05543968555789898, 0.016156, 0.01384, 0.076603, 0.061600348751588534, 0.018815, 0.016026, 0.084584, 0.06775944114895142, 0.021692, 0.01832, 0.09261, 0.07392010434264097, 0.024792, 0.020715, 0.100676, 0.08007919674000387, 0.028123, 0.023201, 0.108787, 0.08623985993369342, 0.031696, 0.025765, 0.116965, 0.09240052312738298, 0.03552, 0.028397, 0.125209, 0.09855961552474586, 0.039608, 0.03109, 0.133515, 0.10472027871843542, 0.04383, 0.03383, 0.141886, 0.1108793711157983, 0.048062, 0.036607, 0.150327, 0.11704003430948787, 0.05232, 0.039407, 0.158841, 0.1231991267068504, 0.056615, 0.04216, 0.167446, 0.12935978990053995, 0.060949, 0.044794, 0.176129, 0.13552045309422986, 0.06533, 0.047318, 0.184892, 0.1416795454915924, 0.069764, 0.049726, 0.193735, 0.14784020868528194, 0.074257, 0.052017, 0.20266, 0.15399930108264484, 0.078815, 0.054184, 0.211667, 0.1601599642763344, 0.083446, 0.056225, 0.220755, 0.16631905667369729, 0.088155, 0.058133, 0.229922, 0.17247971986738683, 0.092949, 0.059904, 0.239164, 0.17863881226474973, 0.097833, 0.061531, 0.248477, 0.18479947545843928, 0.102815, 0.06301, 0.257854, 0.19096013865212885, 0.107899, 0.064335, 0.267289, 0.19711923104949172, 0.113094, 0.065492, 0.276784, 0.20327989424318127, 0.118405, 0.066479, 0.286321, 0.20943898664054383, 0.123833, 0.067295, 0.295879, 0.21559964983423338, 0.12938, 0.067935, 0.305443, 0.22175874223159628, 0.135053, 0.068391, 0.315, 0.22791940542528583, 0.140858, 0.068654, 0.324538, 0.23408006861897537, 0.146785, 0.068738, 0.334011, 0.24023916101633827, 0.152839, 0.068637, 0.343404, 0.24639982421002782, 0.159018, 0.068354, 0.352688, 0.2525589166073907, 0.165308, 0.067911, 0.361816, 0.25871957980108024, 0.171713, 0.067305, 0.370771, 0.26487867219844313, 0.178212, 0.066576, 0.379497, 0.2710393353921327, 0.184801, 0.065732, 0.387973, 0.2771999985858223, 0.19146, 0.064818, 0.396152, 0.2833590909831848, 0.198177, 0.063862, 0.404009, 0.2895197541768747, 0.204935, 0.062907, 0.411514, 0.29567884657423726, 0.211718, 0.061992, 0.418647, 0.3018395097679268, 0.218512, 0.061158, 0.425392, 0.3079986021652897, 0.225302, 0.060445, 0.431742, 0.31415926535897926, 0.232077, 0.059889, 0.437695, 0.3203199285526688, 0.238826, 0.059517, 0.443256, 0.3264790209500317, 0.245543, 0.059352, 0.448436, 0.33263968414372125, 0.25222, 0.059415, 0.453248, 0.33879877654108415, 0.258857, 0.059706, 0.45771, 0.34495943973477367, 0.265447, 0.060237, 0.46184, 0.35111853213213656, 0.271994, 0.060994, 0.46566, 0.35727919532582614, 0.278493, 0.061978, 0.46919, 0.36343985851951566, 0.284951, 0.063168, 0.472451, 0.36959895091687855, 0.291366, 0.064553, 0.475462, 0.37575961411056813, 0.29774, 0.066117, 0.478243, 0.38191870650793064, 0.304081, 0.067835, 0.480812, 0.3880793697016202, 0.310382, 0.069702, 0.483186, 0.3942384620989831, 0.316654, 0.07169, 0.48538, 0.4003991252926727, 0.322899, 0.073782, 0.487408, 0.4065597884863622, 0.329114, 0.075972, 0.489287, 0.4127188808837251, 0.335308, 0.078236, 0.491024, 0.4188795440774147, 0.341482, 0.080564, 0.492631, 0.4250386364747776, 0.347636, 0.082946, 0.494121, 0.4311992996684671, 0.353773, 0.085373, 0.495501, 0.43735839206582966, 0.359898, 0.087831, 0.496778, 0.44351905525951957, 0.366012, 0.090314, 0.49796, 0.4496797184532091, 0.372116, 0.092816, 0.499053, 0.45583881085057165, 0.378211, 0.095332, 0.500067, 0.46199947404426156, 0.384299, 0.097855, 0.501002, 0.46815856644162407, 0.390384, 0.100379, 0.501864, 0.47431922963531364, 0.396467, 0.102902, 0.502658, 0.48047832203267654, 0.402548, 0.10542, 0.503386, 0.48663898522636606, 0.408629, 0.10793, 0.504052, 0.49279807762372896, 0.414709, 0.110431, 0.504662, 0.49895874081741853, 0.420791, 0.11292, 0.505215, 0.5051194040111081, 0.426877, 0.115395, 0.505714, 0.511278496408471, 0.432967, 0.117855, 0.50616, 0.5174391596021605, 0.439062, 0.120298, 0.506555, 0.5235982519995234, 0.445163, 0.122724, 0.506901, 0.5297589151932129, 0.451271, 0.125132, 0.507198, 0.5359180075905755, 0.457386, 0.127522, 0.507448, 0.5420786707842654, 0.463508, 0.129893, 0.507652, 0.5482393339779547, 0.46964, 0.132245, 0.507809, 0.5543984263753176, 0.47578, 0.134577, 0.507921, 0.560559089569007, 0.481929, 0.136891, 0.507989, 0.5667181819663699, 0.488088, 0.139186, 0.508011, 0.5728788451600595, 0.494258, 0.141462, 0.507988, 0.5790379375574224, 0.500438, 0.143719, 0.50792, 0.585198600751112, 0.506629, 0.145958, 0.507806, 0.5913592639448015, 0.512831, 0.148179, 0.507648, 0.5975183563421644, 0.519045, 0.150383, 0.507443, 0.6036790195358539, 0.52527, 0.152569, 0.507192, 0.6098381119332165, 0.531507, 0.154739, 0.506895, 0.6159987751269064, 0.537755, 0.156894, 0.506551, 0.6221578675242689, 0.544015, 0.159033, 0.506159, 0.6283185307179585, 0.550287, 0.161158, 0.505719, 0.6344791939116484, 0.556571, 0.163269, 0.50523, 0.6406382863090109, 0.562866, 0.165368, 0.504692, 0.6467989495027004, 0.569172, 0.167454, 0.504105, 0.6529580419000633, 0.57549, 0.16953, 0.503466, 0.6591187050937529, 0.581819, 0.171596, 0.502777, 0.6652777974911158, 0.588158, 0.173652, 0.502035, 0.6714384606848054, 0.594508, 0.175701, 0.501241, 0.677599123878495, 0.600868, 0.177743, 0.500394, 0.6837582162758579, 0.607238, 0.179779, 0.499492, 0.6899188794695473, 0.613617, 0.181811, 0.498536, 0.6960779718669102, 0.620005, 0.18384, 0.497524, 0.7022386350605995, 0.626401, 0.185867, 0.496456, 0.7083977274579624, 0.632805, 0.187893, 0.495332, 0.7145583906516519, 0.639216, 0.189921, 0.49415, 0.7207190538453414, 0.645633, 0.191952, 0.49291, 0.7268781462427043, 0.652056, 0.193986, 0.491611, 0.7330388094363939, 0.658483, 0.196027, 0.490253, 0.7391979018337568, 0.664915, 0.198075, 0.488836, 0.7453585650274464, 0.671349, 0.200133, 0.487358, 0.7515176574248092, 0.677786, 0.202203, 0.485819, 0.7576783206184988, 0.684224, 0.204286, 0.484219, 0.7638389838121884, 0.690661, 0.206384, 0.482558, 0.7699980762095512, 0.697098, 0.208501, 0.480835, 0.7761587394032408, 0.703532, 0.210638, 0.479049, 0.7823178318006033, 0.709962, 0.212797, 0.477201, 0.7884784949942932, 0.716387, 0.214982, 0.47529, 0.7946375873916561, 0.722805, 0.217194, 0.473316, 0.8007982505853454, 0.729216, 0.219437, 0.471279, 0.8069573429827082, 0.735616, 0.221713, 0.46918, 0.8131180061763977, 0.742004, 0.224025, 0.467018, 0.8192786693700873, 0.748378, 0.226377, 0.464794, 0.8254377617674502, 0.754737, 0.228772, 0.462509, 0.8315984249611398, 0.761077, 0.231214, 0.460162, 0.8377575173585023, 0.767398, 0.233705, 0.457755, 0.8439181805521923, 0.773695, 0.236249, 0.455289, 0.8500772729495552, 0.779968, 0.238851, 0.452765, 0.8562379361432443, 0.786212, 0.241514, 0.450184, 0.8623985993369342, 0.792427, 0.244242, 0.447543, 0.8685576917342971, 0.798608, 0.24704, 0.444848, 0.8747183549279863, 0.804752, 0.249911, 0.442102, 0.8808774473253492, 0.810855, 0.252861, 0.439305, 0.8870381105190391, 0.816914, 0.255895, 0.436461, 0.8931972029164016, 0.822926, 0.259016, 0.433573, 0.8993578661100912, 0.828886, 0.262229, 0.430644, 0.9055185293037807, 0.834791, 0.26554, 0.427671, 0.9116776217011436, 0.840636, 0.268953, 0.424666, 0.9178382848948332, 0.846416, 0.272473, 0.421631, 0.9239973772921961, 0.852126, 0.276106, 0.418573, 0.9301580404858857, 0.857763, 0.279857, 0.415496, 0.9363171328832481, 0.86332, 0.283729, 0.412403, 0.942477796076938, 0.868793, 0.287728, 0.409303, 0.9486384592706273, 0.874176, 0.291859, 0.406205, 0.9547975516679902, 0.879464, 0.296125, 0.403118, 0.9609582148616801, 0.884651, 0.30053, 0.400047, 0.9671173072590427, 0.889731, 0.305079, 0.397002, 0.9732779704527321, 0.8947, 0.309773, 0.393995, 0.979437062850095, 0.899552, 0.314616, 0.391037, 0.9855977260437846, 0.904281, 0.31961, 0.388137, 0.9917583892374742, 0.908884, 0.324755, 0.385308, 0.9979174816348371, 0.913354, 0.330052, 0.382563, 1.0040781448285265, 0.917689, 0.3355, 0.379915, 1.0102372372258892, 0.921884, 0.341098, 0.377376, 1.0163979004195791, 0.925937, 0.346844, 0.374959, 1.022556992816942, 0.929845, 0.352734, 0.372677, 1.0287176560106313, 0.933606, 0.358764, 0.370541, 1.034878319204321, 0.937221, 0.364929, 0.368567, 1.0410374116016838, 0.940687, 0.371224, 0.366762, 1.047198074795373, 0.944006, 0.377643, 0.365136, 1.053357167192736, 0.94718, 0.384178, 0.363701, 1.0595178303864259, 0.95021, 0.39082, 0.362468, 1.0656769227837886, 0.953099, 0.397563, 0.361438, 1.071837585977478, 0.955849, 0.4044, 0.360619, 1.0779982491711675, 0.958464, 0.411324, 0.360014, 1.0841573415685304, 0.960949, 0.418323, 0.35963, 1.09031800476222, 0.96331, 0.42539, 0.359469, 1.096477097159583, 0.965549, 0.432519, 0.359529, 1.1026377603532722, 0.967671, 0.439703, 0.35981, 1.108796852750635, 0.96968, 0.446936, 0.360311, 1.114957515944325, 0.971582, 0.45421, 0.36103, 1.121116608341688, 0.973381, 0.46152, 0.361965, 1.127277271535377, 0.975082, 0.468861, 0.363111, 1.1334379347290668, 0.97669, 0.476226, 0.364466, 1.1395970271264295, 0.97821, 0.483612, 0.366025, 1.145757690320119, 0.979645, 0.491014, 0.367783, 1.1519167827174819, 0.981, 0.498428, 0.369734, 1.1580774459111713, 0.982279, 0.505851, 0.371874, 1.164236538308534, 0.983485, 0.51328, 0.374198, 1.170397201502224, 0.984622, 0.520713, 0.376698, 1.1765578646959134, 0.985693, 0.528148, 0.379371, 1.182716957093276, 0.9867, 0.535582, 0.38221, 1.188877620286966, 0.987646, 0.543015, 0.38521, 1.1950367126843289, 0.988533, 0.550446, 0.388365, 1.201197375878018, 0.989363, 0.557873, 0.391671, 1.2073564682753808, 0.990138, 0.565296, 0.395122, 1.2135171314690707, 0.990871, 0.572706, 0.398714, 1.21967779466276, 0.991558, 0.580107, 0.402441, 1.2258368870601228, 0.992196, 0.587502, 0.406299, 1.2319975502538125, 0.992785, 0.594891, 0.410283, 1.2381566426511754, 0.993326, 0.602275, 0.41439, 1.2443173058448649, 0.993834, 0.609644, 0.418613, 1.2504763982422278, 0.994309, 0.616999, 0.42295, 1.2566370614359172, 0.994738, 0.62435, 0.427397, 1.262797724629607, 0.995122, 0.631696, 0.431951, 1.2689568170269698, 0.99548, 0.639027, 0.436607, 1.275117480220659, 0.99581, 0.646344, 0.441361, 1.2812765726180217, 0.996096, 0.653659, 0.446213, 1.2874372358117117, 0.996341, 0.660969, 0.45116, 1.2935963282090743, 0.99658, 0.668256, 0.456192, 1.2997569914027638, 0.996775, 0.675541, 0.461314, 1.3059176545964537, 0.996925, 0.682828, 0.466526, 1.3120767469938164, 0.997077, 0.690088, 0.471811, 1.3182374101875058, 0.997186, 0.697349, 0.477182, 1.3243965025848687, 0.997254, 0.704611, 0.482635, 1.3305571657785582, 0.997325, 0.711848, 0.488154, 1.3367162581759209, 0.997351, 0.719089, 0.493755, 1.3428769213696108, 0.997351, 0.726324, 0.499428, 1.3490375845633003, 0.997341, 0.733545, 0.505167, 1.355196676960663, 0.997285, 0.740772, 0.510983, 1.3613573401543528, 0.997228, 0.747981, 0.516859, 1.3675164325517157, 0.997138, 0.75519, 0.522806, 1.3736770957454048, 0.997019, 0.762398, 0.528821, 1.3798361881427677, 0.996898, 0.769591, 0.534892, 1.3859968513364576, 0.996727, 0.776795, 0.541039, 1.3921575145301468, 0.996571, 0.783977, 0.547233, 1.3983166069275097, 0.996369, 0.791167, 0.553499, 1.4044772701211992, 0.996162, 0.798348, 0.55982, 1.410636362518562, 0.995932, 0.805527, 0.566202, 1.4167970257122517, 0.99568, 0.812706, 0.572645, 1.4229561181096146, 0.995424, 0.819875, 0.57914, 1.4291167813033039, 0.995131, 0.827052, 0.585701, 1.4352758737006668, 0.994851, 0.834213, 0.592307, 1.4414365368943567, 0.994524, 0.841387, 0.598983, 1.4475972000880457, 0.994222, 0.84854, 0.605696, 1.4537562924854086, 0.993866, 0.855711, 0.612482, 1.4599169556790985, 0.993545, 0.862859, 0.619299, 1.4660760480764612, 0.99317, 0.870024, 0.626189, 1.4722367112701507, 0.992831, 0.877168, 0.633109, 1.4783958036675136, 0.99244, 0.88433, 0.640099, 1.4845564668612032, 0.992089, 0.89147, 0.647116, 1.4907171300548927, 0.991688, 0.898627, 0.654202, 1.4968762224522556, 0.991332, 0.905763, 0.661309, 1.503036885645945, 0.99093, 0.912915, 0.668481, 1.5091959780433077, 0.99057, 0.920049, 0.675675, 1.5153566412369976, 0.990175, 0.927196, 0.682926, 1.5215157336343605, 0.989815, 0.934329, 0.690198, 1.5276763968280496, 0.989434, 0.94147, 0.697519, 1.5338370600217395, 0.989077, 0.948604, 0.704863, 1.5399961524191024, 0.988717, 0.955742, 0.712242, 1.5461568156127916, 0.988367, 0.962878, 0.719649, 1.5523159080101545, 0.988033, 0.970012, 0.727077, 1.5584765712038442, 0.987691, 0.977154, 0.734536, 1.564635663601207, 0.987387, 0.984288, 0.742002, 1.5707963267948966, 0.987053, 0.991438, 0.749504]
+anglesLUT.UseLogScale = 0
+anglesLUT.ColorSpace = 'RGB'
+anglesLUT.UseBelowRangeColor = 0
+anglesLUT.BelowRangeColor = [0.0, 0.0, 0.0]
+anglesLUT.UseAboveRangeColor = 0
+anglesLUT.AboveRangeColor = [0.5, 0.5, 0.5]
+anglesLUT.NanColor = [0.0, 1.0, 0.0]
+anglesLUT.NanOpacity = 1.0
+anglesLUT.Discretize = 1
+anglesLUT.NumberOfTableValues = 256
+anglesLUT.ScalarRangeInitialized = 1.0
+anglesLUT.HSVWrap = 0
+anglesLUT.VectorComponent = 0
+anglesLUT.VectorMode = 'Magnitude'
+anglesLUT.AllowDuplicateScalars = 1
+anglesLUT.Annotations = []
+anglesLUT.ActiveAnnotatedValues = []
+anglesLUT.IndexedColors = []
+anglesLUT.IndexedOpacities = []
+
+# get opacity transfer function/opacity map for 'angles'
+anglesPWF = GetOpacityTransferFunction('angles')
+anglesPWF.Points = [0.0, 0.8602941036224365, 0.5, 0.0, 0.33589355211648514, 0.625, 0.5, 0.0, 0.5680552565701278, 0.6691176295280457, 0.5, 0.0, 1.0471451860825267, 0.8676470518112183, 0.5, 0.0, 1.15586894219242, 0.6911764740943909, 0.5, 0.0, 1.5707963267948966, 0.8014705777168274, 0.5, 0.0]
+anglesPWF.AllowDuplicateScalars = 1
+anglesPWF.UseLogScale = 0
+anglesPWF.ScalarRangeInitialized = 1
+
+# Apply a preset using its name. Note this may not work as expected when presets have duplicate names.
+anglesLUT.ApplyPreset('Cool to Warm', True)
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.XTitle = '$\\theta_\\rho$'
+renderView1.AxesGrid.ZTitle = '$\\theta'
+renderView1.AxesGrid.XTitleFontSize = 20
+renderView1.AxesGrid.ZTitleFontSize = 20
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.ShowGrid = 1
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.FacesToRender = 2
+renderView1.AxesGrid.CullFrontface = 0
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.Visibility = 1
+renderView1.AxesGrid.DataScale = [0.1, 1.0, 1.0]
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.ZTitle = '$\\theta$'
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.XTitleFontSize = 22
+renderView1.AxesGrid.ZTitleFontSize = 22
+renderView1.AxesGrid.ZAxisUseCustomLabels = 1
+renderView1.AxesGrid.ZAxisLabels = [0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.ShowEdges = 0
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.XAxisPrecision = 1
+renderView1.AxesGrid.XAxisLabels = [0.0]
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.ZAxisPrecision = 1
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.XAxisUseCustomLabels = 1
+renderView1.AxesGrid.XAxisLabels = [0.0, -5.0, -4.0, -3.0, -2.0, -1.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0]
+
+# Properties modified on renderView1.AxesGrid
+renderView1.AxesGrid.XLabelFontSize = 16
+renderView1.AxesGrid.ZLabelFontSize = 16
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['', '']
+anglesLUT.IndexedColors = [1.0, 1.0, 1.0]
+anglesLUT.IndexedOpacities = [1.0]
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '']
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '', '']
+anglesLUT.IndexedColors = [1.0, 1.0, 1.0, 1.0, 0.0, 0.0]
+anglesLUT.IndexedOpacities = [1.0, 1.0]
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '1.5707963267948966', '']
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '1.5707963267948966', '$\\pi/2$', '', '']
+anglesLUT.IndexedColors = [1.0, 1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0]
+anglesLUT.IndexedOpacities = [1.0, 1.0, 1.0]
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '1.5707963267948966', '$\\pi/2$', '0,785398163', '']
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '1.5707963267948966', '$\\pi/2$', '0,785398163', '$\\pi/4$']
+
+# Properties modified on anglesLUT
+anglesLUT.Annotations = ['0', '0', '1.5707963267948966', '$\\pi/2$', '0.785398163', '$\\pi/4$']
+
+# get color legend/bar for anglesLUT in view renderView1
+anglesLUTColorBar = GetScalarBar(anglesLUT, renderView1)
+anglesLUTColorBar.AutoOrient = 1
+anglesLUTColorBar.Orientation = 'Vertical'
+anglesLUTColorBar.WindowLocation = 'LowerRightCorner'
+anglesLUTColorBar.Position = [0.89, 0.02]
+anglesLUTColorBar.Title = 'angles'
+anglesLUTColorBar.ComponentTitle = ''
+anglesLUTColorBar.TitleJustification = 'Centered'
+anglesLUTColorBar.HorizontalTitle = 0
+anglesLUTColorBar.TitleColor = [1.0, 1.0, 1.0]
+anglesLUTColorBar.TitleOpacity = 1.0
+anglesLUTColorBar.TitleFontFamily = 'Arial'
+anglesLUTColorBar.TitleFontFile = ''
+anglesLUTColorBar.TitleBold = 0
+anglesLUTColorBar.TitleItalic = 0
+anglesLUTColorBar.TitleShadow = 0
+anglesLUTColorBar.TitleFontSize = 16
+anglesLUTColorBar.LabelColor = [1.0, 1.0, 1.0]
+anglesLUTColorBar.LabelOpacity = 1.0
+anglesLUTColorBar.LabelFontFamily = 'Arial'
+anglesLUTColorBar.LabelFontFile = ''
+anglesLUTColorBar.LabelBold = 0
+anglesLUTColorBar.LabelItalic = 0
+anglesLUTColorBar.LabelShadow = 0
+anglesLUTColorBar.LabelFontSize = 16
+anglesLUTColorBar.AutomaticLabelFormat = 1
+anglesLUTColorBar.LabelFormat = '%-#6.3g'
+anglesLUTColorBar.DrawTickMarks = 1
+anglesLUTColorBar.DrawTickLabels = 1
+anglesLUTColorBar.UseCustomLabels = 0
+anglesLUTColorBar.CustomLabels = []
+anglesLUTColorBar.AddRangeLabels = 1
+anglesLUTColorBar.RangeLabelFormat = '%-#6.1e'
+anglesLUTColorBar.DrawAnnotations = 1
+anglesLUTColorBar.AddRangeAnnotations = 0
+anglesLUTColorBar.AutomaticAnnotations = 0
+anglesLUTColorBar.DrawNanAnnotation = 0
+anglesLUTColorBar.NanAnnotation = 'NaN'
+anglesLUTColorBar.TextPosition = 'Ticks right/top, annotations left/bottom'
+anglesLUTColorBar.ReverseLegend = 0
+anglesLUTColorBar.ScalarBarThickness = 16
+anglesLUTColorBar.ScalarBarLength = 0.33
+
+# Properties modified on anglesLUTColorBar
+anglesLUTColorBar.DrawTickLabels = 0
+anglesLUTColorBar.AddRangeLabels = 0
+
+# Properties modified on anglesLUTColorBar
+anglesLUTColorBar.Title = 'angle $\\angle$'
+anglesLUTColorBar.HorizontalTitle = 1
+anglesLUTColorBar.TitleFontSize = 20
+
+# Properties modified on anglesLUTColorBar
+anglesLUTColorBar.LabelFontSize = 20
+
+# Properties modified on anglesLUTColorBar
+anglesLUTColorBar.DrawTickMarks = 0
+
+#### saving camera placements for all active views
+
+# current camera placement for renderView1
+renderView1.InteractionMode = '2D'
+renderView1.CameraPosition = [0.0, 7.295131721662212, 0.5]
+renderView1.CameraFocalPoint = [0.0, 10.0, 0.5]
+renderView1.CameraViewUp = [0.0, 0.0, 1.0]
+renderView1.CameraParallelScale = 0.7000714249274855
+
+#### uncomment the following to render all views
+# RenderAllViews()
+# alternatively, if you want to write images, you can use SaveScreenshot(...).
+# ----------------- EXPORT -----------------
+# save screenshot
+# SaveScreenshot('/home/klaus/Desktop/PhaseDiagramPlot.png', renderView1, ImageResolution=[1257, 934])
+# SaveScreenshot('/home/klaus/Desktop/PhaseDiagramPlot.png', renderView1, ImageResolution=[2514, 1868])
+# SaveScreenshot('/home/klaus/Desktop/PhaseDiagramPlot.png', renderView1, ImageResolution=[1063, 657])
+
+
+
+ SaveScreenshot('/home/klaus/Desktop/2DPhaseDiagram_1', renderView1, ImageResolution=[2126, 1314])
+# export view
+# ExportView('/home/klaus/Desktop/PhaseDiagramPlot.pdf', view=renderView1)
diff --git a/src/PhaseDiagram.py b/src/PhaseDiagram.py
index 2c7917d04e58fdf433d6799ee047edc395100d91..14011a66e8c5ef4cdb5ba115d0ff4986053549f1 100644
--- a/src/PhaseDiagram.py
+++ b/src/PhaseDiagram.py
@@ -110,7 +110,8 @@ print("Path: ", path)
# -------------------------- Input Parameters --------------------
-mu1 = 10.0 # TODO : here must be the same values as in the Parset for computeMuGamma
+# mu1 = 10.0 # TODO : here must be the same values as in the Parset for computeMuGamma
+mu1 = 1.0
rho1 = 1.0
alpha = 2.0
beta = 2.0
@@ -221,20 +222,21 @@ SamplePoints_3D = 300 # Number of sample points in each direction
# SamplePoints_3D = 100 # Number of sample points in each direction
# SamplePoints_3D = 200 # Number of sample points in each direction
# SamplePoints_3D = 400 # Number of sample points in each direction
-SamplePoints_2D = 300 # Number of sample points in each direction
-
+# SamplePoints_2D = 7500 # Number of sample points in each direction
+SamplePoints_2D = 40 # 4000 # Number of sample points in each direction
if make_3D_PhaseDiagram:
alphas_ = np.linspace(-20, 20, SamplePoints_3D)
+
# betas_ = np.linspace(0.01,40.01,SamplePoints_3D) # Full Range
- # betas_ = np.linspace(0.01,1.01,SamplePoints_3D) # weird part
- betas_ = np.linspace(1,40.01,SamplePoints_3D) #TEST !!!!! For Beta <1 weird tings happen...
+ betas_ = np.linspace(0.01,1.0,SamplePoints_3D) # weird part
+ # betas_ = np.linspace(1,40.01,SamplePoints_3D) #TEST !!!!! For Beta <1 weird tings happen...
thetas_ = np.linspace(0.01,0.99,SamplePoints_3D)
-
+ print('betas:', betas_)
# TEST :
# alphas_ = np.linspace(-40, 40, SamplePoints_3D)
@@ -262,15 +264,17 @@ if make_3D_PhaseDiagram:
print('Written to VTK-File:', VTKOutputName )
if make_2D_PhaseDiagram:
- alphas_ = np.linspace(-20, 20, SamplePoints_2D)
+ # alphas_ = np.linspace(-20, 20, SamplePoints_2D)
+ # alphas_ = np.linspace(0, 1, SamplePoints_2D)
thetas_ = np.linspace(0.01,0.99,SamplePoints_2D)
+ alphas_ = np.linspace(-5, 5, SamplePoints_2D)
# alphas_ = np.linspace(-5, 15, SamplePoints_2D)
# thetas_ = np.linspace(0.05,0.25,SamplePoints_2D)
# good range:
- alphas_ = np.linspace(9, 10, SamplePoints_2D)
- thetas_ = np.linspace(0.075,0.14,SamplePoints_2D)
+ # alphas_ = np.linspace(9, 10, SamplePoints_2D)
+ # thetas_ = np.linspace(0.075,0.14,SamplePoints_2D)
# range used:
# alphas_ = np.linspace(8, 10, SamplePoints_2D)
@@ -280,7 +284,14 @@ if make_2D_PhaseDiagram:
# thetas_ = np.linspace(0.05,0.2,SamplePoints_2D)
# betas_ = np.linspace(0.01,40.01,1)
#fix to one value:
- betas_ = 2.0;
+ # betas_ = 2.0;
+ betas_ = 10.0;
+
+ # TEST
+ # alphas_ = np.linspace(-8, 8, SamplePoints_2D)
+ # thetas_ = np.linspace(0.01,0.99,SamplePoints_2D)
+ # betas_ = 1.0; #TEST Problem: disvison by zero if alpha = 9, theta = 0.1 !
+ # betas_ = 0.9;
# betas_ = 0.5; #TEST!!!
alphas, betas, thetas = np.meshgrid(alphas_, betas_, thetas_, indexing='ij')
@@ -308,6 +319,7 @@ if make_2D_PhaseDiagram:
# --- Write to VTK
# VTKOutputName = + path + "./PhaseDiagram2DNEW"
+
GammaString = str(gamma)
VTKOutputName = "outputs/PhaseDiagram2D" + "Gamma_" + GammaString
gridToVTK(VTKOutputName , alphas, betas, thetas, pointData = {'Type': Types, 'angles': angles, 'curvature': curvature} )
diff --git a/src/PhaseDiagram_PlotScript.py b/src/PhaseDiagram_PlotScript.py
index 6884e2abd7f5f215e44fbed202839a0b107e4d7b..acfe626f71e5d41d7de97b01870afe2be902b296 100644
--- a/src/PhaseDiagram_PlotScript.py
+++ b/src/PhaseDiagram_PlotScript.py
@@ -17,7 +17,7 @@ case = 1
# case = 2
#
-curvature = 0
+curvature = 1
# create a new 'XML Structured Grid Reader'
@@ -27,8 +27,11 @@ curvature = 0
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinityBetaSmallerOne.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0300sp.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity_greaterOne.vts'])
-phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0_greaterOne.vts'])
+# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0_greaterOne.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0.vts'])
+# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity_300SpMu1.vts'])
+phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0_300spmu1.vts'])
+
phaseDiagram3DGammainfinityvts.PointArrayStatus = ['Type', 'angles', 'curvature']
@@ -351,15 +354,23 @@ anglesLUT.EnableOpacityMapping = 1
renderView1.AxesGrid.Visibility = 1
# Properties modified on renderView1.AxesGrid
-renderView1.AxesGrid.XTitle = '$\\theta_\\rho$'
-renderView1.AxesGrid.YTitle = '$\\theta_\\mu$'
-renderView1.AxesGrid.ZTitle = '$\\theta$'
-renderView1.AxesGrid.XTitleFontSize = 16
-renderView1.AxesGrid.YTitleFontSize = 16
-renderView1.AxesGrid.ZTitleFontSize = 16
-renderView1.AxesGrid.XTitleFontSize = 24
-renderView1.AxesGrid.YTitleFontSize = 24
-renderView1.AxesGrid.ZTitleFontSize = 24
+renderView1.AxesGrid.XTitle = ' $\\theta_\\rho$ '
+renderView1.AxesGrid.YTitle = ' $\\theta_\\mu$ '
+renderView1.AxesGrid.ZTitle = ' $\\theta$ '
+
+# renderView1.AxesGrid.XTitleFontSize = 45 #24
+# renderView1.AxesGrid.YTitleFontSize = 45
+# renderView1.AxesGrid.ZTitleFontSize = 45
+# renderView1.AxesGrid.XLabelFontSize = 30 #default:12
+# renderView1.AxesGrid.YLabelFontSize = 30
+# renderView1.AxesGrid.ZLabelFontSize = 30
+
+renderView1.AxesGrid.XTitleFontSize = 65 #24
+renderView1.AxesGrid.YTitleFontSize = 65
+renderView1.AxesGrid.ZTitleFontSize = 65
+renderView1.AxesGrid.XLabelFontSize = 45 #default:12
+renderView1.AxesGrid.YLabelFontSize = 45
+renderView1.AxesGrid.ZLabelFontSize = 45
# ---- Switch of/on OrientationAXES:
# renderView1.OrientationAxesVisibility = 1
@@ -390,15 +401,8 @@ renderView1.AxesGrid.ShowGrid = 1
# Properties modified on renderView1.AxesGrid
# --- Change FontSize of Values on Axis:
-renderView1.AxesGrid.XLabelFontSize = 12 #default:12
-renderView1.AxesGrid.YLabelFontSize = 12
-renderView1.AxesGrid.ZLabelFontSize = 12
-renderView1.AxesGrid.XLabelFontSize = 10 #default:12
-renderView1.AxesGrid.YLabelFontSize = 10
-renderView1.AxesGrid.ZLabelFontSize = 10
-renderView1.AxesGrid.XLabelFontSize = 16 #default:12
-renderView1.AxesGrid.YLabelFontSize = 16
-renderView1.AxesGrid.ZLabelFontSize = 16
+
+
# # --- Change Color of Values on Axis:
## -- Green:
# # renderView1.AxesGrid.XLabelColor = [0.5411764705882353, 0.8862745098039215, 0.20392156862745098]
@@ -425,13 +429,14 @@ renderView1.AxesGrid.ZAxisUseCustomLabels = 1
renderView1.AxesGrid.YAxisUseCustomLabels = 1
renderView1.AxesGrid.XAxisLabels = [-20, -15.0, -10.0, -5.0, 0.0, 5.0, 10.0, 15.0, 20.0]
# renderView1.AxesGrid.YAxisLabels = [0 , 10.0, 20.0, 30.0, 40.0]
-renderView1.AxesGrid.YAxisLabels = [1 ,5, 10.0,15, 20.0,25, 30.0,35, 40.0]
+# renderView1.AxesGrid.YAxisLabels = [1 ,5, 10.0,15, 20.0,25, 30.0,35, 40.0]
+renderView1.AxesGrid.YAxisLabels = [1 , 10, 20, 30, 40]
# renderView1.AxesGrid.YAxisLabels = [0 ,5.0, 10.0, 15.0, 20.0, 25.0, 30.0,35.0, 40.0]
renderView1.AxesGrid.ZAxisLabels = [ 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]
## ---- SET Center-Axis visible:
-renderView1.CenterAxesVisibility = 1
+renderView1.CenterAxesVisibility = 0
# ----------------- CAMERA POSITION -----------------
# current camera placement for renderView1
@@ -510,7 +515,7 @@ anglesLUTColorBar = GetScalarBar(anglesLUT, renderView1)
# Properties modified on anglesLUTColorBar
anglesLUTColorBar.WindowLocation = 'LowerRightCorner'
anglesLUTColorBar.HorizontalTitle = 1
-anglesLUTColorBar.Title = 'angle $\\angle$'
+anglesLUTColorBar.Title = 'angle $\\alpha$' #'angle $\\angle$'
anglesLUTColorBar.TextPosition = 'Ticks left/bottom, annotations right/top'
anglesLUTColorBar.TitleJustification = 'Left'
# anglesLUTColorBar.TitleJustification = 'Centered'
@@ -523,8 +528,8 @@ anglesLUTColorBar.AutomaticAnnotations = 0
anglesLUTColorBar.AddRangeLabels = 0
anglesLUTColorBar.DrawTickMarks = 0
anglesLUTColorBar.DrawTickLabels = 0
-anglesLUTColorBar.TitleFontSize = 20
-anglesLUTColorBar.LabelFontSize = 16
+anglesLUTColorBar.TitleFontSize = 24
+anglesLUTColorBar.LabelFontSize = 20
# anglesLUTColorBar.AddRangeAnnotations = 1
## -- Add Annotations :
@@ -801,6 +806,167 @@ if curvature == 1:
curvatureLUT.ApplyPreset('Cold and Hot', True)
+
+ #NEW:
+
+ # set active source
+ SetActiveSource(transform1)
+
+ # get color transfer function/color map for 'curvature'
+ curvatureLUT = GetColorTransferFunction('curvature')
+ curvatureLUT.RGBPoints = [-1187.69685, 0.0, 1.0, 1.0, -118.22955000000002, 0.0, 0.0, 1.0, 0.6001499999999851, 0.0, 0.0, 0.501960784314, 119.42984999999999, 1.0, 0.0, 0.0, 1188.89715, 1.0, 1.0, 0.0]
+ curvatureLUT.ColorSpace = 'RGB'
+ curvatureLUT.ScalarRangeInitialized = 1.0
+
+ # get opacity transfer function/opacity map for 'curvature'
+ curvaturePWF = GetOpacityTransferFunction('curvature')
+ curvaturePWF.Points = [-1187.69685, 0.8602941036224365, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+ curvaturePWF.ScalarRangeInitialized = 1
+
+ # Properties modified on curvatureLUT
+ curvatureLUT.EnableOpacityMapping = 1
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.6911764740943909, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.6691176295280457, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.6617646813392639, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.6323529481887817, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5073529481887817, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.47058823704719543, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.4779411852359772, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.4852941334247589, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -328.2367841697693, 0.6691176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -231.0803985595703, 0.5441176295280457, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -246.02752685546875, 0.5588235259056091, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -253.5010986328125, 0.5588235259056091, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -679.4943625218383, 0.625, 0.5, 0.0, -253.5010986328125, 0.5661764740943909, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # Properties modified on curvaturePWF
+ curvaturePWF.Points = [-1187.69685, 0.5, 0.5, 0.0, -694.4414672851562, 0.5147058963775635, 0.5, 0.0, -253.5010986328125, 0.5661764740943909, 0.5, 0.0, 396.61992227722794, 0.8676470518112183, 0.5, 0.0, 561.1175226475239, 0.6911764740943909, 0.5, 0.0, 1188.89715, 0.8014705777168274, 0.5, 0.0]
+
+ # create a new 'Contour'
+ contour1 = Contour(Input=transform1)
+ contour1.ContourBy = ['POINTS', 'Type']
+ contour1.Isosurfaces = [2.0]
+ contour1.PointMergeMethod = 'Uniform Binning'
+
+ # Properties modified on contour1
+ contour1.ContourBy = ['POINTS', 'curvature']
+ contour1.Isosurfaces = [0.0]
+
+ # show data in view
+ contour1Display = Show(contour1, renderView1)
+
+ # trace defaults for the display properties.
+ contour1Display.Representation = 'Surface'
+ contour1Display.ColorArrayName = ['POINTS', 'curvature']
+ contour1Display.LookupTable = curvatureLUT
+ contour1Display.OSPRayScaleArray = 'curvature'
+ contour1Display.OSPRayScaleFunction = 'PiecewiseFunction'
+ contour1Display.SelectOrientationVectors = 'None'
+ contour1Display.ScaleFactor = 0.09752500000000001
+ contour1Display.SelectScaleArray = 'curvature'
+ contour1Display.GlyphType = 'Arrow'
+ contour1Display.GlyphTableIndexArray = 'curvature'
+ contour1Display.GaussianRadius = 0.00487625
+ contour1Display.SetScaleArray = ['POINTS', 'curvature']
+ contour1Display.ScaleTransferFunction = 'PiecewiseFunction'
+ contour1Display.OpacityArray = ['POINTS', 'curvature']
+ contour1Display.OpacityTransferFunction = 'PiecewiseFunction'
+ contour1Display.DataAxesGrid = 'GridAxesRepresentation'
+ contour1Display.PolarAxes = 'PolarAxesRepresentation'
+
+ # init the 'PiecewiseFunction' selected for 'OSPRayScaleFunction'
+ contour1Display.OSPRayScaleFunction.Points = [1.0, 0.8602941036224365, 0.5, 0.0, 1.213836476688135, 0.625, 0.5, 0.0, 1.3616352081298828, 0.6691176295280457, 0.5, 0.0, 1.6666333299996667, 0.8676470518112183, 0.5, 0.0, 1.7358490228652954, 0.6911764740943909, 0.5, 0.0, 2.0, 0.8014705777168274, 0.5, 0.0]
+
+ # init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
+ contour1Display.ScaleTransferFunction.Points = [0.0, 0.8602941036224365, 0.5, 0.0, 2.5142493840581125e-39, 0.625, 0.5, 0.0, 4.252039284300157e-39, 0.6691176295280457, 0.5, 0.0, 7.838150278676329e-39, 0.8676470518112183, 0.5, 0.0, 8.651975477491057e-39, 0.6911764740943909, 0.5, 0.0, 1.1757813367477812e-38, 0.8014705777168274, 0.5, 0.0]
+
+ # init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
+ contour1Display.OpacityTransferFunction.Points = [0.0, 0.8602941036224365, 0.5, 0.0, 2.5142493840581125e-39, 0.625, 0.5, 0.0, 4.252039284300157e-39, 0.6691176295280457, 0.5, 0.0, 7.838150278676329e-39, 0.8676470518112183, 0.5, 0.0, 8.651975477491057e-39, 0.6911764740943909, 0.5, 0.0, 1.1757813367477812e-38, 0.8014705777168274, 0.5, 0.0]
+
+ # show color bar/color legend
+ contour1Display.SetScalarBarVisibility(renderView1, True)
+
+ # find source
+ xMLStructuredGridReader1 = FindSource('XMLStructuredGridReader1')
+
+ # update the view to ensure updated data information
+ renderView1.Update()
+
+ # set scalar coloring using an separate color/opacity maps
+ ColorBy(contour1Display, ('POINTS', 'curvature'), True)
+
+ # Hide the scalar bar for this color map if no visible data is colored by it.
+ HideScalarBarIfNotNeeded(curvatureLUT, renderView1)
+
+ # rescale color and/or opacity maps used to include current data range
+ contour1Display.RescaleTransferFunctionToDataRange(True, False)
+
+ # show color bar/color legend
+ contour1Display.SetScalarBarVisibility(renderView1, True)
+
+ # get separate color transfer function/color map for 'curvature'
+ separate_contour1Display_curvatureLUT = GetColorTransferFunction('curvature', contour1Display, separate=True)
+ separate_contour1Display_curvatureLUT.RGBPoints = [0.0, 0.001462, 0.000466, 0.013866, 4.611414402724776e-41, 0.002258, 0.001295, 0.018331, 9.2216530241129e-41, 0.003279, 0.002305, 0.023708, 1.3833067426837675e-40, 0.004512, 0.00349, 0.029965, 1.84433060482258e-40, 0.00595, 0.004843, 0.03713, 2.3054720450950578e-40, 0.007588, 0.006356, 0.044973, 2.766495907233844e-40, 0.009426, 0.008022, 0.052844, 3.2276373475063216e-40, 0.011465, 0.009828, 0.06075, 3.688778787778825e-40, 0.013708, 0.011771, 0.068667, 4.1498026499176114e-40, 0.016156, 0.01384, 0.076603, 4.6109440901900895e-40, 0.018815, 0.016026, 0.084584, 5.071967952328902e-40, 0.021692, 0.01832, 0.09261, 5.533109392601379e-40, 0.024792, 0.020715, 0.100676, 5.994133254740192e-40, 0.028123, 0.023201, 0.108787, 6.455274695012669e-40, 0.031696, 0.025765, 0.116965, 6.9164161352851465e-40, 0.03552, 0.028397, 0.125209, 7.37743999742396e-40, 0.039608, 0.03109, 0.133515, 7.838581437696437e-40, 0.04383, 0.03383, 0.141886, 8.299605299835249e-40, 0.048062, 0.036607, 0.150327, 8.760746740107726e-40, 0.05232, 0.039407, 0.158841, 9.221770602246514e-40, 0.056615, 0.04216, 0.167446, 9.682912042518991e-40, 0.060949, 0.044794, 0.176129, 1.0144053482791495e-39, 0.06533, 0.047318, 0.184892, 1.060507734493028e-39, 0.069764, 0.049726, 0.193735, 1.1066218785202758e-39, 0.074257, 0.052017, 0.20266, 1.152724264734157e-39, 0.078815, 0.054184, 0.211667, 1.1988384087614049e-39, 0.083446, 0.056225, 0.220755, 1.2449407949752861e-39, 0.088155, 0.058133, 0.229922, 1.2910549390025338e-39, 0.092949, 0.059904, 0.239164, 1.337157325216415e-39, 0.097833, 0.061531, 0.248477, 1.3832714692436628e-39, 0.102815, 0.06301, 0.257854, 1.4293856132709107e-39, 0.107899, 0.064335, 0.267289, 1.475487999484792e-39, 0.113094, 0.065492, 0.276784, 1.5216021435120396e-39, 0.118405, 0.066479, 0.286321, 1.5677045297259182e-39, 0.123833, 0.067295, 0.295879, 1.613818673753166e-39, 0.12938, 0.067935, 0.305443, 1.6599210599670472e-39, 0.135053, 0.068391, 0.315, 1.706035203994295e-39, 0.140858, 0.068654, 0.324538, 1.7521493480215426e-39, 0.146785, 0.068738, 0.334011, 1.7982517342354242e-39, 0.152839, 0.068637, 0.343404, 1.844365878262672e-39, 0.159018, 0.068354, 0.352688, 1.890468264476553e-39, 0.165308, 0.067911, 0.361816, 1.936582408503801e-39, 0.171713, 0.067305, 0.370771, 1.982684794717682e-39, 0.178212, 0.066576, 0.379497, 2.0287989387449298e-39, 0.184801, 0.065732, 0.387973, 2.0749130827721775e-39, 0.19146, 0.064818, 0.396152, 2.121015468986056e-39, 0.198177, 0.063862, 0.404009, 2.1671296130133065e-39, 0.204935, 0.062907, 0.411514, 2.213231999227185e-39, 0.211718, 0.061992, 0.418647, 2.2593461432544328e-39, 0.218512, 0.061158, 0.425392, 2.305448529468314e-39, 0.225302, 0.060445, 0.431742, 2.351562673495562e-39, 0.232077, 0.059889, 0.437695, 2.3976768175228098e-39, 0.238826, 0.059517, 0.443256, 2.443779203736691e-39, 0.245543, 0.059352, 0.448436, 2.4898933477639387e-39, 0.25222, 0.059415, 0.453248, 2.53599573397782e-39, 0.258857, 0.059706, 0.45771, 2.5821098780050677e-39, 0.265447, 0.060237, 0.46184, 2.628212264218949e-39, 0.271994, 0.060994, 0.46566, 2.6743264082461966e-39, 0.278493, 0.061978, 0.46919, 2.7204405522734444e-39, 0.284951, 0.063168, 0.472451, 2.7665429384873256e-39, 0.291366, 0.064553, 0.475462, 2.8126570825145736e-39, 0.29774, 0.066117, 0.478243, 2.8587594687284522e-39, 0.304081, 0.067835, 0.480812, 2.9048736127557e-39, 0.310382, 0.069702, 0.483186, 2.950975998969581e-39, 0.316654, 0.07169, 0.48538, 2.9970901429968286e-39, 0.322899, 0.073782, 0.487408, 3.0432042870240766e-39, 0.329114, 0.075972, 0.489287, 3.089306673237958e-39, 0.335308, 0.078236, 0.491024, 3.135420817265206e-39, 0.341482, 0.080564, 0.492631, 3.181523203479087e-39, 0.347636, 0.082946, 0.494121, 3.2276373475063345e-39, 0.353773, 0.085373, 0.495501, 3.273739733720213e-39, 0.359898, 0.087831, 0.496778, 3.319853877747464e-39, 0.366012, 0.090314, 0.49796, 3.365968021774711e-39, 0.372116, 0.092816, 0.499053, 3.41207040798859e-39, 0.378211, 0.095332, 0.500067, 3.4581845520158405e-39, 0.384299, 0.097855, 0.501002, 3.504286938229719e-39, 0.390384, 0.100379, 0.501864, 3.5504010822569665e-39, 0.396467, 0.102902, 0.502658, 3.5965034684708484e-39, 0.402548, 0.10542, 0.503386, 3.642617612498096e-39, 0.408629, 0.10793, 0.504052, 3.688719998711977e-39, 0.414709, 0.110431, 0.504662, 3.734834142739225e-39, 0.420791, 0.11292, 0.505215, 3.7809482867664724e-39, 0.426877, 0.115395, 0.505714, 3.8270506729803536e-39, 0.432967, 0.117855, 0.50616, 3.873164817007602e-39, 0.439062, 0.120298, 0.506555, 3.919267203221483e-39, 0.445163, 0.122724, 0.506901, 3.96538134724873e-39, 0.451271, 0.125132, 0.507198, 4.011483733462609e-39, 0.457386, 0.127522, 0.507448, 4.0575978774898596e-39, 0.463508, 0.129893, 0.507652, 4.103712021517105e-39, 0.46964, 0.132245, 0.507809, 4.149814407730986e-39, 0.47578, 0.134577, 0.507921, 4.1959285517582336e-39, 0.481929, 0.136891, 0.507989, 4.242030937972115e-39, 0.488088, 0.139186, 0.508011, 4.288145081999363e-39, 0.494258, 0.141462, 0.507988, 4.334247468213244e-39, 0.500438, 0.143719, 0.50792, 4.3803616122404915e-39, 0.506629, 0.145958, 0.507806, 4.4264757562677396e-39, 0.512831, 0.148179, 0.507648, 4.472578142481621e-39, 0.519045, 0.150383, 0.507443, 4.518692286508868e-39, 0.52527, 0.152569, 0.507192, 4.564794672722747e-39, 0.531507, 0.154739, 0.506895, 4.6109088167499975e-39, 0.537755, 0.156894, 0.506551, 4.657011202963876e-39, 0.544015, 0.159033, 0.506159, 4.703125346991124e-39, 0.550287, 0.161158, 0.505719, 4.749239491018374e-39, 0.556571, 0.163269, 0.50523, 4.795341877232253e-39, 0.562866, 0.165368, 0.504692, 4.841456021259501e-39, 0.569172, 0.167454, 0.504105, 4.887558407473382e-39, 0.57549, 0.16953, 0.503466, 4.9336725515006294e-39, 0.581819, 0.171596, 0.502777, 4.9797749377145106e-39, 0.588158, 0.173652, 0.502035, 5.025889081741759e-39, 0.594508, 0.175701, 0.501241, 5.072003225769006e-39, 0.600868, 0.177743, 0.500394, 5.118105611982888e-39, 0.607238, 0.179779, 0.499492, 5.1642197560101354e-39, 0.613617, 0.181811, 0.498536, 5.2103221422240166e-39, 0.620005, 0.18384, 0.497524, 5.256436286251262e-39, 0.626401, 0.185867, 0.496456, 5.302538672465143e-39, 0.632805, 0.187893, 0.495332, 5.3486528164923907e-39, 0.639216, 0.189921, 0.49415, 5.394766960519639e-39, 0.645633, 0.191952, 0.49291, 5.44086934673352e-39, 0.652056, 0.193986, 0.491611, 5.486983490760767e-39, 0.658483, 0.196027, 0.490253, 5.5330858769746485e-39, 0.664915, 0.198075, 0.488836, 5.5792000210018966e-39, 0.671349, 0.200133, 0.487358, 5.625302407215778e-39, 0.677786, 0.202203, 0.485819, 5.671416551243026e-39, 0.684224, 0.204286, 0.484219, 5.717530695270273e-39, 0.690661, 0.206384, 0.482558, 5.7636330814841545e-39, 0.697098, 0.208501, 0.480835, 5.8097472255114025e-39, 0.703532, 0.210638, 0.479049, 5.855849611725281e-39, 0.709962, 0.212797, 0.477201, 5.901963755752531e-39, 0.716387, 0.214982, 0.47529, 5.948066141966412e-39, 0.722805, 0.217194, 0.473316, 5.994180285993657e-39, 0.729216, 0.219437, 0.471279, 6.040282672207538e-39, 0.735616, 0.221713, 0.46918, 6.086396816234787e-39, 0.742004, 0.224025, 0.467018, 6.132510960262034e-39, 0.748378, 0.226377, 0.464794, 6.178613346475916e-39, 0.754737, 0.228772, 0.462509, 6.224727490503163e-39, 0.761077, 0.231214, 0.460162, 6.270829876717042e-39, 0.767398, 0.233705, 0.457755, 6.316944020744293e-39, 0.773695, 0.236249, 0.455289, 6.363046406958174e-39, 0.779968, 0.238851, 0.452765, 6.409160550985419e-39, 0.786212, 0.241514, 0.450184, 6.455274695012669e-39, 0.792427, 0.244242, 0.447543, 6.50137708122655e-39, 0.798608, 0.24704, 0.444848, 6.547491225253795e-39, 0.804752, 0.249911, 0.442102, 6.593593611467676e-39, 0.810855, 0.252861, 0.439305, 6.639707755494928e-39, 0.816914, 0.255895, 0.436461, 6.685810141708806e-39, 0.822926, 0.259016, 0.433573, 6.731924285736054e-39, 0.828886, 0.262229, 0.430644, 6.778038429763301e-39, 0.834791, 0.26554, 0.427671, 6.824140815977182e-39, 0.840636, 0.268953, 0.424666, 6.870254960004431e-39, 0.846416, 0.272473, 0.421631, 6.916357346218312e-39, 0.852126, 0.276106, 0.418573, 6.96247149024556e-39, 0.857763, 0.279857, 0.415496, 7.008573876459438e-39, 0.86332, 0.283729, 0.412403, 7.054688020486688e-39, 0.868793, 0.287728, 0.409303, 7.100802164513933e-39, 0.874176, 0.291859, 0.406205, 7.146904550727814e-39, 0.879464, 0.296125, 0.403118, 7.193018694755065e-39, 0.884651, 0.30053, 0.400047, 7.239121080968944e-39, 0.889731, 0.305079, 0.397002, 7.285235224996191e-39, 0.8947, 0.309773, 0.393995, 7.331337611210073e-39, 0.899552, 0.314616, 0.391037, 7.37745175523732e-39, 0.904281, 0.31961, 0.388137, 7.423565899264569e-39, 0.908884, 0.324755, 0.385308, 7.46966828547845e-39, 0.913354, 0.330052, 0.382563, 7.515782429505697e-39, 0.917689, 0.3355, 0.379915, 7.561884815719576e-39, 0.921884, 0.341098, 0.377376, 7.607998959746826e-39, 0.925937, 0.346844, 0.374959, 7.654101345960707e-39, 0.929845, 0.352734, 0.372677, 7.700215489987953e-39, 0.933606, 0.358764, 0.370541, 7.746329634015203e-39, 0.937221, 0.364929, 0.368567, 7.792432020229085e-39, 0.940687, 0.371224, 0.366762, 7.83854616425633e-39, 0.944006, 0.377643, 0.365136, 7.88464855047021e-39, 0.94718, 0.384178, 0.363701, 7.93076269449746e-39, 0.95021, 0.39082, 0.362468, 7.976865080711339e-39, 0.953099, 0.397563, 0.361438, 8.022979224738588e-39, 0.955849, 0.4044, 0.360619, 8.069093368765835e-39, 0.958464, 0.411324, 0.360014, 8.115195754979717e-39, 0.960949, 0.418323, 0.35963, 8.161309899006964e-39, 0.96331, 0.42539, 0.359469, 8.207412285220845e-39, 0.965549, 0.432519, 0.359529, 8.253526429248091e-39, 0.967671, 0.439703, 0.35981, 8.299628815461972e-39, 0.96968, 0.446936, 0.360311, 8.345742959489222e-39, 0.971582, 0.45421, 0.36103, 8.391845345703104e-39, 0.973381, 0.46152, 0.361965, 8.437959489730349e-39, 0.975082, 0.468861, 0.363111, 8.484073633757599e-39, 0.97669, 0.476226, 0.364466, 8.530176019971477e-39, 0.97821, 0.483612, 0.366025, 8.576290163998726e-39, 0.979645, 0.491014, 0.367783, 8.622392550212607e-39, 0.981, 0.498428, 0.369734, 8.668506694239854e-39, 0.982279, 0.505851, 0.371874, 8.714609080453733e-39, 0.983485, 0.51328, 0.374198, 8.760723224480983e-39, 0.984622, 0.520713, 0.376698, 8.806837368508232e-39, 0.985693, 0.528148, 0.379371, 8.85293975472211e-39, 0.9867, 0.535582, 0.38221, 8.89905389874936e-39, 0.987646, 0.543015, 0.38521, 8.945156284963242e-39, 0.988533, 0.550446, 0.388365, 8.991270428990486e-39, 0.989363, 0.557873, 0.391671, 9.037372815204368e-39, 0.990138, 0.565296, 0.395122, 9.083486959231618e-39, 0.990871, 0.572706, 0.398714, 9.129601103258864e-39, 0.991558, 0.580107, 0.402441, 9.175703489472745e-39, 0.992196, 0.587502, 0.406299, 9.221817633499992e-39, 0.992785, 0.594891, 0.410283, 9.267920019713874e-39, 0.993326, 0.602275, 0.41439, 9.314034163741121e-39, 0.993834, 0.609644, 0.418613, 9.360136549955002e-39, 0.994309, 0.616999, 0.42295, 9.406250693982251e-39, 0.994738, 0.62435, 0.427397, 9.452364838009498e-39, 0.995122, 0.631696, 0.431951, 9.49846722422338e-39, 0.99548, 0.639027, 0.436607, 9.544581368250624e-39, 0.99581, 0.646344, 0.441361, 9.590683754464506e-39, 0.996096, 0.653659, 0.446213, 9.636797898491756e-39, 0.996341, 0.660969, 0.45116, 9.682900284705634e-39, 0.99658, 0.668256, 0.456192, 9.729014428732883e-39, 0.996775, 0.675541, 0.461314, 9.775128572760133e-39, 0.996925, 0.682828, 0.466526, 9.821230958974011e-39, 0.997077, 0.690088, 0.471811, 9.867345103001259e-39, 0.997186, 0.697349, 0.477182, 9.91344748921514e-39, 0.997254, 0.704611, 0.482635, 9.959561633242389e-39, 0.997325, 0.711848, 0.488154, 1.0005664019456267e-38, 0.997351, 0.719089, 0.493755, 1.0051778163483517e-38, 0.997351, 0.726324, 0.499428, 1.0097892307510765e-38, 0.997341, 0.733545, 0.505167, 1.0143994693724643e-38, 0.997285, 0.740772, 0.510983, 1.0190108837751893e-38, 0.997228, 0.747981, 0.516859, 1.0236211223965776e-38, 0.997138, 0.75519, 0.522806, 1.0282325367993021e-38, 0.997019, 0.762398, 0.528821, 1.0328427754206902e-38, 0.996898, 0.769591, 0.534892, 1.0374541898234152e-38, 0.996727, 0.776795, 0.541039, 1.0420656042261397e-38, 0.996571, 0.783977, 0.547233, 1.0466758428475278e-38, 0.996369, 0.791167, 0.553499, 1.0512872572502527e-38, 0.996162, 0.798348, 0.55982, 1.0558974958716408e-38, 0.995932, 0.805527, 0.566202, 1.0605089102743655e-38, 0.99568, 0.812706, 0.572645, 1.0651191488957537e-38, 0.995424, 0.819875, 0.57914, 1.0697305632984781e-38, 0.995131, 0.827052, 0.585701, 1.0743408019198663e-38, 0.994851, 0.834213, 0.592307, 1.0789522163225914e-38, 0.994524, 0.841387, 0.598983, 1.0835636307253159e-38, 0.994222, 0.84854, 0.605696, 1.088173869346704e-38, 0.993866, 0.855711, 0.612482, 1.092785283749429e-38, 0.993545, 0.862859, 0.619299, 1.0973955223708168e-38, 0.99317, 0.870024, 0.626189, 1.1020069367735416e-38, 0.992831, 0.877168, 0.633109, 1.1066171753949297e-38, 0.99244, 0.88433, 0.640099, 1.1112285897976546e-38, 0.992089, 0.89147, 0.647116, 1.1158400042003793e-38, 0.991688, 0.898627, 0.654202, 1.1204502428217674e-38, 0.991332, 0.905763, 0.661309, 1.1250616572244922e-38, 0.99093, 0.912915, 0.668481, 1.12967189584588e-38, 0.99057, 0.920049, 0.675675, 1.1342833102486052e-38, 0.990175, 0.927196, 0.682926, 1.1388935488699933e-38, 0.989815, 0.934329, 0.690198, 1.1435049632727178e-38, 0.989434, 0.94147, 0.697519, 1.1481163776754428e-38, 0.989077, 0.948604, 0.704863, 1.1527266162968309e-38, 0.988717, 0.955742, 0.712242, 1.1573380306995554e-38, 0.988367, 0.962878, 0.719649, 1.1619482693209435e-38, 0.988033, 0.970012, 0.727077, 1.1665596837236684e-38, 0.987691, 0.977154, 0.734536, 1.1711699223450565e-38, 0.987387, 0.984288, 0.742002, 1.1757813367477812e-38, 0.987053, 0.991438, 0.749504]
+ separate_contour1Display_curvatureLUT.ColorSpace = 'RGB'
+ separate_contour1Display_curvatureLUT.NanColor = [0.0, 1.0, 0.0]
+ separate_contour1Display_curvatureLUT.ScalarRangeInitialized = 1.0
+
+ # get separate opacity transfer function/opacity map for 'curvature'
+ separate_contour1Display_curvaturePWF = GetOpacityTransferFunction('curvature', contour1Display, separate=True)
+ separate_contour1Display_curvaturePWF.Points = [0.0, 0.8602941036224365, 0.5, 0.0, 2.5142493840581125e-39, 0.625, 0.5, 0.0, 4.252039284300157e-39, 0.6691176295280457, 0.5, 0.0, 7.838150278676329e-39, 0.8676470518112183, 0.5, 0.0, 8.651975477491057e-39, 0.6911764740943909, 0.5, 0.0, 1.1757813367477812e-38, 0.8014705777168274, 0.5, 0.0]
+ separate_contour1Display_curvaturePWF.ScalarRangeInitialized = 1
+
+ # Apply a preset using its name. Note this may not work as expected when presets have duplicate names.
+ separate_contour1Display_curvatureLUT.ApplyPreset('X Ray', True)
+
+ # hide color bar/color legend
+ contour1Display.SetScalarBarVisibility(renderView1, False)
+
+ # hide data in view
+ Hide(line1, renderView1)
+
+
+
+ # get color legend/bar for curvatureLUT in view renderView1
+ curvatureLUTColorBar = GetScalarBar(curvatureLUT, renderView1)
+ curvatureLUTColorBar.Title = 'curvature'
+ curvatureLUTColorBar.ComponentTitle = ''
+
+ # Properties modified on curvatureLUTColorBar
+ curvatureLUTColorBar.Title = '$ \\kappa $'
+ curvatureLUTColorBar.HorizontalTitle = 1
+ curvatureLUTColorBar.TitleFontSize = 24
+ curvatureLUTColorBar.LabelFontSize = 20
+
+ transform1Display.Opacity = 0.01 #overall Opacity
+
+
# update the view to ensure updated data information
renderView1.Update()
diff --git a/src/PhaseDiagram_PlotScript_BetaSmallerOne.py b/src/PhaseDiagram_PlotScript_BetaSmallerOne.py
index ff25144f6ab38aef408d19a0da40511428def9d4..01f93eb18728f634274114ab6d037be04f31716b 100644
--- a/src/PhaseDiagram_PlotScript_BetaSmallerOne.py
+++ b/src/PhaseDiagram_PlotScript_BetaSmallerOne.py
@@ -3,6 +3,7 @@
#### import the simple module from the paraview
from paraview.simple import *
+from paraview.servermanager import *
#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()
@@ -14,18 +15,18 @@ paraview.simple._DisableFirstRenderCameraReset()
case = 1
case = 2
#
-# contour = 1
+contour = 0
# create a new 'XML Structured Grid Reader'
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity.vts'])
-phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity_100SP.vts'])
+# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity_100SP.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity300sp.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinityBetaSmallerOne.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0BetaSmallerOne.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0300sp.vts'])
# phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGamma0.vts'])
-
+phaseDiagram3DGammainfinityvts = XMLStructuredGridReader(FileName=['/home/klaus/Desktop/DUNE/dune-microstructure/outputs/PhaseDiagram3DGammainfinity_300sp_smallerOne_mu1.vts'])
phaseDiagram3DGammainfinityvts.PointArrayStatus = ['Type', 'angles', 'curvature']
@@ -347,9 +348,9 @@ anglesLUT.EnableOpacityMapping = 1
renderView1.AxesGrid.Visibility = 1
# Properties modified on renderView1.AxesGrid
-renderView1.AxesGrid.XTitle = '$\\theta_\\rho$'
-renderView1.AxesGrid.YTitle = '$\\theta_\\mu$'
-renderView1.AxesGrid.ZTitle = '$\\theta$'
+renderView1.AxesGrid.XTitle = ' $\\theta_\\rho$ '
+renderView1.AxesGrid.YTitle = ' $\\theta_\\mu$ '
+renderView1.AxesGrid.ZTitle = ' $\\theta$ '
renderView1.AxesGrid.XTitleFontSize = 16
renderView1.AxesGrid.YTitleFontSize = 16
renderView1.AxesGrid.ZTitleFontSize = 16
@@ -386,15 +387,25 @@ renderView1.AxesGrid.ShowGrid = 1
# Properties modified on renderView1.AxesGrid
# --- Change FontSize of Values on Axis:
-renderView1.AxesGrid.XLabelFontSize = 12 #default:12
-renderView1.AxesGrid.YLabelFontSize = 12
-renderView1.AxesGrid.ZLabelFontSize = 12
-renderView1.AxesGrid.XLabelFontSize = 10 #default:12
-renderView1.AxesGrid.YLabelFontSize = 10
-renderView1.AxesGrid.ZLabelFontSize = 10
-renderView1.AxesGrid.XLabelFontSize = 16 #default:12
-renderView1.AxesGrid.YLabelFontSize = 16
-renderView1.AxesGrid.ZLabelFontSize = 16
+# renderView1.AxesGrid.XLabelFontSize = 12 #default:12
+# renderView1.AxesGrid.YLabelFontSize = 12
+# renderView1.AxesGrid.ZLabelFontSize = 12
+# renderView1.AxesGrid.XLabelFontSize = 10 #default:12
+# renderView1.AxesGrid.YLabelFontSize = 10
+# renderView1.AxesGrid.ZLabelFontSize = 10
+# renderView1.AxesGrid.XLabelFontSize = 16 #default:12
+# renderView1.AxesGrid.YLabelFontSize = 16
+# renderView1.AxesGrid.ZLabelFontSize = 16
+
+
+
+
+renderView1.AxesGrid.XTitleFontSize = 65 #24
+renderView1.AxesGrid.YTitleFontSize = 65
+renderView1.AxesGrid.ZTitleFontSize = 65
+renderView1.AxesGrid.XLabelFontSize = 45 #default:12
+renderView1.AxesGrid.YLabelFontSize = 45
+renderView1.AxesGrid.ZLabelFontSize = 45
# # --- Change Color of Values on Axis:
## -- Green:
# # renderView1.AxesGrid.XLabelColor = [0.5411764705882353, 0.8862745098039215, 0.20392156862745098]
@@ -423,12 +434,13 @@ renderView1.AxesGrid.YAxisUseCustomLabels = 1
renderView1.AxesGrid.XAxisLabels = [-20, -15.0, -10.0, -5.0, 0.0, 5.0, 10.0, 15.0, 20.0]
# renderView1.AxesGrid.YAxisLabels = [0 , 10.0, 20.0, 30.0, 40.0]
# renderView1.AxesGrid.YAxisLabels = [0 ,5.0, 10.0, 15.0, 20.0, 25.0, 30.0,35.0, 40.0]
-renderView1.AxesGrid.YAxisLabels = [0 , 0.1, 0.2, 0.3 ,0.4, 0.5 ,0.6,0.7,0.8,0.9,1.0]
+# renderView1.AxesGrid.YAxisLabels = [0 , 0.1, 0.2, 0.3 ,0.4, 0.5 ,0.6,0.7,0.8,0.9,1.0]
+renderView1.AxesGrid.YAxisLabels = [0 , 0.2, 0.4 ,0.6,0.8,1.0]
renderView1.AxesGrid.ZAxisLabels = [ 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]
## ---- SET Center-Axis visible:
-renderView1.CenterAxesVisibility = 1
+renderView1.CenterAxesVisibility = 0
# ----------------- CAMERA POSITION -----------------
# current camera placement for renderView1
@@ -507,7 +519,8 @@ anglesLUTColorBar = GetScalarBar(anglesLUT, renderView1)
# Properties modified on anglesLUTColorBar
anglesLUTColorBar.WindowLocation = 'LowerRightCorner'
anglesLUTColorBar.HorizontalTitle = 1
-anglesLUTColorBar.Title = 'angle $\\angle$'
+anglesLUTColorBar.Title = 'angle $\\alpha$'
+# anglesLUTColorBar.Title = 'angle $\\angle$'
anglesLUTColorBar.TextPosition = 'Ticks left/bottom, annotations right/top'
anglesLUTColorBar.TitleJustification = 'Left'
# anglesLUTColorBar.TitleJustification = 'Centered'
@@ -520,8 +533,10 @@ anglesLUTColorBar.AutomaticAnnotations = 0
anglesLUTColorBar.AddRangeLabels = 0
anglesLUTColorBar.DrawTickMarks = 0
anglesLUTColorBar.DrawTickLabels = 0
-anglesLUTColorBar.TitleFontSize = 20
-anglesLUTColorBar.LabelFontSize = 16
+# anglesLUTColorBar.TitleFontSize = 20
+# anglesLUTColorBar.LabelFontSize = 16
+anglesLUTColorBar.TitleFontSize = 24
+anglesLUTColorBar.LabelFontSize = 20
# anglesLUTColorBar.AddRangeAnnotations = 1
## -- Add Annotations :
diff --git a/src/Plot-Angle-Gamma.py b/src/Plot-Angle-Gamma.py
index ea2007763fa5a3fcc86948d2f7297deae61444d9..643667099b66e20397637b20a7d673e57a790302 100644
--- a/src/Plot-Angle-Gamma.py
+++ b/src/Plot-Angle-Gamma.py
@@ -95,7 +95,8 @@ print('---- Input parameters: -----')
alpha = 10.0
# beta = 40.0
# theta = 0.125
-mu1 = 10.0
+# mu1 = 10.0
+mu1 = 1.0
rho1 = 1.0
# alpha = 2.0
beta = 2.0 #5.0
@@ -293,7 +294,8 @@ ax.plot(Gamma_Values, angles, 'royalblue', zorder=3, )
# ax.set_xlabel(r"$q_3(\gamma)$")
ax.set_xlabel(r"$\gamma$")
-ax.set_ylabel(r"angle $\angle$")
+# ax.set_ylabel(r"angle $\angle$")
+ax.set_ylabel(r"angle $\alpha$")
# plt.xlabel("$q_3$")
# plt.xlabel("$\gamma$")
diff --git a/src/Plot-Angle-q3.py b/src/Plot-Angle-q3.py
index 78e58067cf2c63624d6a29e52f205cea402f82b9..4fc1bd7e493b9bd6e30f727c06d4acc6fe5cda85 100644
--- a/src/Plot-Angle-q3.py
+++ b/src/Plot-Angle-q3.py
@@ -95,7 +95,8 @@ print('---- Input parameters: -----')
alpha = 10.0
# beta = 40.0
# theta = 0.125
-mu1 = 10.0
+# mu1 = 10.0
+mu1 = 1.0
rho1 = 1.0
# alpha = 2.0
beta = 2.0 #5.0
@@ -176,7 +177,7 @@ print('----------------------------')
gamma_min = 0.01
gamma_max = 3.0
-Gamma_Values = np.linspace(gamma_min, gamma_max, num=200) # TODO variable Input Parameters...alpha,beta...
+Gamma_Values = np.linspace(gamma_min, gamma_max, num=400) # TODO variable Input Parameters...alpha,beta...
print('(Input) Gamma_Values:', Gamma_Values)
# mu_gamma = []
@@ -262,14 +263,14 @@ fig = plt.figure()
ax = plt.axes((0.15,0.21 ,0.8,0.75))
ax.tick_params(axis='x',which='major', direction='out',pad=5)
ax.tick_params(axis='y',which='major', length=3, width=1, direction='out',pad=3)
-ax.xaxis.set_major_locator(MultipleLocator(0.02))
-ax.xaxis.set_minor_locator(MultipleLocator(0.01))
-ax.yaxis.set_major_locator(plt.MultipleLocator(np.pi / 8))
+# ax.xaxis.set_major_locator(MultipleLocator(0.1))
+# ax.xaxis.set_minor_locator(MultipleLocator(0.1))
+# ax.yaxis.set_major_locator(plt.MultipleLocator(np.pi / 8))
ax.yaxis.set_minor_locator(plt.MultipleLocator(np.pi / 16))
ax.yaxis.set_major_formatter(plt.FuncFormatter(format_func))
ax.grid(True,which='major',axis='both',alpha=0.3)
-ax.set(xlim=(1.760, 1.880), ylim=(-0.1, np.pi/4.0))
-
+# ax.set(xlim=(1.760, 1.880), ylim=(-0.1, np.pi/4.0))
+ax.set(ylim=(-0.1, np.pi/4.0))
# # plt.title(r'angle$-\mu_\gamma(\gamma)$-Plot')
# # plt.title(r'angle$-\gamma$-Plot')
@@ -292,12 +293,15 @@ ax.set(xlim=(1.760, 1.880), ylim=(-0.1, np.pi/4.0))
# ax.plot(muGammas, angles)
# ax.scatter(muGammas, angles)
-ax.plot(muGammas, angles, 'royalblue', zorder=3, )
+
+# print('muGammas:', muGammas)
+# print('angles:', angles)
+ax.plot(muGammas, angles, 'royalblue', zorder=3 )
# ax.scatter(Gamma_Values, angles)
ax.set_xlabel(r"$q_3(\gamma)$")
# ax.set_xlabel(r"$\gamma$")
-ax.set_ylabel(r"angle $\angle$")
-
+# ax.set_ylabel(r"angle $\angle$")
+ax.set_ylabel(r"angle $\alpha$")
# plt.xlabel("$q_3$")
# plt.xlabel("$\gamma$")
# plt.ylabel("angle")
diff --git a/src/Plot-Curvature-Gamma.py b/src/Plot-Curvature-Gamma.py
index 12184027f449ec451e934463b09c14c2077a1a38..b3302fec5be9eb386172ade14cbee095ca5adf24 100644
--- a/src/Plot-Curvature-Gamma.py
+++ b/src/Plot-Curvature-Gamma.py
@@ -95,7 +95,8 @@ print('---- Input parameters: -----')
alpha = 10.0
# beta = 40.0
# theta = 0.125
-mu1 = 10.0
+# mu1 = 10.0
+mu1 = 1.0
rho1 = 1.0
# alpha = 2.0
beta = 2.0 #5.0
diff --git a/src/Plot-Curvature-q3.py b/src/Plot-Curvature-q3.py
index e1088e04caef229b488f01e282631306efa84e27..95f8f5a1e6b0d5f8c4b32b9813ca6c06e1c3f2e9 100644
--- a/src/Plot-Curvature-q3.py
+++ b/src/Plot-Curvature-q3.py
@@ -95,7 +95,7 @@ print('---- Input parameters: -----')
alpha = 10.0
# beta = 40.0
# theta = 0.125
-mu1 = 10.0
+mu1 = 1.0
rho1 = 1.0
# alpha = 2.0
beta = 2.0 #5.0
@@ -267,14 +267,14 @@ ax = plt.axes((0.15,0.18 ,0.8,0.75))
ax.tick_params(axis='x',which='major', direction='out',pad=1)
# ax.tick_params(axis='y',which='major', length=3, width=1, direction='out',pad=3)
ax.tick_params(axis='y',which='major', length=3, width=1, direction='out',pad=1) # changed pad = distance to title to 1 here!
-ax.xaxis.set_major_locator(MultipleLocator(0.02))
-ax.xaxis.set_minor_locator(MultipleLocator(0.01))
+# ax.xaxis.set_major_locator(MultipleLocator(0.02))
+# ax.xaxis.set_minor_locator(MultipleLocator(0.01))
# ax.yaxis.set_major_locator(plt.MultipleLocator(np.pi / 8))
# ax.yaxis.set_minor_locator(plt.MultipleLocator(np.pi / 16))
# ax.yaxis.set_major_formatter(plt.FuncFormatter(format_func))
ax.grid(True,which='major',axis='both',alpha=0.3)
# ax.set(xlim=(1.760, 1.880), ylim=(-0.1, np.pi/4.0))
-ax.set(xlim=(1.760, 1.880))
+# ax.set(xlim=(1.760, 1.880))
# # plt.title(r'angle$-\mu_\gamma(\gamma)$-Plot')
diff --git a/src/Plot_Prestrain_Lemma1.4.py b/src/Plot_Prestrain_Lemma1.4.py
index adef5d8e9a2800f34f8cddba3601947c9f9b74f0..839cbe0b3f82ad72831fe4fae53bcef2b9da86c8 100644
--- a/src/Plot_Prestrain_Lemma1.4.py
+++ b/src/Plot_Prestrain_Lemma1.4.py
@@ -60,7 +60,7 @@ print("Path: ", path)
print('---- Input parameters: -----')
-mu1 = 10.0
+mu1 = 1.0
# mu1 = 10.0
# lambda1 = 10.0
rho1 = 1.0
diff --git a/src/Plot_aStar.py b/src/Plot_aStar.py
new file mode 100644
index 0000000000000000000000000000000000000000..5aed6de966a32b464614fd136ad4a2aa19e6f2aa
--- /dev/null
+++ b/src/Plot_aStar.py
@@ -0,0 +1,769 @@
+import numpy as np
+import matplotlib.pyplot as plt
+import sympy as sym
+import math
+import os
+import subprocess
+import fileinput
+import re
+import matlab.engine
+from HelperFunctions import *
+from ClassifyMin import *
+
+import matplotlib.ticker as tickers
+import matplotlib as mpl
+from matplotlib.ticker import MultipleLocator,FormatStrFormatter,MaxNLocator
+import pandas as pd
+
+# import tikzplotlib
+# # from pylab import *
+# from tikzplotlib import save as tikz_save
+
+
+# Needed ?
+mpl.use('pdf')
+
+# from subprocess import Popen, PIPE
+#import sys
+
+###################### makePlot.py #########################
+# Generalized Plot-Script giving the option to define
+# quantity of interest and the parameter it depends on
+# to create a plot
+#
+# Input: Define y & x for "x-y plot" as Strings
+# - Run the 'Cell-Problem' for the different Parameter-Points
+# (alternatively run 'Compute_MuGamma' if quantity of interest
+# is q3=muGamma for a significant Speedup)
+
+###########################################################
+
+
+
+# figsize argument takes inputs in inches
+# and we have the width of our document in pts.
+# To set the figure size we construct a function
+# to convert from pts to inches and to determine
+# an aesthetic figure height using the golden ratio:
+# def set_size(width, fraction=1):
+# """Set figure dimensions to avoid scaling in LaTeX.
+#
+# Parameters
+# ----------
+# width: float
+# Document textwidth or columnwidth in pts
+# fraction: float, optional
+# Fraction of the width which you wish the figure to occupy
+#
+# Returns
+# -------
+# fig_dim: tuple
+# Dimensions of figure in inches
+# """
+# # Width of figure (in pts)
+# fig_width_pt = width * fraction
+#
+# # Convert from pt to inches
+# inches_per_pt = 1 / 72.27
+#
+# # Golden ratio to set aesthetic figure height
+# # https://disq.us/p/2940ij3
+# golden_ratio = (5**.5 - 1) / 2
+#
+# # Figure width in inches
+# fig_width_in = fig_width_pt * inches_per_pt
+# # Figure height in inches
+# fig_height_in = fig_width_in * golden_ratio
+#
+# fig_dim = (fig_width_in, fig_height_in)
+#
+# return fig_dim
+#
+
+
+
+def format_func(value, tick_number):
+ # # find number of multiples of pi/2
+ # N = int(np.round(2 * value / np.pi))
+ # if N == 0:
+ # return "0"
+ # elif N == 1:
+ # return r"$\pi/2$"
+ # elif N == 2:
+ # return r"$\pi$"
+ # elif N % 2 > 0:
+ # return r"${0}\pi/2$".format(N)
+ # else:
+ # return r"${0}\pi$".format(N // 2)
+ # find number of multiples of pi/2
+ N = int(np.round(4 * value / np.pi))
+ if N == 0:
+ return "0"
+ elif N == 1:
+ return r"$\pi/4$"
+ elif N == 2:
+ return r"$\pi/2$"
+ elif N % 2 > 0:
+ return r"${0}\pi/2$".format(N)
+ else:
+ return r"${0}\pi$".format(N // 2)
+
+
+
+
+
+def find_nearest(array, value):
+ array = np.asarray(array)
+ idx = (np.abs(array - value)).argmin()
+ return array[idx]
+
+
+def find_nearestIdx(array, value):
+ array = np.asarray(array)
+ idx = (np.abs(array - value)).argmin()
+ return idx
+
+
+
+# TODO
+# - Fallunterscheidung (Speedup) falls gesuchter value mu_gamma = q3
+# - Also Add option to plot Minimization Output
+
+
+# ----- Setup Paths -----
+# InputFile = "/inputs/cellsolver.parset"
+# OutputFile = "/outputs/output.txt"
+
+InputFile = "/inputs/computeMuGamma.parset"
+OutputFile = "/outputs/outputMuGamma.txt"
+
+# path = os.getcwd()
+# InputFilePath = os.getcwd()+InputFile
+# OutputFilePath = os.getcwd()+OutputFile
+# --------- Run from src folder:
+path_parent = os.path.dirname(os.getcwd())
+os.chdir(path_parent)
+path = os.getcwd()
+print(path)
+InputFilePath = os.getcwd()+InputFile
+OutputFilePath = os.getcwd()+OutputFile
+print("InputFilepath: ", InputFilePath)
+print("OutputFilepath: ", OutputFilePath)
+print("Path: ", path)
+
+#---------------------------------------------------------------
+
+print('---- Input parameters: -----')
+mu1 = 10.0
+# lambda1 = 10.0
+rho1 = 1.0
+alpha = 5.0
+beta = 10.0
+theta = 1.0/4.0
+
+lambda1 = 0.0
+gamma = 1.0/4.0
+
+gamma = 'infinity' #Elliptic Setting
+# gamma = '0' #Hyperbolic Setting
+# gamma = 0.5
+
+
+print('mu1: ', mu1)
+print('rho1: ', rho1)
+print('alpha: ', alpha)
+print('beta: ', beta)
+print('theta: ', theta)
+print('gamma:', gamma)
+print('----------------------------')
+
+
+# TODO? : Ask User for Input ...
+# function = input("Enter value you want to plot (y-value):\n")
+# print(f'You entered {function}')
+# parameter = input("Enter Parameter this value depends on (x-value) :\n")
+# print(f'You entered {parameter}')
+
+# Add Option to change NumberOfElements used for computation of Cell-Problem
+
+
+# --- Define Quantity of interest:
+# Options: 'q1', 'q2', 'q3', 'q12' ,'q21', 'q31', 'q13' , 'q23', 'q32' , 'b1', 'b2' ,'b3'
+# TODO: EXTRA (MInimization Output) 'Minimizer (norm?)' 'angle', 'type', 'curvature'
+# yName = 'q12'
+# # yName = 'b1'
+# yName = 'q3'
+# yName = 'angle'
+# yName = 'curvature'
+yName = 'MinVec'
+
+# --- Define Parameter this function/quantity depends on:
+# Options: mu1 ,lambda1, rho1 , alpha, beta, theta, gamma
+# xName = 'theta'
+# xName = 'gamma'
+# xName = 'lambda1'
+xName = 'theta'
+
+
+# --- define Interval of x-va1ues:
+# xmin = 0.15
+xmin = 0.01
+xmax = 0.41
+
+# xmin = 0.18 #Achtung bei manchen werten von theta ist integration in ComputeMuGama/Cell_problem schlecht!
+# xmax = 0.41 # Materialfunktion muss von Gitter aufgelöst werden
+ # müssen vielfache von (1/2^i) sein wobei i integer
+
+
+# xmin = 0.18 #Achtung bei manchen werten von theta ist integration in ComputeMuGama/Cell_problem schlecht!
+# xmax = 0.23
+
+
+
+
+# xmin = 0.01
+# xmax = 3.0
+numPoints = 70
+# numPoints = 50
+X_Values = np.linspace(xmin, xmax, num=numPoints)
+print(X_Values)
+
+
+Y_Values = []
+
+
+
+
+
+
+
+
+for theta in X_Values:
+
+ print('Situation of Lemma1.4')
+ q12 = 0.0
+ q1 = (1.0/6.0)*harmonicMean(mu1, beta, theta)
+ q2 = (1.0/6.0)*arithmeticMean(mu1, beta, theta)
+ b1 = prestrain_b1(rho1, beta, alpha,theta)
+ b2 = prestrain_b2(rho1, beta, alpha,theta)
+ b3 = 0.0
+ # if gamma == '0':
+ # q3 = q2
+ # if gamma == 'infinity':
+ # q3 = q1
+ q3 = GetMuGamma(beta,theta,gamma,mu1,rho1,InputFilePath ,OutputFilePath)
+
+
+ if yName == 'q1': # TODO: Better use dictionary?...
+ print('q1 used')
+ Y_Values.append(q1)
+ elif yName =='q2':
+ print('q2 used')
+ Y_Values.append(q2)
+ elif yName =='q3':
+ print('q3 used')
+ Y_Values.append(q3)
+ elif yName =='q12':
+ print('q12 used')
+ Y_Values.append(q12)
+ elif yName =='b1':
+ print('b1 used')
+ Y_Values.append(b1)
+ elif yName =='b2':
+ print('b2 used')
+ Y_Values.append(b2)
+ elif yName =='b3':
+ print('b3 used')
+ Y_Values.append(b3)
+ elif yName == 'angle' or yName =='type' or yName =='curvature' or yName =='MinVec':
+ G, angle, Type, curvature = classifyMin_ana(alpha,beta,theta, q3, mu1, rho1)
+ if yName =='angle':
+ print('angle used')
+ Y_Values.append(angle)
+ if yName =='type':
+ print('angle used')
+ Y_Values.append(type)
+ if yName =='curvature':
+ print('angle used')
+ Y_Values.append(curvature)
+ if yName =='MinVec':
+ print('MinVec used')
+ Y_Values.append(G)
+
+
+print("(Output) Values of " + yName + ": ", Y_Values)
+
+
+# idx = find_nearestIdx(Y_Values, 0)
+# print(' Idx of value closest to 0', idx)
+# ValueClose = Y_Values[idx]
+# print('GammaValue(Idx) with mu_gamma closest to q_3^*', ValueClose)
+#
+#
+#
+# # Find Indices where the difference between the next one is larger than epsilon...
+# jump_idx = []
+# jump_xValues = []
+# jump_yValues = []
+# tmp = X_Values[0]
+# for idx, x in enumerate(X_Values):
+# print(idx, x)
+# if idx > 0:
+# if abs(Y_Values[idx]-Y_Values[idx-1]) > 1:
+# print('jump candidate')
+# jump_idx.append(idx)
+# jump_xValues.append(x)
+# jump_yValues.append(Y_Values[idx])
+#
+#
+#
+
+
+#
+#
+# print("Jump Indices", jump_idx)
+# print("Jump X-values:", jump_xValues)
+# print("Jump Y-values:", jump_yValues)
+#
+# y_plotValues = [Y_Values[0]]
+# x_plotValues = [X_Values[0]]
+# # y_plotValues.extend(jump_yValues)
+# for i in jump_idx:
+# y_plotValues.extend([Y_Values[i-1], Y_Values[i]])
+# x_plotValues.extend([X_Values[i-1], X_Values[i]])
+#
+#
+# y_plotValues.append(Y_Values[-1])
+# # x_plotValues = [X_Values[0]]
+# # x_plotValues.extend(jump_xValues)
+# x_plotValues.append(X_Values[-1])
+#
+#
+# print("y_plotValues:", y_plotValues)
+# print("x_plotValues:", x_plotValues)
+# Y_Values[np.diff(y) >= 0.5] = np.nan
+
+
+#get values bigger than jump position
+# gamma = infty
+# x_rest = X_Values[X_Values>x_plotValues[1]]
+# Y_Values = np.array(Y_Values) #convert the np array
+# y_rest = Y_Values[X_Values>x_plotValues[1]]
+#
+#
+# # gamma = 0
+# x_rest = X_Values[X_Values>x_plotValues[3]]
+# Y_Values = np.array(Y_Values) #convert the np array
+# y_rest = Y_Values[X_Values>x_plotValues[3]]
+
+# gamma between
+# Y_Values = np.array(Y_Values) #convert the np array
+# X_Values = np.array(X_Values) #convert the np array
+#
+# x_one = X_Values[X_Values>x_plotValues[3]]
+# # ax.scatter(X_Values, Y_Values)
+# y_rest = Y_Values[X_Values>x_plotValues[3]]
+# ax.plot(X_Values[X_Values>0.135], Y_Values[X_Values<0.135])
+#
+#
+#
+
+
+# y_rest = Y_Values[np.nonzero(X_Values>x_plotValues[1]]
+# print('X_Values:', X_Values)
+# print('Y_Values:', Y_Values)
+# print('x_rest:', x_rest)
+# print('y_rest:', y_rest)
+# print('np.nonzero(X_Values>x_plotValues[1]', np.nonzero(X_Values>x_plotValues[1]) )
+
+
+
+
+# --- Convert to numpy array
+Y_Values = np.array(Y_Values)
+X_Values = np.array(X_Values)
+
+Y_arr = np.asarray(Y_Values, dtype=float)
+X_Values = np.asarray(X_Values, dtype=float)
+
+
+print('X_Values:', X_Values)
+print('Y_arr:', Y_arr)
+# ---------------- Create Plot -------------------
+
+#--- change plot style: SEABORN
+# plt.style.use("seaborn-paper")
+
+
+#--- Adjust gobal matplotlib variables
+# mpl.rcParams['pdf.fonttype'] = 42
+# mpl.rcParams['ps.fonttype'] = 42
+mpl.rcParams['text.usetex'] = True
+mpl.rcParams["font.family"] = "serif"
+mpl.rcParams["font.size"] = "9"
+
+
+# plt.rc('font', family='serif', serif='Times')
+# plt.rc('font', family='serif')
+# # plt.rc('text', usetex=True) #also works...
+# plt.rc('xtick', labelsize=8)
+# plt.rc('ytick', labelsize=8)
+# plt.rc('axes', labelsize=8)
+
+
+
+
+
+#---- Scale Figure apropriately to fit tex-File Width
+# width = 452.9679
+
+# width as measured in inkscape
+width = 6.28 *0.5
+height = width / 1.618
+
+#setup canvas first
+fig = plt.figure() #main
+# fig, ax = plt.subplots()
+# fig, (ax, ax2) = plt.subplots(ncols=2)
+# fig,axes = plt.subplots(nrows=1,ncols=2,figsize=(width,height)) # more than one plot
+
+
+# fig.subplots_adjust(left=.15, bottom=.16, right=.99, top=.97) #TEST
+
+
+# TEST
+# mpl.rcParams['figure.figsize'] = (width+0.1,height+0.1)
+# fig = plt.figure(figsize=(width+0.1,height+0.1))
+
+
+# mpl.rcParams['figure.figsize'] = (width,height)
+# fig = plt.figure(figsize=(10,6)) # default is [6.4,4.8] 6.4 is the width, 4.8 is the height
+# fig = plt.figure(figsize=(width,height)) # default is [6.4,4.8] 6.4 is the width, 4.8 is the height
+# fig = plt.figure(figsize=set_size(width))
+# fig = plt.subplots(1, 1, figsize=set_size(width))
+
+# --- To create a figure half the width of your document:#
+# fig = plt.figure(figsize=set_size(width, fraction=0.5))
+
+
+
+#--- You must select the correct size of the plot in advance
+# fig.set_size_inches(3.54,3.54)
+
+ax = plt.axes((0.1,0.1,0.8,0.8))
+# ax = plt.axes((0.1,0.1,0.5,0.8))
+# ax = plt.axes((0.1,0.1,1,1))
+# ax = plt.axes()
+
+# ax.spines['right'].set_visible(False)
+# ax.spines['left'].set_visible(False)
+# ax.spines['bottom'].set_visible(False)
+# ax.spines['top'].set_visible(False)
+# ax.tick_params(axis='x',which='major',direction='out',length=10,width=5,color='red',pad=15,labelsize=15,labelcolor='green',
+# labelrotation=15)
+# ax.tick_params(axis='x',which='major', direction='out',pad=5,labelsize=10)
+# ax.tick_params(axis='y',which='major', length=5, width=1, direction='out',pad=5,labelsize=10)
+
+
+# ax.tick_params(axis='x',which='major', direction='out',pad=3)
+# ax.tick_params(axis='y',which='major', length=3, width=1, direction='out',pad=3)
+# ax.xaxis.set_major_locator(MultipleLocator(0.05))
+# ax.xaxis.set_minor_locator(MultipleLocator(0.025))
+
+
+#---- print data-types
+# print(ax.xaxis.get_major_locator())
+# print(ax.xaxis.get_minor_locator())
+# print(ax.xaxis.get_major_formatter())
+# print(ax.xaxis.get_minor_formatter())
+
+#---- Hide Ticks or Labels
+# ax.yaxis.set_major_locator(plt.NullLocator())
+# ax.xaxis.set_major_formatter(plt.NullFormatter())
+
+#---- Reducing or Increasing the Number of Ticks
+# ax.xaxis.set_major_locator(plt.MaxNLocator(3))
+# ax.yaxis.set_major_locator(plt.MaxNLocator(3))
+
+
+#----- Fancy Tick Formats
+# ax.yaxis.set_major_locator(plt.MultipleLocator(np.pi / 4))
+# ax.yaxis.set_minor_locator(plt.MultipleLocator(np.pi / 12))
+#
+#
+# # ax.set_yticks([0, np.pi/8, np.pi/4 ])
+#
+# ax.yaxis.set_major_formatter(plt.FuncFormatter(format_func))
+
+
+
+# --- manually change ticks&labels:
+# ax.set_xticks([0.2,1])
+# ax.set_xticklabels(['pos1','pos2'])
+
+# ax.set_yticks([0, np.pi/8, np.pi/4 ])
+# labels = ['$0$',r'$\pi/8$', r'$\pi/4$']
+# ax.set_yticklabels(labels)
+
+# a=ax.yaxis.get_major_locator()
+# b=ax.yaxis.get_major_formatter()
+# c = ax.get_xticks()
+# d = ax.get_xticklabels()
+# print('xticks:',c)
+# print('xticklabels:',d)
+#
+# ax.grid(True,which='major',axis='both',alpha=0.3)
+
+# ax.plot(Y_arr[:,0], Y_arr[:,1] , 'royalblue')
+
+print('Y_arr[:,0]:', Y_arr[:,0])
+
+print('Y_arr[:,1]:', Y_arr[:,1])
+
+ax.plot(Y_arr[:,0], Y_arr[:,1] , 'royalblue', # data
+marker='o', # each marker will be rendered as a circle
+markersize=2, # marker size
+markerfacecolor='orange', # marker facecolor
+markeredgecolor='black', # marker edgecolor
+markeredgewidth=0.5, # marker edge width
+# linestyle='--', # line style will be dash line
+linewidth=1,
+zorder = 3) # line width
+# plt.figure()
+
+
+#--- Coordinate Axes:
+ax.spines.left.set_position('zero')
+ax.spines.right.set_color('none')
+ax.spines.bottom.set_position('zero')
+ax.spines.top.set_color('none')
+ax.xaxis.set_ticks_position('bottom')
+ax.yaxis.set_ticks_position('left')
+
+ax.set(xlim=(-25, 15), ylim=(-3, 3))
+
+#-- Decorate the spins
+arrow_length = 8 # In points
+# X-axis arrow
+ax.annotate('x', xy=(1, 0), xycoords=('axes fraction', 'data'),
+ xytext=(arrow_length, 0), textcoords='offset points',
+ ha='left', va='center',
+ arrowprops=dict(arrowstyle='<|-', fc='black'))
+
+# Y-axis arrow
+ax.annotate('y', xy=(0, 1), xycoords=('data', 'axes fraction'),
+ xytext=(0, arrow_length), textcoords='offset points',
+ ha='center', va='bottom',
+ arrowprops=dict(arrowstyle='<|-', fc='black'))
+
+
+
+# ax.scatter(Y_arr[21,0],Y_arr[21,1], s=6, marker='o', cmap=None, norm=None, facecolor = 'forestgreen',
+# edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=None, zorder=5)
+
+# ax.text(Y_arr[21,0]-0.25 , Y_arr[21,1]+0.15, r"$1$", size=4, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5))
+# ax.text(Y_arr[21,0] , Y_arr[21,1], r"$1$", size=2, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5), zorder=5)
+
+ax.scatter(Y_arr[21,0] , Y_arr[21,1], s=4, marker='o', cmap=None, norm=None, facecolor = 'forestgreen',
+ edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=0.5, zorder=5)
+
+ax.scatter(Y_arr[31,0] , Y_arr[31,1], s=4, marker='o', cmap=None, norm=None, facecolor = 'forestgreen',
+ edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=0.5, zorder=5)
+
+ax.scatter(Y_arr[40,0] , Y_arr[40,1], s=4, marker='o', cmap=None, norm=None, facecolor = 'forestgreen',
+ edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=0.5, zorder=5)
+
+ax.annotate( 1 , (Y_arr[21,0] , Y_arr[21,1]), xytext=(Y_arr[21,0]-0.25 , Y_arr[21,1]+0.5),
+ bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5),
+ arrowprops = dict(arrowstyle="simple",color='blue', linewidth=0.1), fontsize=4)
+
+ax.annotate( 2 , (Y_arr[31,0] , Y_arr[31,1]), xytext=(Y_arr[31,0]+2 , Y_arr[31,1]-0.05),
+ bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5),
+ arrowprops = dict(arrowstyle="simple",color='blue', linewidth=0.1), fontsize=4)
+
+ax.annotate( 3 , (Y_arr[40,0] , Y_arr[40,1]), xytext=(Y_arr[40,0]-0.25 , Y_arr[40,1]+0.5),
+ bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5),
+ arrowprops = dict(arrowstyle="simple",color='blue', linewidth=0.1), fontsize=4)
+
+
+ # arrowprops = dict(arrowstyle="simple",color='blue', linewidth=0.1, shrink=0.05), fontsize=4)
+# f,ax=plt.subplots(1)
+
+# plt.title(r''+ yName + '-Plot')
+# plt.plot(X_Values, Y_Values,linewidth=2, '.k')
+# plt.plot(X_Values, Y_Values,'.k',markersize=1)
+# plt.plot(X_Values, Y_Values,'.',markersize=0.8)
+
+# plt.plot(X_Values, Y_Values)
+
+# ax.plot([[0],X_Values[-1]], [Y_Values[0],Y_Values[-1]])
+
+
+
+# Gamma = '0'
+# ax.plot([x_plotValues[0],x_plotValues[1]], [y_plotValues[0],y_plotValues[1]] , 'b')
+#
+# ax.plot([x_plotValues[1],x_plotValues[3]], [y_plotValues[2],y_plotValues[3]] , 'b')
+#
+# ax.plot(x_rest, y_rest, 'b')
+
+
+# Gamma between
+
+# x jump values (gamma 0): [0.13606060606060608, 0.21090909090909093]
+
+# ax.plot([[0,jump_xValues[0]], [0, 0]] , 'b')
+# ax.plot([jump_xValues[0],xmin], [y_plotValues[2],y_plotValues[2]] , 'b')
+
+# ax.plot([[0,0.13606060606060608], [0, 0]] , 'b')
+# ax.plot([[0.13606060606060608,xmin], [(math.pi/2),(math.pi/2)]], 'b')
+
+# jump_xValues[0]
+
+
+
+# --- leave out jumps:
+# ax.scatter(X_Values, Y_Values)
+
+
+
+
+# # --- leave out jumps:
+# if gamma == 'infinity':
+# ax.plot(X_Values[X_Values>=jump_xValues[0]], Y_Values[X_Values>=jump_xValues[0]] , 'royalblue')
+# ax.plot(X_Values[X_Values<jump_xValues[0]], Y_Values[X_Values<jump_xValues[0]], 'royalblue')
+# # ax.plot(X_Values[X_Values>=jump_xValues[0]], Y_Values[X_Values>=jump_xValues[0]])
+# # ax.plot(X_Values[X_Values<jump_xValues[0]], Y_Values[X_Values<jump_xValues[0]])
+
+
+
+
+# ax.plot(X_Values[X_Values>0.136], Y_Values[X_Values>0.136])
+# ax.plot(X_Values[X_Values<0.135], Y_Values[X_Values<0.135])
+# ax.scatter(X_Values, Y_Values)
+# ax.plot(X_Values, Y_Values)
+
+# plt.plot(x_plotValues, y_plotValues,'.')
+# plt.scatter(X_Values, Y_Values, alpha=0.3)
+# plt.scatter(X_Values, Y_Values)
+# plt.plot(X_Values, Y_Values,'.')
+# plt.plot([X_Values[0],X_Values[-1]], [Y_Values[0],Y_Values[-1]])
+# plt.axis([0, 6, 0, 20])
+
+# ax.set_xlabel(r"volume fraction $\theta$", size=11)
+# ax.set_ylabel(r"angle $\angle$", size=11)
+# ax.set_xlabel(r"volume fraction $\theta$")
+# ax.set_ylabel(r"angle $\angle$")
+# ax.set_ylabel(r"$a^*$")
+# plt.ylabel('$\kappa$')
+
+# ax.yaxis.set_major_formatter(ticker.FormatStrFormatter('%g $\pi$'))
+# ax.yaxis.set_major_locator(ticker.MultipleLocator(base=0.1))
+
+
+
+
+# Plot every other line.. not the jumps...
+
+# if gamma == '0':
+# tmp = 1
+# for idx, x in enumerate(x_plotValues):
+# if idx > 0 and tmp == 1:
+# # plt.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]] )
+# ax.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]], 'royalblue' )
+# tmp = 0
+# else:
+# tmp = 1
+
+
+
+
+# for x in jump_xValues:
+# plt.axvline(x,ymin=0, ymax= 1, color = 'orange',alpha=0.5, linestyle = 'dashed', linewidth=1)
+
+
+
+ # plt.axvline(x,ymin=0, ymax= 1, color = 'orange',alpha=0.5, linestyle = 'dashed', label=r'$\theta_*$')
+
+# plt.axvline(x_plotValues[1],ymin=0, ymax= 1, color = 'g',alpha=0.5, linestyle = 'dashed')
+
+# plt.axhline(y = 1.90476, color = 'b', linestyle = ':', label='$q_1$')
+# plt.axhline(y = 2.08333, color = 'r', linestyle = 'dashed', label='$q_2$')
+# plt.legend()
+
+
+# -- SETUP LEGEND
+# ax.legend(prop={'size': 11})
+# ax.legend()
+
+# ------------------ SAVE FIGURE
+# tikzplotlib.save("TesTout.tex")
+# plt.close()
+# mpl.rcParams.update(mpl.rcParamsDefault)
+
+# plt.savefig("graph.pdf",
+# #This is simple recomendation for publication plots
+# dpi=1000,
+# # Plot will be occupy a maximum of available space
+# bbox_inches='tight',
+# )
+# plt.savefig("graph.pdf")
+
+
+#
+# # Find transition point
+# lastIdx = len(Y_Values)-1
+#
+# for idx, y in enumerate(Y_Values):
+# if idx != lastIdx:
+# if abs(y-0) < 0.01 and abs(Y_Values[idx+1] - 0) > 0.05:
+# transition_point1 = X_Values[idx+1]
+# print('transition point1:', transition_point1 )
+# if abs(y-0.5*np.pi) < 0.01 and abs(Y_Values[idx+1] -0.5*np.pi)>0.01:
+# transition_point2 = X_Values[idx]
+# print('transition point2:', transition_point2 )
+# if abs(y-0) > 0.01 and abs(Y_Values[idx+1] - 0) < 0.01:
+# transition_point3 = X_Values[idx+1]
+# print('transition point3:', transition_point3 )
+#
+# # Add transition Points:
+# if gamma == '0':
+# ax.scatter([transition_point1, transition_point2],[np.pi/2,np.pi/2],s=6, marker='o', cmap=None, norm=None, facecolor = 'black',
+# edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=None, zorder=3)
+#
+# ax.text(transition_point1-0.02 , np.pi/2-0.02, r"$1$", size=6, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5)
+# )
+#
+# ax.text(transition_point2+0.012 , np.pi/2-0.02, r"$2$", size=6, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5)
+# )
+# else:
+# ax.scatter([transition_point1, transition_point2, transition_point3 ],[np.pi/2,np.pi/2,0 ],s=6, marker='o', cmap=None, norm=None, facecolor = 'black',
+# edgecolor = 'black', vmin=None, vmax=None, alpha=None, linewidths=None, zorder=3)
+#
+# ax.text(transition_point1-0.02 , np.pi/2-0.02, r"$1$", size=6, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5)
+# )
+#
+# ax.text(transition_point2 +0.015 , np.pi/2-0.02, r"$2$", size=6, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5)
+# )
+#
+# ax.text(transition_point3 +0.005 , 0+0.06, r"$3$", size=6, bbox=dict(boxstyle="circle",facecolor='white', alpha=1.0, pad=0.1, linewidth=0.5)
+# )
+
+
+fig.set_size_inches(width, height)
+fig.savefig('Plot-aStar.pdf')
+
+
+
+
+# tikz_save('someplot.tex', figureheight='5cm', figurewidth='9cm')
+
+# tikz_save('fig.tikz',
+# figureheight = '\\figureheight',
+# figurewidth = '\\figurewidth')
+
+# ----------------------------------------
+
+
+plt.show()
+# #---------------------------------------------------------------
diff --git a/src/plot-q3-gamma.py b/src/plot-q3-gamma.py
index 616f2f6723fbd136cff5bb5c5aa3cce184f34f99..daad8f02254c785e0fb89a310a3a4d501cc00251 100644
--- a/src/plot-q3-gamma.py
+++ b/src/plot-q3-gamma.py
@@ -52,7 +52,7 @@ print("Path: ", path)
#---------------------------------------------------------------
print('---- Input parameters: -----')
-mu1 = 10.0
+mu1 = 1.0
rho1 = 1.0
# alpha = 2.8
alpha = 2.0
@@ -128,7 +128,8 @@ width = 6.28 *0.5
height = width / 1.618
fig = plt.figure()
# ax = plt.axes((0.15,0.18,0.8,0.8))
-ax = plt.axes((0.15,0.21 ,0.8,0.75))
+# ax = plt.axes((0.15,0.21 ,0.8,0.75))
+ax = plt.axes((0.18,0.21 ,0.8,0.75))
ax.tick_params(axis='x',which='major', direction='out',pad=3)
ax.tick_params(axis='y',which='major', length=3, width=1, direction='out',pad=3)
ax.xaxis.set_major_locator(MultipleLocator(0.5))
@@ -162,8 +163,10 @@ ax.set_ylabel(r"$q_3(\gamma)$")
# plt.axhline(y = 1.90476, color = 'forestgreen', linestyle = ':', label='$q_1$')
# plt.axhline(y = 2.08333, color = 'darkorange', linestyle = 'dashed', label='$q_2$')
-plt.axhline(y = q1, color = 'forestgreen', linestyle = 'dashed', label='$q_1$')
-plt.axhline(y = q2, color = 'darkorange', linestyle = 'dashed', label='$q_2$')
+# plt.axhline(y = q1, color = 'forestgreen', linestyle = 'dashed', label='$q_1$')
+# plt.axhline(y = q2, color = 'darkorange', linestyle = 'dashed', label='$q_2$')
+plt.axhline(y = q1, color = 'forestgreen', linestyle = 'dashed', label=r'$\langle \mu \rangle $')
+plt.axhline(y = q2, color = 'darkorange', linestyle = 'dashed', label='$\overline{\mu}$')
ax.legend(loc = 'center right')
# plt.legend()