Step 21: Plot the dynamic results

Now use the programs and methods previously presented in Step 14 to display maps of every aspect of your dynamic solution. Files used for plotting will include the finite element grid (.feg) file and input parameter (.in) file as before, and the dynamic aspects will be computed from the v_____.out file (nodal velocities) and f_____.out file (nodal reaction forces) produced by Shells in Step 20. In particular, the strain rates will be computed by combining the information in the .feg and v_____.out files, and then the stresses will be computed by using the rheologies from the .in file. Therefore, it is important that the .in file used for plotting should be the same as the one that was used to run Shells. (Personally, I like to keep a separate input parameter file for each dynamic experiment.)

New maps available:

 surface velocity

 velocity discontinuities across faults (proportional to fault slip rates)

 strain rates (in the continuum elements between the faults)

 vertically-integrated stress anomaly tensors

 most-compressive horizontal principal stress directions

 shear tractions on the base of the lithosphere
(not shown for this model because they were uniformly zero)

 forces on boundary nodes (from boundary conditions)