diff --git a/doc/Section_commands.html b/doc/Section_commands.html
index d3268a895efa5457649cc87f979ef2c83ac1dafe..5d56e8872428250cccf5b638b57ee54f803bc316 100644
--- a/doc/Section_commands.html
+++ b/doc/Section_commands.html
@@ -327,8 +327,8 @@ of each style or click on the style itself for a full description:
 <TR ALIGN="center"><TD ><A HREF = "fix_nve_asphere.html">nve/asphere</A></TD><TD ><A HREF = "fix_nve_limit.html">nve/limit</A></TD><TD ><A HREF = "fix_nve_noforce.html">nve/noforce</A></TD><TD ><A HREF = "fix_nve_sphere.html">nve/sphere</A></TD><TD ><A HREF = "fix_nvt.html">nvt</A></TD><TD ><A HREF = "fix_nvt_asphere.html">nvt/asphere</A></TD><TD ><A HREF = "fix_nvt_sllod.html">nvt/sllod</A></TD><TD ><A HREF = "fix_nvt_sphere.html">nvt/sphere</A></TD></TR>
 <TR ALIGN="center"><TD ><A HREF = "fix_orient_fcc.html">orient/fcc</A></TD><TD ><A HREF = "fix_planeforce.html">planeforce</A></TD><TD ><A HREF = "fix_poems.html">poems</A></TD><TD ><A HREF = "fix_pour.html">pour</A></TD><TD ><A HREF = "fix_press_berendsen.html">press/berendsen</A></TD><TD ><A HREF = "fix_print.html">print</A></TD><TD ><A HREF = "fix_rdf.html">rdf</A></TD><TD ><A HREF = "fix_reax_bonds.html">reax/bonds</A></TD></TR>
 <TR ALIGN="center"><TD ><A HREF = "fix_recenter.html">recenter</A></TD><TD ><A HREF = "fix_rigid.html">rigid</A></TD><TD ><A HREF = "fix_setforce.html">setforce</A></TD><TD ><A HREF = "fix_shake.html">shake</A></TD><TD ><A HREF = "fix_spring.html">spring</A></TD><TD ><A HREF = "fix_spring_rg.html">spring/rg</A></TD><TD ><A HREF = "fix_spring_self.html">spring/self</A></TD><TD ><A HREF = "fix_temp_berendsen.html">temp/berendsen</A></TD></TR>
-<TR ALIGN="center"><TD ><A HREF = "fix_temp_rescale.html">temp/rescale</A></TD><TD ><A HREF = "fix_thermal_conductivity.html">thermal/conductivity</A></TD><TD ><A HREF = "fix_tmd.html">tmd</A></TD><TD ><A HREF = "fix_ttm.html">ttm</A></TD><TD ><A HREF = "fix_viscosity.html">viscosity</A></TD><TD ><A HREF = "fix_viscous.html">viscous</A></TD><TD ><A HREF = "fix_wall_colloid.html">wall/colloid</A></TD><TD ><A HREF = "fix_wall_gran.html">wall/gran</A></TD></TR>
-<TR ALIGN="center"><TD ><A HREF = "fix_wall_lj126.html">wall/lj126</A></TD><TD ><A HREF = "fix_wall_lj93.html">wall/lj93</A></TD><TD ><A HREF = "fix_wall_reflect.html">wall/reflect</A></TD><TD ><A HREF = "fix_wiggle.html">wiggle</A> 
+<TR ALIGN="center"><TD ><A HREF = "fix_temp_rescale.html">temp/rescale</A></TD><TD ><A HREF = "fix_thermal_conductivity.html">thermal/conductivity</A></TD><TD ><A HREF = "fix_tmd.html">tmd</A></TD><TD ><A HREF = "fix_ttm.html">ttm</A></TD><TD ><A HREF = "fix_viscosity.html">viscosity</A></TD><TD ><A HREF = "fix_viscous.html">viscous</A></TD><TD ><A HREF = "fix_wall.html">wall/colloid</A></TD><TD ><A HREF = "fix_wall_gran.html">wall/gran</A></TD></TR>
+<TR ALIGN="center"><TD ><A HREF = "fix_wall.html">wall/lj126</A></TD><TD ><A HREF = "fix_wall.html">wall/lj93</A></TD><TD ><A HREF = "fix_wall_reflect.html">wall/reflect</A></TD><TD ><A HREF = "fix_wiggle.html">wiggle</A> 
 </TD></TR></TABLE></DIV>
 
 <P>These are fix styles contributed by users, which can be used if
