Changes between Version 82 and Version 83 of VisIt

Timestamp:

10/23/12 12:43:37 (14 years ago)

Author:

Erica Kaminski

Comment:

—

VisIt

@@ -95 +95 @@
 
 # Redraw
+DrawPlots()
+}}}
+[[CollapsibleEnd]]
+[[CollapsibleStart(Average quantities over concentric spherical shells)]]
+{{{
+ActivateDatabase("localhost:/media/grassdata/erica/from_bluehive/FinalRunsNewest/BP/chombo00046.hdf")
+filename = "average_vrad_00046.curve"
+DeleteActivePlots()
+AddPlot("Pseudocolor", "vrad", 1, 1)
+SetActivePlots(0)
+
+# inner_sphere():
+AddOperator("Clip", 1)
+ClipAtts = ClipAttributes()
+ClipAtts.quality = ClipAtts.Fast # Fast, Accurate
+ClipAtts.funcType = ClipAtts.Sphere # Plane, Sphere
+ClipAtts.plane1Status = 1
+ClipAtts.plane2Status = 0
+ClipAtts.plane3Status = 0
+ClipAtts.plane1Origin = (0, 0, 0)
+ClipAtts.plane2Origin = (0, 0, 0)
+ClipAtts.plane3Origin = (0, 0, 0)
+ClipAtts.plane1Normal = (1, 0, 0)
+ClipAtts.plane2Normal = (0, 1, 0)
+ClipAtts.plane3Normal = (0, 0, 1)
+ClipAtts.planeInverse = 0
+ClipAtts.planeToolControlledClipPlane = ClipAtts.Plane1 # None, Plane1, Plane2, Plane3
+ClipAtts.center = (0, 0, 0)
+ClipAtts.radius = 0.06
+ClipAtts.sphereInverse = 1
+SetOperatorOptions(ClipAtts, 1)
+DrawPlots()
+f = open(filename, "w")
+f.write("# average_vrad\n")
+#Query, write average:
+Query("Average Value")
+t2 = GetQueryOutputValue()
+str = "%25.15e %25.15e\n" %(ClipAtts.radius, t2)
+f.write(str)
+RemoveAllOperators(1)
+for i in range(48):
+     dx=0.03
+     #inner sphere:
+     AddOperator("Clip", 1)
+     ClipAtts = ClipAttributes()
+     ClipAtts.quality = ClipAtts.Fast # Fast, Accurate
+     ClipAtts.funcType = ClipAtts.Sphere # Plane, Sphere
+     ClipAtts.plane1Status = 1
+     ClipAtts.plane2Status = 0
+     ClipAtts.plane3Status = 0
+     ClipAtts.plane1Origin = (0, 0, 0)
+     ClipAtts.plane2Origin = (0, 0, 0)
+     ClipAtts.plane3Origin = (0, 0, 0)
+     ClipAtts.plane1Normal = (1, 0, 0)
+     ClipAtts.plane2Normal = (0, 1, 0)
+     ClipAtts.plane3Normal = (0, 0, 1)
+     ClipAtts.planeInverse = 0
+     ClipAtts.planeToolControlledClipPlane = ClipAtts.Plane1 # None, Plane1, Plane2, Plane3
+     ClipAtts.center = (0, 0, 0)
+     ClipAtts.radius = 0.06 + i*dx
+     ClipAtts.sphereInverse = 0
+     SetOperatorOptions(ClipAtts, 1)
+     #Outer-sphere:
+     AddOperator("Clip", 1)
+     ClipAtts = ClipAttributes()
+     ClipAtts.quality = ClipAtts.Fast # Fast, Accurate
+     ClipAtts.funcType = ClipAtts.Sphere # Plane, Sphere
+     ClipAtts.plane1Status = 1
+     ClipAtts.plane2Status = 0
+     ClipAtts.plane3Status = 0
+     ClipAtts.plane1Origin = (0, 0, 0)
+     ClipAtts.plane2Origin = (0, 0, 0)
+     ClipAtts.plane3Origin = (0, 0, 0)
+     ClipAtts.plane1Normal = (1, 0, 0)
+     ClipAtts.plane2Normal = (0, 1, 0)
+     ClipAtts.plane3Normal = (0, 0, 1)
+     ClipAtts.planeInverse = 0
+     ClipAtts.planeToolControlledClipPlane = ClipAtts.Plane1 # None, Plane1, Plane2, Plane3
+     ClipAtts.center = (0, 0, 0)
+     ClipAtts.radius = 0.06 + (i+1)*dx
+     ClipAtts.sphereInverse = 1
+     SetOperatorOptions(ClipAtts, 1)
+     DrawPlots()
+     Query("Average Value")
+     t2 = GetQueryOutputValue()
+     str = "%25.15e %25.15e\n" %(ClipAtts.radius, t2)
+     #Not sure how to differentiate clips..., would like the outer radius
+     f.write(str)
+     RemoveAllOperators(1)
+f.close()
+DeleteAllPlots()
+OpenDatabase(filename)
+AddPlot("Curve", "average_vrad")
 DrawPlots()
 }}}