MS_save

MS_SAVE - Write a set of elastic constants from a file.
MS_save( fname, C, rho, ... )
Usage:
    MS_save(fname, C, rho)
         Write the elastic constants, C, and density, rho, to file
         'fname'. C must be provided in GPa and rho must be provided in
         kg/m3. See file format below.
MS_save(fname, C, rho, ..., '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).
MS_save(fname, C, rho, ..., 'eunit',unit)
    Specify elasticity unit; this is used to convert from GPa (msat's
    default operating unit)
       'GPa' - elastic constants are written in GPa (default)
       'Mbar' - elastic constants are written in Mbar
       'Pa' - elastic constants are written in Pascals
       'bar' - elastic constants are written in bar
    Note: when using a density normalised tensor, specify the unit
    *before* normalisation (i.e., a unit of pressure).
MS_save(fname, C, rho, ..., 'dunit',unit)
    Specify density unit; this is used to convert from kg/m3 (msat's
    default operating unit)
       'kgm3' - density is written in kg/m3 (default)
       'gcc' - density is written g/cc
MS_save(fname, C, rho, ..., 'Aij')
    Elasticities are written in density normalised form. This option
    requires a density to be specified. The provided elastic
    constants are first converted to the unit specified by the
    'eunit' argument (if provided) and the density converted to the
    unit specified in the by the 'dunit' argument (if provided) then
    the elastic constants are divided by the density before being
    written to the file.
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