diff --git a/src/smooth_pressures.f90 b/src/smooth_pressures.f90
index b655c25c2df72ea93a566c2f8de1028e5446cfb0..06697f950347c9e428dc7f248aa61af58e73d52d 100644
--- a/src/smooth_pressures.f90
+++ b/src/smooth_pressures.f90
@@ -79,6 +79,10 @@ allocate (field(osolve%nnode))
 field=0
 countnode=0
 
+! A few comments for reference by dwhipp - 12/09 
+! This loop goes through each element, loops through the associated nodes and 
+! adds the elemental pressure to field (the nodal sum). countnode tracks the 
+! number of times a pressure value is contributed 
 do i=1,osolve%nleaves
    do j=1,8
       k=osolve%icon(j,i)
@@ -87,12 +91,16 @@ do i=1,osolve%nleaves
    end do
 end do
 
+! This loop goes through each node and divides the field value from above by 
+! countnode, to yield a nodal value that averages the elemental pressures from 
+! each element that shares that node 
 do i=1,osolve%nnode
    if (countnode(i).gt.0.d0) then
       field(i)  = field(i) / countnode(i)
    end if
 end do
 
+! This stores the averaged nodal pressures
 ov%temporary_nodal_pressure=field
 
 do ie=1,osolve%nleaves