MS_LOAD - Load a set of elastic constants from a file.

[ ... ] = 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)
       'ematrix' : Format read by EMATRIX (D. Mainprice).

[C,...] = MS_load(fname,...,'symmetry', mode)
    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
    Note: when using a density normalised tensor, specify the unit
    *before* normalisation (i.e., a unit of pressure).

[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. That is, the values of the elastic
    constants are read from the file and multiplied by the value of the
    density read from the file before any unit conversion is performed.
    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:
-----------------------