diff --git a/doc/fix_wall.html b/doc/fix_wall.html
index 7b05bd158c4fc1ca161de82c52b6615fca9abdc1..c0638b4e6dfaf0f5e0d51eea5a7efeb40e4c1da8 100644
--- a/doc/fix_wall.html
+++ b/doc/fix_wall.html
@@ -50,33 +50,34 @@ fix zwalls all wall/colloid zlo 0.0 1.0 1.0 0.858 zhi 40.0 1.0 1.0 0.858
 </P>
 <P>Bound the simulation domain on one or more of its faces with a flat
 wall that interacts with the atoms in the group by generating a force
-on the atom in a direction perpendicular to the wall.  The energy E of
+on the atom in a direction perpendicular to the wall.  The energy of
 wall-particle interactions depends on the style.
 </P>
-<P>For style <I>wall/lj93</I>, the energy is given by the 9/3 potential:
+<P>For style <I>wall/lj93</I>, the energy E is given by the 9/3 potential:
 </P>
 <CENTER><IMG SRC = "Eqs/fix_wall_lj93.jpg">
 </CENTER>
-<P>For style <I>wall/lj126</I>, the energy is given by the 12/6 potential:
+<P>For style <I>wall/lj126</I>, the energy E is given by the 12/6 potential:
 </P>
 <CENTER><IMG SRC = "Eqs/pair_lj.jpg">
 </CENTER>
-<P>For style <I>wall/colloid</I>, the energy is given by an integrated form of
+<P>For style <I>wall/colloid</I>, the energy E is given by an integrated form of
 the <A HREF = "pair_colloid.html">pair_style colloid</A> potential:
 </P>
 <CENTER><IMG SRC = "Eqs/fix_wall_colloid.jpg">
 </CENTER>
-<P>In all cases, <I>r</I> is the distance from the particle to the wall
-<I>coord</I>, and Rc is the <I>cutoff</I> distance at which the particle and
-wall no longer interact, as specified in the command.  The energy of
-the wall potential is shifted so that the energy of the wall-particle
-interaction is 0.0 at the cutoff distance.
+<P>In all cases, <I>r</I> is the distance from the particle to the wall at
+position <I>coord</I>, and Rc is the <I>cutoff</I> distance at which the
+particle and wall no longer interact.  The energy of the wall
+potential is shifted so that the wall-particle interaction energy is
+0.0 at the cutoff distance.
 </P>
 <P>For the <I>wall/lj93</I> and <I>wall/lj126</I> styles, <I>epsilon</I> and <I>sigma</I> are
 the usual Lennard-Jones parameters, which determine the strength and
 size of the particle as it interacts with the wall.  Note that this
-<I>sigma</I> may be different than any <I>sigma</I> values defined for a pair
-style that computes particle-particle interactions.
+<I>epsilon</I> and <I>sigma</I> may be different than any <I>epsilon</I> or <I>sigma</I>
+values defined for a pair style that computes particle-particle
+interactions.
 </P>
 <P>The <I>wall/lj93</I> interaction is derived by integrating over a 3d
 half-lattice of Lennard-Jones 12/6 particles.  The <I>wall/lj126</I>
@@ -92,18 +93,18 @@ center to the wall.
 </P>
 <P>The <I>wall/colloid</I> interaction is derived by integrating over
 constituent LJ particles of size <I>sigma</I> within the colloid particle
-and a 3d half-lattice of Lennard-Jones 12/6 particles with the same
-<I>sigma</I> in the wall.
+and a 3d half-lattice of Lennard-Jones 12/6 particles of size <I>sigma</I>
+in the wall.
 </P>
 <P>If the <I>vel</I> keyword is specified, the position of all walls will move
 during the simulation in a perpendicular direction, based on their
-initial position (coord), the specified velocity (vel), and the time
+initial <I>coord</I> position, the specified velocity <I>vel</I>, and the time
 elapsed since the beginning of the simulation.  See the note below
 about making a wall move continuously across multiple runs.  A
 positive velocity means each wall moves inward, towards the center of
-the box.  I.e. an <I>xlo</I> wall will move in the +x direction; an <I>xhi</I>
-wall will move in the -x direction.  A negative velocity means each
-wall moves outward, away from the center of the box.  If you want
+the box.  I.e. an <I>xlo</I> wall will move in the +x direction and an
+<I>xhi</I> wall will move in the -x direction.  A negative velocity means
+each wall moves outward, away from the center of the box.  If you want
 different walls to move with different velocities, then you need to
 use multiple fix wall commands.
 </P>
