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