/42TestingPoints/SinglePhysics/LinearProblems/DynamicProblems/Stokes/3D/HEX/CubicVelocityQuadraticPressure/LEVEL_1/output/visualise.com

 1#Read in the sequence of nodal positions.
 2for $i (1..3)
 3  {
 4	 $filename = sprintf("./output/TIME_STEP_%04d.exnode", $i);
 5	 print "Reading $filename time $i\n";
 6	 gfx read node "$filename" time $i;
 7  }
 8
 9#Read in the element description
10gfx read elem ./output/TIME_STEP_0001.exelem;
11gfx create window 1
12
13gfx cre spectrum flow
14gfx modify spectrum flow clear overwrite_colour;
15gfx modify spectrum flow linear reverse range 1 1.25 extend_above extend_below rainbow colour_range 0 1 component 2;
16
17gfx def field vector_field coord rectangular_cartesian component general.1 general.2 general.3
18
19gfx cre spectrum pressure
20gfx modify spectrum pressure linear reverse range -1 1 extend_above extend_below rainbow colour_range 0 1 component 4;
21
22gfx modify g_element OpenCMISS node_points glyph arrow_solid general size "1.0*1.0*1.0" centre 0,0,0 select_on material default selected_material default_selected data vector_field orientation vector_field scale_factors "1.0*1.0*1.0" spectrum flow
23
24gfx modify window 1 background colour 1 1 1
25
26gfx define faces egroup OpenCMISS
27gfx modify g_element OpenCMISS surfaces select_on material default selected_material default_selected data general spectrum pressure
28
29#Set the timekeeper playing
30gfx timekeeper default play speed 1 skip;
31gfx create time_editor
32
33gfx edit scene
34gfx edit spectrum