ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@13104 f3b2605a-c512-4ea7-a41b-209d697bcdaa

Browse Source
This commit is contained in:
sjplimp
2015-02-13 20:28:48 +00:00
parent b8541b66ba
commit f28b0794ce
20 changed files with 103 additions and 109 deletions
Download Patch File Download Diff File Expand all files Collapse all files

66
doc/Section_howto.html
View File

@ -2172,7 +2172,7 @@ commands, to calculate various properties of a system:
overviewed. Then some examples are given of how to compute different overviewed. Then some examples are given of how to compute different
properties with chunk commands. properties with chunk commands.
</P> </P>
<H5><A HREF = "compute_chunk_atom.html">Compute chunk/atom</A> command <H5><A HREF = "compute_chunk_atom.html">Compute chunk/atom</A> command:
</H5> </H5>
<P>This compute can assign atoms to chunks of various styles. Only atoms <P>This compute can assign atoms to chunks of various styles. Only atoms
in the specified group and optional specified region are assigned to a in the specified group and optional specified region are assigned to a
@ -2213,7 +2213,7 @@ do this. You can also define a <A HREF = "variable.html">per-atom variable</A>
the input script that uses a formula to generate a chunk ID for each the input script that uses a formula to generate a chunk ID for each
atom. atom.
</P> </P>
<H5><A HREF = "fix_ave_chunk_atom.html">Fix ave/chunk</A> command <H5><A HREF = "fix_ave_chunk_atom.html">Fix ave/chunk</A> command:
</H5> </H5>
<P>This fix takes the ID of a <A HREF = "compute_chunk_atom.html">compute <P>This fix takes the ID of a <A HREF = "compute_chunk_atom.html">compute
chunk/atom</A> command as input. For each chunk, chunk/atom</A> command as input. For each chunk,
@ -2230,7 +2230,7 @@ values to be time-averaged in various ways, and output to a file. The
fix produces a global array as output with one row of values per fix produces a global array as output with one row of values per
chunk. chunk.
</P> </P>
<H5>Compute */chunk commands <H5>Compute */chunk commands:
</H5> </H5>
<P>Currently the following computes operate on chunks of atoms to produce <P>Currently the following computes operate on chunks of atoms to produce
per-chunk values. per-chunk values.
@ -2279,42 +2279,40 @@ largest cluster or fastest diffusing molecule.
<P>Here are eaxmples using chunk commands to calculate various <P>Here are eaxmples using chunk commands to calculate various
properties: properties:
</P> </P>
<P>(1) Mimic the deprecated fix ave/spatial command, to average atom <P>(1) Average velocity in each of 1000 2d spatial bins:
velocity in each spatial bin:
</P> </P>
<P>Old command: <PRE>compute cc1 all chunk/atom bin/2d x 0.0 0.1 y lower 0.01 units reduced
fix 1 all ave/chunk 100 10 1000 cc1 vx vy file tmp.out
</PRE>
<P>(2) Temperature in each spatial bin, after subtracting a flow
velocity:
</P> </P>
<P>fix ave/spatial <PRE>compute cc1 all chunk/atom bin/2d x 0.0 0.1 y lower 0.1 units reduced
compute vbias all temp/profile 1 0 0 y 10
fix 1 all ave/chunk 100 10 1000 cc1 temp bias vbias file tmp.out
</PRE>
<P>(3) Center of mass of each molecule:
</P> </P>
<P>New commands: <PRE>compute cc1 all chunk/atom molecule
compute myChunk all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE>
<P>(4) Total force on each molecule and ave/max across all molecules:
</P> </P>
<P>compute chunk/atom <PRE>compute cc1 all chunk/atom molecule
fix ave/chunk fix 1 all ave/chunk 1000 1 1000 cc1 fx fy fz file tmp.out
</P> variable xave equal ave(f_1<B>2</B>)
<P>(2) Mimincing the deprecated compute msd/molecule command variable xmax equal max(f_1<B>2</B>)
to compute the mean-squared displacement of each molecule: thermo 1000
</P> thermo_style custom step temp v_xave v_xmax
<P>Old commands: </PRE>
</P> <P>(5) Histogram of cluster sizes:
<P>compute molecule/msd
fix ave/time file msd.molecule
</P>
<P>New commands:
</P>
<P>compute chunk/atom style
compute molecule/chunk
fix ave/time file msd.molecule
</P>
<P>(3) print or time ave total force per molecule
average across molecules (variable special functions)
</P>
<P>(4) histogram of cluster sizes (use fix ave/histo)
</P>
<P>(5) count of # of atoms with each coord #
don't need chunking, just fix ave/histo on coord/atom
</P>
<P>(6) ave temperature per bin
</P> </P>
<PRE>compute cluster all cluster/atom 1.0
compute cc1 all chunk/atom c_cluster compress yes
compute size all property/chunk cc1 count
fix 1 all ave/histo 100 1 100 0 20 20 c_size mode vector ave running beyond ignore file tmp.histo
</PRE>
<HR> <HR>
<HR> <HR>

