\(\renewcommand{\AA}{\text{Å}}\)

3.17. Thermodynamic output options

The Thermo class computes and prints thermodynamic information to the screen and log file; see the files thermo.cpp and thermo.h.

There are four styles defined in thermo.cpp: “one”, “multi”, “yaml”, and “custom”. The “custom” style allows the user to explicitly list keywords for individual quantities to print when thermodynamic output is generated. The others have a fixed list of keywords. See the thermo_style command for a list of available quantities. The formatting of the “custom” style defaults to the “one” style, but can be adapted using thermo_modify line.

The thermo styles (one, multi, etc) are defined by lists of keywords with associated formats for integer and floating point numbers and identified by an enumerator constant. Adding a new style thus mostly requires defining a new list of keywords and the associated formats and then inserting the required output processing where the enumerators are identified. Search for the word “CUSTOMIZATION” with references to “thermo style” in the thermo.cpp file to see the locations where code will need to be added. The member function Thermo::header() prints output at the very beginning of a thermodynamic output block and can be used to print column headers or other front matter. The member function Thermo::footer() prints output at the end of a thermodynamic output block. The formatting of the output is done by assembling a “line” (which may span multiple lines if the style inserts newline characters (”n” as in the “multi” style).

New thermodynamic keywords can also be added to thermo.cpp to compute new quantities for output. Search for the word “CUSTOMIZATION” with references to “keyword” in thermo.cpp to see the several locations where code will need to be added. Effectively, you need to define a member function that computes the property, add an if statement in Thermo::parse_fields() where the corresponding header string for the keyword and the function pointer is registered by calling the Thermo::addfield() method, and add an if statement in Thermo::evaluate_keyword() which is called from the Variable class when a thermo keyword is encountered.

Note

The third argument to Thermo::addfield() is a flag indicating whether the function for the keyword computes a floating point (FLOAT), regular integer (INT), or big integer (BIGINT) value. This information is used for formatting the thermodynamic output. Inside the function the result must then be stored either in the dvalue, ivalue or bivalue member variable, respectively.

Since the thermo_style custom command allows to use output of quantities calculated by fixes, computes, and variables, it may often be simpler to compute what you wish via one of those constructs, rather than by adding a new keyword to the thermo_style command.