From fa66305dead7ecc875f8fcf568b0c6342c3837cc Mon Sep 17 00:00:00 2001
From: sjplimp <sjplimp@f3b2605a-c512-4ea7-a41b-209d697bcdaa>
Date: Wed, 21 Apr 2010 20:19:43 +0000
Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@4017
 f3b2605a-c512-4ea7-a41b-209d697bcdaa

---
 doc/fix_heat.html | 33 ++++++++++++++++++++++++---------
 doc/fix_heat.txt  | 28 ++++++++++++++++++----------
 2 files changed, 42 insertions(+), 19 deletions(-)

diff --git a/doc/fix_heat.html b/doc/fix_heat.html
index 53e15feab1..5daf47d5cd 100644
--- a/doc/fix_heat.html
+++ b/doc/fix_heat.html
@@ -15,26 +15,41 @@
 </P>
 <PRE>fix ID group-ID heat N eflux 
 </PRE>
-<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command
-<LI>heat = style name of this fix command
-<LI>N = add/subtract heat every this many timesteps
+<UL><LI>ID, group-ID are documented in <A HREF = "fix.html">fix</A> command 
+
+<LI>heat = style name of this fix command 
+
+<LI>N = add/subtract heat every this many timesteps 
+
 <LI>eflux = rate of heat addition or subtraction (energy/time units) 
+
+<LI>zero or more keyword/value pairs may be appended to args 
+
+<LI>keyword = <I>region</I>  
+
+<PRE>  <I>region</I> value = region-ID
+    region-ID = ID of region atoms must be in to have added force 
+</PRE>
+
 </UL>
 <P><B>Examples:</B>
 </P>
 <PRE>fix 3 qin heat 1 1.0
-fix 4 qout heat 1 -1.0 
+fix 4 qout heat 1 -1.0 region top 
 </PRE>
 <P><B>Description:</B>
 </P>
 <P>Add non-translational kinetic energy (heat) to a group of atoms such
 that their aggregate momentum is conserved.  Two of these fixes can be
 used to establish a temperature gradient across a simulation domain by
-adding heat to one group of atoms (hot reservoir) and subtracting heat
-from another (cold reservoir).  E.g. a simulation sampling from the
-McDLT ensemble.  Note that the fix is applied to a group of atoms, not
-a geometric region, so that the same set of atoms is affected wherever
-they may move to.
+adding heat (energy) to one group of atoms (hot reservoir) and
+subtracting heat from another (cold reservoir).  E.g. a simulation
+sampling from the McDLT ensemble.
+</P>
+<P>If the <I>region</I> keyword is used, the atom must be in both the group
+and the specified geometric <A HREF = "region.html">region</A> in order to have
+energy added or subtracted to it.  If not specified, then the atoms in
+the group are affected wherever they may move to.
 </P>
 <P>Heat addition/subtraction is performed every N timesteps.  The <I>eflux</I>
 parameter determines the change in aggregate energy of the entire
diff --git a/doc/fix_heat.txt b/doc/fix_heat.txt
index 2f9e937892..2a494d5e66 100644
--- a/doc/fix_heat.txt
+++ b/doc/fix_heat.txt
@@ -12,26 +12,34 @@ fix heat command :h3
 
 fix ID group-ID heat N eflux :pre
 
-ID, group-ID are documented in "fix"_fix.html command
-heat = style name of this fix command
-N = add/subtract heat every this many timesteps
-eflux = rate of heat addition or subtraction (energy/time units) :ul
+ID, group-ID are documented in "fix"_fix.html command :ulb,l
+heat = style name of this fix command :l
+N = add/subtract heat every this many timesteps :l
+eflux = rate of heat addition or subtraction (energy/time units) :l
+zero or more keyword/value pairs may be appended to args :l
+keyword = {region}  :l
+  {region} value = region-ID
+    region-ID = ID of region atoms must be in to have added force :pre
+:ule
 
 [Examples:]
 
 fix 3 qin heat 1 1.0
-fix 4 qout heat 1 -1.0 :pre
+fix 4 qout heat 1 -1.0 region top :pre
 
 [Description:]
 
 Add non-translational kinetic energy (heat) to a group of atoms such
 that their aggregate momentum is conserved.  Two of these fixes can be
 used to establish a temperature gradient across a simulation domain by
-adding heat to one group of atoms (hot reservoir) and subtracting heat
-from another (cold reservoir).  E.g. a simulation sampling from the
-McDLT ensemble.  Note that the fix is applied to a group of atoms, not
-a geometric region, so that the same set of atoms is affected wherever
-they may move to.
+adding heat (energy) to one group of atoms (hot reservoir) and
+subtracting heat from another (cold reservoir).  E.g. a simulation
+sampling from the McDLT ensemble.
+
+If the {region} keyword is used, the atom must be in both the group
+and the specified geometric "region"_region.html in order to have
+energy added or subtracted to it.  If not specified, then the atoms in
+the group are affected wherever they may move to.
 
 Heat addition/subtraction is performed every N timesteps.  The {eflux}
 parameter determines the change in aggregate energy of the entire