56
doc/Section_howto.txt
View File

@ -2159,7 +2159,7 @@ Here, each of the 3 kinds of chunk-related commands is briefly
overviewed. Then some examples are given of how to compute different overviewed. Then some examples are given of how to compute different
properties with chunk commands. properties with chunk commands.
"Compute chunk/atom"_compute_chunk_atom.html command :h5 "Compute chunk/atom"_compute_chunk_atom.html command: :h5
This compute can assign atoms to chunks of various styles. Only atoms This compute can assign atoms to chunks of various styles. Only atoms
in the specified group and optional specified region are assigned to a in the specified group and optional specified region are assigned to a
@ -2198,7 +2198,7 @@ do this. You can also define a "per-atom variable"_variable.html in
the input script that uses a formula to generate a chunk ID for each the input script that uses a formula to generate a chunk ID for each
atom. atom.
"Fix ave/chunk"_fix_ave_chunk_atom.html command :h5 "Fix ave/chunk"_fix_ave_chunk_atom.html command: :h5
This fix takes the ID of a "compute This fix takes the ID of a "compute
chunk/atom"_compute_chunk_atom.html command as input. For each chunk, chunk/atom"_compute_chunk_atom.html command as input. For each chunk,
@ -2215,7 +2215,7 @@ values to be time-averaged in various ways, and output to a file. The
fix produces a global array as output with one row of values per fix produces a global array as output with one row of values per
chunk. chunk.
Compute */chunk commands :h5 Compute */chunk commands: :h5
Currently the following computes operate on chunks of atoms to produce Currently the following computes operate on chunks of atoms to produce
per-chunk values. per-chunk values.
@ -2264,41 +2264,39 @@ Example calculations with chunks :h5
Here are eaxmples using chunk commands to calculate various Here are eaxmples using chunk commands to calculate various
properties: properties:
(1) Mimic the deprecated fix ave/spatial command, to average atom (1) Average velocity in each of 1000 2d spatial bins:
velocity in each spatial bin:
Old command: compute cc1 all chunk/atom bin/2d x 0.0 0.1 y lower 0.01 units reduced
fix 1 all ave/chunk 100 10 1000 cc1 vx vy file tmp.out :pre
fix ave/spatial (2) Temperature in each spatial bin, after subtracting a flow
velocity:
New commands: compute cc1 all chunk/atom bin/2d x 0.0 0.1 y lower 0.1 units reduced
compute vbias all temp/profile 1 0 0 y 10
fix 1 all ave/chunk 100 10 1000 cc1 temp bias vbias file tmp.out :pre
compute chunk/atom (3) Center of mass of each molecule:
fix ave/chunk
(2) Mimincing the deprecated compute msd/molecule command compute cc1 all chunk/atom molecule
to compute the mean-squared displacement of each molecule: compute myChunk all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre
Old commands: (4) Total force on each molecule and ave/max across all molecules:
compute molecule/msd compute cc1 all chunk/atom molecule
fix ave/time file msd.molecule fix 1 all ave/chunk 1000 1 1000 cc1 fx fy fz file tmp.out
variable xave equal ave(f_1[2])
variable xmax equal max(f_1[2])
thermo 1000
thermo_style custom step temp v_xave v_xmax :pre
New commands: (5) Histogram of cluster sizes:
compute chunk/atom style compute cluster all cluster/atom 1.0
compute molecule/chunk compute cc1 all chunk/atom c_cluster compress yes
fix ave/time file msd.molecule compute size all property/chunk cc1 count
fix 1 all ave/histo 100 1 100 0 20 20 c_size mode vector ave running beyond ignore file tmp.histo :pre
(3) print or time ave total force per molecule
average across molecules (variable special functions)
(4) histogram of cluster sizes (use fix ave/histo)
(5) count of # of atoms with each coord #
don't need chunking, just fix ave/histo on coord/atom
(6) ave temperature per bin
:line :line
:line :line