@@ -113,13 +114,13 @@ direction, according to the equation:
 </P>
 <PRE>position = pos0 + A sin (omega * delta) 
 </PRE>
-<P>where <I>pos0</I> is the position at the time the beginning of the
-simulation, <I>A</I> is the <I>amplitude</I>, <I>omega</I> is 2 PI / <I>period</I>, and
-<I>delta</I> is the time elapsed since the beginning of the simulation.
-See the note below about making a wall oscillate continuously across
-multiple runs.  A positive amplitude means each wall initially moves
-inward, towards the center of the box.  I.e. an <I>xlo</I> wall will move
-initially in the +x direction; an <I>xhi</I> wall will move in the -x
+<P>where <I>pos0</I> is the position at the beginning of the simulation, <I>A</I>
+is the <I>amplitude</I>, <I>omega</I> is 2 PI / <I>period</I>, and <I>delta</I> is the
+time elapsed since the beginning of the simulation.  See the note
+below about making a wall oscillate continuously across multiple runs.
+A positive amplitude means each wall initially moves inward, towards
+the center of the box.  I.e. an <I>xlo</I> wall will move initially in the
++x direction and an <I>xhi</I> wall will move initially in the -x
 direction.  A negative velocity means each wall moves initially
 outward, away from the center of the box.  If you want different walls
 to oscillate with different amplitudes or periods, then you need to
@@ -135,14 +136,18 @@ fix to add the energy of interaction between atoms and each wall to
 the system's potential energy as part of <A HREF = "thermo_style.html">thermodynamic
 output</A>.
 </P>
-<P>This fix computes a scalar energy and a 3-vector of forces (on the
-walls), which can be accessed by various <A HREF = "Section_howto.html#4_15">output
+<P>This fix computes a scalar energy and a 6-length vector of forces (one
+force magnitude per wall), which can be accessed by various <A HREF = "Section_howto.html#4_15">output
 commands</A>.  The scalar and vector values
 calculated by this fix are "extensive", meaning they scale with the
-number of atoms in the simulation.  Note that if you define multiple
-walls, then the energy and force will be summed over all the walls.
-If you want the energy/force on a per-wall basis, you need to use
-multiple fix wall commands.
+number of atoms in the simulation.  Note that the scalar energy is the
+sum of interactions with all defined walls.  If you want the energy on
+a per-wall basis, you need to use multiple fix wall commands.  The 6
+vector quantities are the force on the <I>xlo</I> wall, the <I>xhi</I> wall,
+<I>ylo</I>, <I>yhi</I>, <I>zlo</I>, <I>zhi</I>.  These values will only be non-zero if the
+corresponding wall is defined.  Note that an outward force on a wall
+will be a negative value for <I>lo</I> walls and a positive value for <I>hi</I>
+walls.
 </P>
 <P>This fix can change the position of each wall, due to the <I>vel</I> or
 <I>wiggle</I> keywords, continuously over multiple runs, using the <I>start</I>
diff --git a/doc/fix_wall.txt b/doc/fix_wall.txt
index 1c4728a6d60e7d3fbfba75d88224b1ddda7f9206..5894ba82565ce2a59611e35009979c219e485980 100644
--- a/doc/fix_wall.txt
+++ b/doc/fix_wall.txt
@@ -126,14 +126,18 @@ fix to add the energy of interaction between atoms and each wall to
 the system's potential energy as part of "thermodynamic
 output"_thermo_style.html.
 
-This fix computes a scalar energy and a 3-vector of forces (on the
-walls), which can be accessed by various "output
+This fix computes a scalar energy and a 6-length vector of forces (one
+force magnitude per wall), which can be accessed by various "output
 commands"_Section_howto.html#4_15.  The scalar and vector values
 calculated by this fix are "extensive", meaning they scale with the
-number of atoms in the simulation.  Note that if you define multiple
-walls, then the energy and force will be summed over all the walls.
-If you want the energy/force on a per-wall basis, you need to use
-multiple fix wall commands.
+number of atoms in the simulation.  Note that the scalar energy is the
+sum of interactions with all defined walls.  If you want the energy on
+a per-wall basis, you need to use multiple fix wall commands.  The 6
+vector quantities are the force on the {xlo} wall, the {xhi} wall,
+{ylo}, {yhi}, {zlo}, {zhi}.  These values will only be non-zero if the
+corresponding wall is defined.  Note that an outward force on a wall
+will be a negative value for {lo} walls and a positive value for {hi}
+walls.
 
 This fix can change the position of each wall, due to the {vel} or
 {wiggle} keywords, continuously over multiple runs, using the {start}