MS_load

MS_LOAD - Load a set of elastic constants from a file.
// Part of MSAT - The Matlab Seismic Anisotropy Toolkit //
[ ... ] = MS_load( fname, ... )
Usage:
    [C] = MS_load(fname)
        Load the elastic constants in file 'fname'. See file format below.
[C,r] = MS_load(fname,...)
    Load the elastic constants and density in file 'fname'. See file
    format below.
[C,...] = MS_load(fname,...,'format',fmt)
    Specify format for file. See below for descriptions. Available options
    are:
       'simple' : Formatted text file (default)
[C,...] = MS_load(fname,...,'symmetry', mode)
    ** some parts not yet implemented **
    Fill out elastic tensor based on symmetry, defined by mode. This can
    take the following values:
       'none' - nothing attempted, unspecified Cijs are zero (default)
       'auto' - assume symmetry based on number of Cijs specified
       'iso' - isotropic (nCij=2) ; C33 and C66 must be specified.
       'hex' - hexagonal (nCij=5) ; C11, C33, C44, C66 and C13 must be
                  specified, x3 is symmetry axis
       'vti' - synonym for hexagonal
       'cubic' - cubic (nCij=3) ; C33, C66 and C12 must be specified
[C,...] = MS_load(fname,...,'eunit',unit)
    Specify elasticity unit; this is used to convert to GPa (msat's
    default operating unit)
       'GPa' - elastic constants are already in GPa (default)
       'Mbar' - elastic constants are in Mbar
       'Pa' - elastic constants are in Pascals
       'bar' - elastic constants are in bar
[C,...] = MS_load(fname,...,'dunit',unit)
    Specify density unit; this is used to convert to kg/m3 (msat's
    default operating unit)
       'kgm3' - density is already in kg/m3 (default)
       'gcc' - density is in g/cc
[C,...] = MS_load(fname,...,'Aij')
    Elasticities are specified in density normalised form. This option
    requires a density to be specified, the elasticities are returned with
    the normalisation reversed (default format for msat). The density
    normalisation is assumed to have been performed using the density
    specified in the file (so, obviously, one is required) without any
    different unit scalings. It is *highly* advisable to check (for
    example using MS_phasevels) that sensible numbers result. The default
    is to assume that the input is not density normalised.
[C,...] = MS_load(fname,...,'force')
    Disable post-load checking. You're on your own, buddy.
Supported file formats:
-----------------------
  • simple (default)

    In this mode MS_load expects an ascii file containing lines of the form:
    I J C(I,J)
    Elasticity tensor components are specified by their indices (i,j). Any
    number other than 1-6 in the i or j column denotes that a density is
    being given.
    So, for olivine, the file might look like:
    1 1 320.5
    1 2 68.1
    1 3 71.6
    ...
    6 6 78.7
    0 0 3325
    Text following a '' in any line is ignored (useful for header lines or
    comments)
  • ematrix

    Elastic constants in the format read by EMATRIX (D. Mainprice).
    See also: MS_LOAD_LIST, MS_ELASTICDB