#Read in the sequence of nodal positions.
for $i (1..3)
  {
	 $filename = sprintf("./output/TIME_STEP_%04d.exnode", $i);
	 print "Reading $filename time $i\n";
	 gfx read node "$filename" time $i;
  }

#Read in the element description
gfx read elem ./output/TIME_STEP_0001.exelem;
gfx create window 1

gfx cre spectrum flow
gfx modify spectrum flow clear overwrite_colour;
gfx modify spectrum flow linear reverse range 1 1.25 extend_above extend_below rainbow colour_range 0 1 component 2;

gfx def field vector_field coord rectangular_cartesian component general.1 general.2

gfx cre spectrum pressure
gfx modify spectrum pressure linear reverse range 5 15 extend_above extend_below rainbow colour_range 0 1 component 4;

gfx 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

gfx modify window 1 background colour 1 1 1

gfx define faces egroup OpenCMISS
gfx modify g_element OpenCMISS surfaces select_on material default selected_material default_selected data general spectrum pressure

#Set the timekeeper playing
gfx timekeeper default play speed 1 skip;
gfx create time_editor

gfx edit scene
gfx edit spectrum
gfx modify window 1 image scene default light_model default;
gfx modify window 1 image add_light default;
gfx modify window 1 layout 2d ortho_axes z -y eye_spacing 0.25 width 512 height 512;
gfx modify window 1 set current_pane 1;
gfx modify window 1 background colour 1 1 1 texture none;
gfx modify window 1 view parallel eye_point 5 5 28.8485 interest_point 5 5 0 up_vector -0 1 -0 view_angle 40 near_clipping_plane 0.288485 far_clipping_plane 103.095 relative_viewport ndc_placement -1 1 2 2 viewport_coordinates 0 0 1 1;
gfx modify window 1 overlay scene none;
gfx modify window 1 set transform_tool current_pane 1 std_view_angle 40 normal_lines no_antialias depth_of_field 0.0 fast_transparency blend_normal;