2
doc/compute_com_chunk.html
View File

@ -63,7 +63,7 @@ image</A> command.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all com/chunk cc1 compute myChunk all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_com_chunk.txt
View File

@ -60,7 +60,7 @@ The simplest way to output the results of the compute com/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all com/chunk cc1 compute myChunk all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_gyration_chunk.html
View File

@ -90,7 +90,7 @@ image</A> command.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all gyration/chunk cc1 compute myChunk all gyration/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_gyration_chunk.txt
View File

@ -81,7 +81,7 @@ The simplest way to output the results of the compute gyration/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all gyration/chunk cc1 compute myChunk all gyration/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_inertia_chunk.html
View File

@ -64,7 +64,7 @@ image</A> command.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all inertia/chunk cc1 compute myChunk all inertia/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_inertia_chunk.txt
View File

@ -61,7 +61,7 @@ The simplest way to output the results of the compute inertia/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all inertia/chunk cc1 compute myChunk all inertia/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_msd_chunk.html
View File

@ -97,7 +97,7 @@ file</A>.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all com/msd cc1 compute myChunk all com/msd cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_msd_chunk.txt
View File

@ -94,7 +94,7 @@ The simplest way to output the results of the compute com/msd
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all com/msd cc1 compute myChunk all com/msd cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_property_chunk.html
View File

@ -92,7 +92,7 @@ will be in unitless reduced units (0-1).
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk1 all property/chunk cc1 compute myChunk1 all property/chunk cc1
compute myChunk2 all com/chunk cc1 compute myChunk2 all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk1 c_myChunk2 file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk1 c_myChunk2 file tmp.out mode vector

2
doc/compute_property_chunk.txt
View File

@ -85,7 +85,7 @@ The simplest way to output the results of the compute property/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk1 all property/chunk cc1 compute myChunk1 all property/chunk cc1
compute myChunk2 all com/chunk cc1 compute myChunk2 all com/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk1 c_myChunk2 file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk1 c_myChunk2 file tmp.out mode vector :pre

2
doc/compute_temp_chunk.html
View File

@ -114,7 +114,7 @@ non-zero chunk IDs.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all temp/chunk cc1 compute myChunk all temp/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_temp_chunk.txt
View File

@ -104,7 +104,7 @@ The simplest way to output the results of the compute temp/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all temp/chunk cc1 compute myChunk all temp/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_torque_chunk.html
View File

@ -64,7 +64,7 @@ image</A> command.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all torque/chunk cc1 compute myChunk all torque/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_torque_chunk.txt
View File

@ -61,7 +61,7 @@ The simplest way to output the results of the compute torque/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all torque/chunk cc1 compute myChunk all torque/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

2
doc/compute_vcm_chunk.html
View File

