OOF.Microstructure.Create_From_ImageFile(filename='MeshSubject_02_2dec.pgm', microstructure_name='MeshSubject_02_2dec.pgm', height=automatic, width=automatic)
OOF.Windows.Graphics.New()
OOF.Graphics_1.Settings.Zoom.Fill_Window()
OOF.File.Load.Data(filename='/Users/ccarter/People/Neri/MeshSubject_02/black.mat')
OOF.File.Load.Data(filename='/Users/ccarter/People/Neri/MeshSubject_02/gray.mat')
OOF.File.Load.Data(filename='/Users/ccarter/People/Neri/MeshSubject_02/white.mat')
OOF.PixelGroup.New(name='black', microstructure='MeshSubject_02_2dec.pgm')
OOF.PixelGroup.New(name='white', microstructure='MeshSubject_02_2dec.pgm')
OOF.PixelGroup.New(name='gray', microstructure='MeshSubject_02_2dec.pgm')
OOF.Graphics_1.Toolbox.Pixel_Select.Color(source='MeshSubject_02_2dec.pgm:MeshSubject_02_2dec.pgm', range=DeltaGray(delta_gray=0.29752099999999998), points=[Point(85.2906,662.235)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Pixel_Select.Color(source='MeshSubject_02_2dec.pgm:MeshSubject_02_2dec.pgm', range=DeltaGray(delta_gray=0.29752099999999998), points=[Point(150.131,634.447)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Pixel_Select.Color(source='MeshSubject_02_2dec.pgm:MeshSubject_02_2dec.pgm', range=DeltaGray(delta_gray=0.29752099999999998), points=[Point(85.2906,658.265)], shift=0, ctrl=0)
OOF.PixelGroup.AddSelection(microstructure='MeshSubject_02_2dec.pgm', group='white')
OOF.Graphics_1.Toolbox.Pixel_Select.Color(source='MeshSubject_02_2dec.pgm:MeshSubject_02_2dec.pgm', range=DeltaGray(delta_gray=0.29752099999999998), points=[Point(142.191,580.193)], shift=0, ctrl=0)
OOF.PixelGroup.AddSelection(microstructure='MeshSubject_02_2dec.pgm', group='black')
OOF.PixelSelection.Add_Group(microstructure='MeshSubject_02_2dec.pgm', group='white')
OOF.PixelSelection.Invert(microstructure='MeshSubject_02_2dec.pgm')
OOF.PixelGroup.AddSelection(microstructure='MeshSubject_02_2dec.pgm', group='gray')
OOF.Material.Assign(material='white', microstructure='MeshSubject_02_2dec.pgm', pixels='white')
OOF.Material.Assign(material='black', microstructure='MeshSubject_02_2dec.pgm', pixels='black')
OOF.Material.Assign(material='gray', microstructure='MeshSubject_02_2dec.pgm', pixels='gray')
OOF.Skeleton.New(name='skeleton', microstructure='MeshSubject_02_2dec.pgm', x_elements=12, y_elements=12, skeleton_geometry=QuadSkeleton(left_right_periodicity=False,top_bottom_periodicity=False))
OOF.Graphics_1.Layer.Select(n=4)
OOF.LayerEditor.LayerSet.Edit(window='Graphics_1', layer_number=4)
OOF.LayerEditor.LayerSet.Replace_Method(layer_number=0, method=SkeletonMaterialDisplay())
OOF.LayerEditor.LayerSet.Send(window='Graphics_1')
OOF.Graphics_1.Toolbox.Pixel_Select.Clear(source='MeshSubject_02_2dec.pgm:MeshSubject_02_2dec.pgm')
OOF.LayerEditor.LayerSet.Edit(window='Graphics_1', layer_number=1)
OOF.LayerEditor.LayerSet.Replace_Method(layer_number=0, method=SkeletonEdgeDisplay(color=HSVColor(hue=114.286,saturation=0.49206299999999997,value=0.57142899999999996),width=0))
OOF.LayerEditor.LayerSet.Send(window='Graphics_1')
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(33.6832,457.129)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(28.3901,402.875)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(95.8767,394.935)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(37.653,323.479)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(89.2604,328.772)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(89.2604,285.104)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(40.2995,291.72)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(31.0366,233.497)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(32.3599,158.07)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(85.2906,106.463)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(160.717,40.2995)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(144.838,94.5535)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(99.8465,42.946)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(36.3297,103.816)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(221.587,93.2302)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(527.262,660.912)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(578.869,712.519)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(627.83,715.166)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(581.516,646.356)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(631.8,649.002)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(529.908,586.809)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(572.253,590.779)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(629.153,590.779)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(377.733,353.914)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(336.711,351.267)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(97.2,693.994)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(98.5233,646.356)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(154.1,649.002)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(155.424,704.58)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(45.5926,696.64)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(35.0064,643.709)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(32.3599,586.809)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(201.738,332.742)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(384.349,106.463)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(574.9,95.8767)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(506.09,94.5535)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(619.891,95.8767)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(639.74,97.2)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(692.67,97.2)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(508.736,201.738)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(531.232,142.191)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(454.482,205.708)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(593.425,332.742)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(581.516,208.354)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(639.74,273.195)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(527.262,443.896)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(457.129,450.512)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(384.349,568.283)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(332.742,643.709)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(144.838,516.676)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(150.131,469.038)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(203.061,471.685)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(213.648,519.322)], shift=1, ctrl=0)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SnapNodes(targets=SnapAll(),criterion=AverageEnergy(alpha=0.29999999999999999)))
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.90000000000000002)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.39583299999999999)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SnapNodes(targets=SnapAll(),criterion=AverageEnergy(alpha=0.29999999999999999)))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SplitQuads(targets=AllElements(),criterion=AverageEnergy(alpha=0.29999999999999999),split_how=GeographicQ2T()))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=MergeTriangles(targets=AllElements(),criterion=AverageEnergy(alpha=0.29999999999999999)))
OOF.Skeleton.PinNodes.Pin_Internal_Boundary_Nodes(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Anneal(targets=AllNodes(),criterion=AverageEnergy(alpha=0.29999999999999999),T=0.0,delta=1.0,iteration=ConditionalIteration(condition=AcceptanceRate(acceptance_rate=3),extra=5,maximum=100)))
OOF.Skeleton.PinNodes.UnpinAll(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.39583299999999999)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.52083299999999999)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.42708299999999999)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.46875)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SnapNodes(targets=SnapAll(),criterion=AverageEnergy(alpha=0.29999999999999999)))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SplitQuads(targets=AllElements(),criterion=AverageEnergy(alpha=0.29999999999999999),split_how=GeographicQ2T()))
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(156.747,589.455)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(171.303,586.809)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(154.1,565.637)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(9.86436,519.322)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(9.86436,417.431)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(28.3901,420.077)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(29.7134,398.905)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(31.0366,375.086)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(52.2089,376.409)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(36.3297,352.591)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(49.5624,355.237)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(35.0064,336.711)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(46.9158,336.711)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(48.2391,193.799)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(72.0579,200.415)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(69.4114,189.829)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(77.351,191.152)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(73.3812,193.799)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(77.351,191.152)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Undo(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Undo(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(328.772,44.2693)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(683.407,352.591)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(265.255,33.6832)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(273.195,266.578)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(380.379,222.91)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(409.491,213.648)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(332.742,199.092)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(318.186,199.092)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(331.418,226.88)], shift=1, ctrl=0)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SnapNodes(targets=SnapAll(),criterion=AverageEnergy(alpha=0.29999999999999999)))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=SplitQuads(targets=AllElements(),criterion=AverageEnergy(alpha=0.29999999999999999),split_how=GeographicQ2T()))
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckHomogeneity(threshold=0.90000000000000002),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.LayerEditor.LayerSet.Edit(window='Graphics_1', layer_number=4)
OOF.LayerEditor.LayerSet.Add_Method(method=SkeletonMaterialDisplay())
OOF.LayerEditor.LayerSet.Send(window='Graphics_1')
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.46875)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.70833299999999999)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.114583)
OOF.Skeleton.PinNodes.Pin_Internal_Boundary_Nodes(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(253.346,576.223)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(253.346,577.546)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(274.518,596.072)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(291.72,580.193)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(274.518,561.667)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(254.669,578.869)], shift=1, ctrl=0)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Skeleton.PinNodes.UnpinAll(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(470.361,189.829)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(486.241,187.182)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(482.271,179.243)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(475.654,195.122)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(470.361,191.152)], shift=1, ctrl=0)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(273.195,398.905)], shift=0, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(273.195,416.107)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(261.285,319.509)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(486.241,700.61)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(169.98,327.449)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(169.98,322.155)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(605.335,499.473)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(606.658,529.908)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(597.395,520.646)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(600.042,490.21)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(604.011,482.271)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Undo(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(611.951,484.917)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(585.486,348.621)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(589.455,348.621)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(687.377,205.708)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(699.287,193.799)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(622.537,27.0668)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(630.477,25.7436)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(53.5322,498.15)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(267.901,548.434)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(255.992,541.818)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(205.708,676.791)], shift=1, ctrl=0)
OOF.Graphics_1.Toolbox.Select_Element.Single_Element(skeleton='MeshSubject_02_2dec.pgm:skeleton', points=[Point(213.648,683.407)], shift=1, ctrl=0)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=Unconditionally(),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.Graphics_1.Toolbox.Select_Element.Clear(skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.File.Save.Skeleton(filename='/Users/ccarter/People/Neri/MeshSubject_02/15481.skel', mode='w', format='ascii', skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.28125)
OOF.ElementSelection.Select_by_Homogeneity(skeleton='MeshSubject_02_2dec.pgm:skeleton', threshold=0.70833299999999999)
OOF.Skeleton.Modify(skeleton='MeshSubject_02_2dec.pgm:skeleton', modifier=Refine(targets=CheckSelectedElements(),criterion=MinimumArea(threshold=24,units='Pixel'),degree=Trisection(rule_set='conservative'),alpha=0.29999999999999999))
OOF.File.Save.Skeleton(filename='/Users/ccarter/People/Neri/MeshSubject_02/33676.skel', mode='w', format='ascii', skeleton='MeshSubject_02_2dec.pgm:skeleton')
OOF.LayerEditor.LayerSet.Edit(window='Graphics_1', layer_number=5)
OOF.Graphics_1.Layer.Hide(n=5)
OOF.Graphics_1.Layer.Hide(n=1)
OOF.Graphics_1.Layer.Select(n=1)
OOF.Mesh.New(name='mesh', skeleton='MeshSubject_02_2dec.pgm:skeleton', element_types=['T3_3', 'Q4_4'])
OOF.Graphics_1.Layer.Hide(n=6)
OOF.Graphics_1.Layer.Select(n=6)
OOF.Subproblem.Field.Define(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Subproblem.Field.Define(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Heat_Eqn)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Heat_Flux)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Mesh.Boundary_Conditions.New(name='bc', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Temperature,field_component='',equation=Heat_Eqn,eqn_component='',profile=ConstantProfile(value=0.0),boundary='left'))
OOF.Mesh.Boundary_Conditions.New(name='bc<2>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Temperature,field_component='',equation=Heat_Eqn,eqn_component='',profile=ConstantProfile(value=1.0),boundary='right'))
OOF.Mesh.Boundary_Conditions.New(name='bc<3>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Displacement,field_component='x',equation=Force_Balance,eqn_component='x',profile=ConstantProfile(value=1.0),boundary='right'))
OOF.Mesh.Boundary_Conditions.New(name='bc<4>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Displacement,field_component='x',equation=Force_Balance,eqn_component='x',profile=ConstantProfile(value=0),boundary='left'))
OOF.Mesh.Boundary_Conditions.Edit(name='bc<3>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Displacement,field_component='x',equation=Force_Balance,eqn_component='x',profile=ConstantProfile(value=3.0),boundary='right'))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=5000,krylov_dimension=500,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Graphics_1.Layer.Show(n=6)
OOF.Graphics_1.Layer.Hide(n=6)
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=200,krylov_dimension=500,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=1000,krylov_dimension=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.File.Save.Mesh(filename='/Users/ccarter/People/Neri/MeshSubject_02/horizontal_loads/33676_converged.msh', mode='w', format='ascii', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Component',component='x',flux=Heat_Flux), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/jx_horiz.csv', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Invariant',invariant=Magnitude(),flux=Heat_Flux), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/jmag_horiz.csv', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Invariant',invariant=MatrixTrace(),flux=Stress), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/stress_horiz.csv', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Field:Component',component='z',field=Displacement_z), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/z_horiz.csv', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Energy',etype='Total'), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/energy_horiz.csv', comment_character='#', separator='comma')
OOF.LayerEditor.LayerSet.Edit(window='Graphics_1', layer_number=6)
OOF.LayerEditor.LayerSet.Replace_Method(layer_number=0, method=FilledContourDisplay(what=getOutput('Flux:Invariant',invariant=MatrixTrace(),flux=Stress),where=getOutput('original'),min=automatic,max=automatic,levels=11,nbins=5,colormap=ThermalMap()))
OOF.LayerEditor.LayerSet.Send(window='Graphics_1')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Strain:Invariant',invariant=MatrixTrace(),type=GeometricStrain()), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/horizontal_loads/strain_horiz.csv', comment_character='#', separator='comma')
OOF.Graphics_1.Layer.Hide(n=2)
OOF.Graphics_1.File.Close()