The bottom part of the input files describes the types of output requested. Here is an example:

bottom part of input file |
explanation |

* Run everything else interactively | comment |

Non-interactive grid specification | specify grid to which output is referred below |

plane | let this grid be a plane |

vector | request vector output (note: xyz format is always output) |

* output coordinate system | specify output coodinate system |

global | choose global rather than in-plane |

* output files to create | select whichever output is required (there are many choices, read on) |

elem | relative displacements on elements |

displacements | absolute displacements on inspection grid |

rbr | rigid-body rotations on inspection grid |

* output file suffix | 4 characters or less |

t | output file suffix (format allows no more than 4 characters) |

* Xo,Yo,Zo (Ref. corner),strike,dip,length,width: | specify inspection plane |

0 55 0 90 0 59 59 | see below |

* number of grid points in strike and dip directions: | actually the number of equal divisions along strike and dip |

59 59 | see below |

In this example, results will be written with reference to coordinates that define a horizontal plane at the surface whose origin in global coordinates lies 55 km north of the coordinate origin, and that is 59 km along strike and dip. The plane is then divided into 59 intervals along both strike and dip, so that quantities (in the case, components of displacements and rbr) will be written at 1 km spacings. This will yield 60 results along both strike and dip, for a total of 3,600 x 3 displacement results. Displacements of all subelements will also be written to elem files.

Output files formats can be 'xyz' (column) or 'vector' (matrix) format, and they are referenced to either the inspection grid or the elements themselves.

'xyz' format lists each component of the calculated quantity in column format next to its location in either global or in-plane coordinates.

The first three columns are the x, y, and z locations in global coordinates (it's usually better or easier to use global coordinates), the next columns are deformation quantities calculated at this location.

The equivalent vector output comprises three matrices, one for each component of displacement (or whatever quantity is being output): disp1, disp2, and disp3, where 1, 2, and 3 are cryptic mnemonics for x, y, and z. Locations are implicit in the position of the matrix elements. Vector output is useful for plotting purposes (using something like MATLAB or Spyglass) and for quick visual inspection of the results.

Remember that the vector files elem.. are written in matrix format (with (1,1) at top-left, etc.), in contrast to their order in the input file. There are some applications in which you might want to put the results of one model run back into another model. A few MATLAB commands to do this one elem block at a time are:

- load elemds01...
- new_elem = reshape(flipud(elemds01), rows x cols, 1)

which produces a column of, in this case, elemds01 components that can be pasted directly into the input file. A more complete matlab m-file to do this for all elem-files can be found in the m-file directory.

A description of the some of the output files and the type of output is given below. Global coordinates are referenced if the quantities have x, y, or z as part of their output name (e.g., Exy, Exx, etc.) In-plane coordinates are referenced if s, d, or n (strike, dip, or normal) are part of the output file names (e.g., Esd, Ess, etc.). To avoid a cluttered look, we refer always to the global coordinates in the list below.

The underlined word or portion of a word is the name of the output file as it should be entered in the input file. For example, if you want the displacement gradient tensor, enter grad, or for rigid-body rotations, enter rbr.

Each filename, in italics below, whether output in
xyz or in vector format carries a suffix (e.g., *strainten*.<suffix>).
For clarity, this is not indicated below.

xyz: *stressten* (x y z sxx txy txz syy tyz szz tmax tmax's rake, tstrike, tdip)

Note: tmax is the max shear stress, tstrike and tdip are shear stresses along strike and dip of the inspection plane.

vector: *Sxx, Sxy, Sxz, Syy, Syz, Szz, Syy, tmax, tmaxa, tmx, tmy*

xyz: *strainten* (x y z Exx Exy Exz Eyy Eyz Ezz Evol)

Note: Evol is the trace of the tensor and is equal to the volume change.

vector: *Exx Exy Exz Eyy Eyz Ezz Evol*

Although we label these quantities as simple strains, they ARE indeed strain tensor components, e.g., *Exy* = (dUx/dy + dUy/dx)/2, etc.

xyz: *displacements* (x y z Ux Uy Uz)

vector: *disp1, disp2, disp3* (1, 2, and 3 refer to x, y, and z)

xyz: *dgten* (x y z dUx/dx, dUx/dy, dUx/dz, dUy/dx, dUy/dy, dUy/dz, dUz/dx, dUz/dy, dUz/dz)

vector: *dU/dx, dU/dy,* etc.

xyz: *strainorient* (x y z Emax, Pmax, Tmax, Eint, Pint, Tint, Emin, Pmin, Tmin

Note: For each principal strain axis, in descending order, given is the magnitude, trend, and plunge of the strain.

vector: *prince_max, plunge_max, trend_max, prince_int, plunge_int, trend_int, prince_min, plunge_min, trend_min*.

xyz: *invariants* (x y z Evol, failure stresses, octahedral shear stress, strain energy density)

vector: *delv, aftershocks, oss, work*

xyz: *rbr* (x y z rbr_x rbr_y rbr_z)

vector: *rbr_deg_x, rbr_deg_y, rbr_deg_z*

xyz: *failure* (x y z strike dip sv_east, sv_north, sv_up)

Note: sv_east, etc. are direction cosines of the slip vector w/r to east (global x-axis), etc.

vector: *strike1, strike2, ce1, cn1, cz1, ce2, cn2, cz2*

Note: *ce1*, etc. are the direction cosines of the slip vector for plane 1w/r to east (global x-axis), etc.

xyz: *elements*

vector: *elemds01, elemdd01, elemdn01, *etc.

Note: ds, dd, and dn refer to displacements of the hanging wall relative to
the footwall in the strike, dip, and normal directions, respectively.
Thus, negative values of *elemds*, for example, mean that the
displacement is dextral or right-lateral; negative values of *elemdd*
refers to a normal fault, etc. The numeral 01 refers to the first
element specified in the input file; 02 would refer to the second, and
so on.

Last modified: 2019-12-24 02:29:41 America/Denver