Added a few comments

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