@ -55,7 +55,7 @@ non-zero chunk IDs.
calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A> calculation to a file is to use the <A HREF = "fix_ave_time.html">fix ave/time</A>
command, for example: command, for example:
</P> </P>
<PRE>compute cc1 chunk/atom molecule <PRE>compute cc1 all chunk/atom molecule
compute myChunk all vcm/chunk cc1 compute myChunk all vcm/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector
</PRE> </PRE>

2
doc/compute_vcm_chunk.txt
View File

@ -52,7 +52,7 @@ The simplest way to output the results of the compute vcm/chunk
calculation to a file is to use the "fix ave/time"_fix_ave_time.html calculation to a file is to use the "fix ave/time"_fix_ave_time.html
command, for example: command, for example:
compute cc1 chunk/atom molecule compute cc1 all chunk/atom molecule
compute myChunk all vcm/chunk cc1 compute myChunk all vcm/chunk cc1
fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre fix 1 all ave/time 100 1 100 c_myChunk file tmp.out mode vector :pre

29
doc/fix_ave_chunk.html
View File

@ -72,7 +72,7 @@
</UL> </UL>
<P><B>Examples:</B> <P><B>Examples:</B>
</P> </P>
<PRE>fix 1 all ave/chunk 10000 1 10000 binchunk c_myCentro ttle1 "My output values" <PRE>fix 1 all ave/chunk 10000 1 10000 binchunk c_myCentro title1 "My output values"
fix 1 flow ave/chunk 100 10 1000 molchunk vx vz norm sample file vel.profile fix 1 flow ave/chunk 100 10 1000 molchunk vx vz norm sample file vel.profile
fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running
fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running
@ -95,11 +95,11 @@ fix 1 flow ave/chunk 100 10 1000 cc1 vx vz norm sample file vel.profile
<P><B>Description:</B> <P><B>Description:</B>
</P> </P>
<P>Use one or more per-atom vectors as inputs every few timesteps, sum <P>Use one or more per-atom vectors as inputs every few timesteps, sum
their values for multiple chunks of atoms, and average the values for the values over the atoms in each chunk at each timestep, then average
each chunk over longer timescales. The resulting chunk averages can the per-chunk values over longer timescales. The resulting chunk
be used by other <A HREF = "Section_howto.html#howto_15">output commands</A> such as averages can be used by other <A HREF = "Section_howto.html#howto_15">output
<A HREF = "thermo_style.html">thermo_style custom</A>, and can also be written to a commands</A> such as <A HREF = "thermo_style.html">thermo_style
file. custom</A>, and can also be written to a file.
</P> </P>
<P>In LAMMPS, chunks are collections of atoms defined by a <A HREF = "compute_chunk_atom.html">compute <P>In LAMMPS, chunks are collections of atoms defined by a <A HREF = "compute_chunk_atom.html">compute
chunk/atom</A> command, which assigns each atom chunk/atom</A> command, which assigns each atom
@ -195,15 +195,14 @@ property/atom</A> command and then specifying
an input value from that compute. an input value from that compute.
</P> </P>
<P>The <I>density/number</I> value means the number density is computed for <P>The <I>density/number</I> value means the number density is computed for
each chunk, i.e. a weighting of 1 for each atom. The <I>density/mass</I> each chunk, i.e. number/volume. The <I>density/mass</I> value means the
value means the mass density is computed for each , i.e. each atom is mass density is computed for each chunk, i.e. total-mass/volume. The
weighted by its mass. The resulting density is normalized by the output values are in units of 1/volume or density (mass/volume). See
volume of the chunk so that units of number/volume or density are the <A HREF = "units.html">units</A> command doc page for the definition of density
output. See the <A HREF = "units.html">units</A> command doc page for the for each choice of units, e.g. gram/cm^3. If the chunks defined by
definition of density for each choice of units, e.g. gram/cm^3. If the <A HREF = "compute_chunk_atom.html">compute chunk/atom</A> command are spatial
the chunks defined by the <A HREF = "compute_chunk_atom.html">compute chunk/atom</A> bins, the volume is the bin volume. Otherwise it is the volume of the
command are spatial bins, the volume is the bin volume. Otherwise it entire simulation box.
is the volume of the entire simulation box.
</P> </P>
<P>The <I>temp</I> value means the temperature is computed for each chunk, by <P>The <I>temp</I> value means the temperature is computed for each chunk, by
the formula KE = DOF/2 k T, where KE = total kinetic energy of the the formula KE = DOF/2 k T, where KE = total kinetic energy of the

