Changes between Version 80 and Version 81 of VisIt

Timestamp:

11/10/11 10:33:28 (13 years ago)

Author:

Jonathan

Comment:

—

VisIt

@@ -269 +269 @@
 }}}
 [[CollapsibleEnd]]
+[[CollapsibleStart(Script for following collapse of uniform sphere)]]
+This script adjusts the window bounds to zoom in on a region at the same rate the region collapses (assuming it follows the solution for a uniform density sphere undergoing free fall collapse)
+{{{
+r_0=60                 # Initial radius of view (usually 150% of clump radius is good)
+r_min=8                # Radius to stop shrinking at (32*levels(MaxLevel)%dx)
+singular_time=.66      #t_ff using clump density in computational units (or by eye)
+center=(240,240)       #Clump location in computational units
+starting_frame=0       #Frame to start from
+final_frame=63         #Frame to stop at
+my_dir = "/grassdata/johannjc/scrambler_111007/modules/Problem/Images/2D" #directory for movie output
+movie_name = "ZoomingRhoTestwSinksA" #base name for movie
+
+
+# ~~~~~~~~~~~~~~~~~  Output file attributes...  ~~~~~~~~~~~~~~~~~~~~~~~~~
+OldSaveWindowAtts = SaveWindowAttributes()
+SaveWindowAtts = SaveWindowAttributes()
+SaveWindowAtts.outputToCurrentDirectory = 0
+SaveWindowAtts.outputDirectory = my_dir
+SaveWindowAtts.fileName = movie_name 
+SaveWindowAtts.family = 0
+SaveWindowAtts.format = SaveWindowAtts.PNG  #BMP, CURVE, JPEG, OBJ, PNG, POSTSCRIPT, PPM, RGB, STL, TIFF, ULTRA, VTK
+SaveWindowAtts.width = 1024
+SaveWindowAtts.height = 1024
+SaveWindowAtts.screenCapture = 1
+SaveWindowAtts.saveTiled = 1
+SaveWindowAtts.quality = 80
+SaveWindowAtts.progressive = 0
+SaveWindowAtts.binary = 0
+SaveWindowAtts.stereo = 0
+SaveWindowAtts.compression = SaveWindowAtts.PackBits
+GlobalAtts = GetGlobalAttributes()
+
+
+a=GetView2D()
+from math import *
+for state in range(starting_frame, final_frame):
+  Query('Time')   #Get current database time
+  tratio=min(GetQueryOutputValue()/singular_time*1.0,1.0) # Calculate t/t_ff
+  rratio=.7178*sqrt(1-tratio**2)+.3068*(1-tratio**2)-.01842 # Approximate fit to free fall radius ratio
+  r=max(r_min,r_0*rratio) # Scale current radius based on rratio keeping it > r_min
+  a.windowCoords=(center[0]-r, center[0]+r, center[1]-r, center[1]+r) # Set window coords to encompass clump circle
+  for wins in GlobalAtts.windows: #In case there are multiple windows
+    SetActiveWindow(wins)
+    SetTimeSliderState(state)
+    SetView2D(a)
+  SaveWindowAtts.fileName = movie_name + "%(frame)04d" % {'frame':state}
+  SetSaveWindowAttributes(SaveWindowAtts)
+  SaveWindow()
+  ClearCache("localhost")
+SetSaveWindowAttributes(OldSaveWindowAtts)
+}}}
+
 [[CollapsibleStart(More complicated movie tracking a particle)]]
 I first had to use some linux-fu to get the particle's xy location that I wanted to track into xloc.dat.  I then used {{{d=f.readline}}} to read each line and {{{c=d.split}}} to parse the line.  I then used {{{c[3]}}} and {{{c[5]}}} as the x and y position of the particle which I then used to set the window coords etc...  Here is an AVI version of the [attachment:sinkzoom2.AVI movie]