From a3c5931a088f3aca307951c48958cd02e7d25fd6 Mon Sep 17 00:00:00 2001
From: Paul Romano <paul.k.romano@gmail.com>
Date: Sun, 21 Apr 2013 22:58:24 -0400
Subject: [PATCH] Fixed bug in use of free gas scattering for H-1.

---
 src/physics.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/physics.F90 b/src/physics.F90
index 76c4141cdd..3000a2dbce 100644
--- a/src/physics.F90
+++ b/src/physics.F90
@@ -790,7 +790,7 @@ contains
     kT = nuc % kT
 
     ! Check if energy is above threshold
-    if (p % E >= FREE_GAS_THRESHOLD * kT) then
+    if (p % E >= FREE_GAS_THRESHOLD * kT .and. nuc % awr > ONE) then
       v_target = ZERO
       return
     end if