29
doc/fix_ave_chunk.txt
View File

@ -58,7 +58,7 @@ keyword = {norm} or {ave} or {bias} or {adof} or {cdof} or {file} or {overwrite}
[Examples:] [Examples:]
fix 1 all ave/chunk 10000 1 10000 binchunk c_myCentro ttle1 "My output values" fix 1 all ave/chunk 10000 1 10000 binchunk c_myCentro title1 "My output values"
fix 1 flow ave/chunk 100 10 1000 molchunk vx vz norm sample file vel.profile fix 1 flow ave/chunk 100 10 1000 molchunk vx vz norm sample file vel.profile
fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running
fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running :pre fix 1 flow ave/chunk 100 5 1000 binchunk density/mass ave running :pre
@ -81,11 +81,11 @@ fix 1 flow ave/chunk 100 10 1000 cc1 vx vz norm sample file vel.profile :pre
[Description:] [Description:]
Use one or more per-atom vectors as inputs every few timesteps, sum Use one or more per-atom vectors as inputs every few timesteps, sum
their values for multiple chunks of atoms, and average the values for the values over the atoms in each chunk at each timestep, then average
each chunk over longer timescales. The resulting chunk averages can the per-chunk values over longer timescales. The resulting chunk
be used by other "output commands"_Section_howto.html#howto_15 such as averages can be used by other "output
"thermo_style custom"_thermo_style.html, and can also be written to a commands"_Section_howto.html#howto_15 such as "thermo_style
file. custom"_thermo_style.html, and can also be written to a file.
In LAMMPS, chunks are collections of atoms defined by a "compute In LAMMPS, chunks are collections of atoms defined by a "compute
chunk/atom"_compute_chunk_atom.html command, which assigns each atom chunk/atom"_compute_chunk_atom.html command, which assigns each atom
@ -181,15 +181,14 @@ property/atom"_compute_property_atom.html command and then specifying
an input value from that compute. an input value from that compute.
The {density/number} value means the number density is computed for The {density/number} value means the number density is computed for
each chunk, i.e. a weighting of 1 for each atom. The {density/mass} each chunk, i.e. number/volume. The {density/mass} value means the
value means the mass density is computed for each , i.e. each atom is mass density is computed for each chunk, i.e. total-mass/volume. The
weighted by its mass. The resulting density is normalized by the output values are in units of 1/volume or density (mass/volume). See
volume of the chunk so that units of number/volume or density are the "units"_units.html command doc page for the definition of density
output. See the "units"_units.html command doc page for the for each choice of units, e.g. gram/cm^3. If the chunks defined by
definition of density for each choice of units, e.g. gram/cm^3. If the "compute chunk/atom"_compute_chunk_atom.html command are spatial
the chunks defined by the "compute chunk/atom"_compute_chunk_atom.html bins, the volume is the bin volume. Otherwise it is the volume of the
command are spatial bins, the volume is the bin volume. Otherwise it entire simulation box.
is the volume of the entire simulation box.
The {temp} value means the temperature is computed for each chunk, by The {temp} value means the temperature is computed for each chunk, by
the formula KE = DOF/2 k T, where KE = total kinetic energy of the the formula KE = DOF/2 k T, where KE = total kinetic energy of the