diff --git a/src/pflotran/richards.F90 b/src/pflotran/richards.F90
index 9156b6a6c..043082ab4 100644
--- a/src/pflotran/richards.F90
+++ b/src/pflotran/richards.F90
@@ -1798,6 +1798,8 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)
   PetscReal :: pressure_min
   PetscReal :: well_inj_water
   PetscReal :: Dq, dphi, v_darcy, ukvr
+  PetscInt  :: count, itest1
+  PetscReal :: qflx_pf_sum, flow_sum1, flow_sum2, ss_flow_fluxes_sum, qsrc_sum

   Mat, parameter :: null_mat = tMat(0)

@@ -1822,6 +1824,14 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)
   ! Source/sink terms -------------------------------------
   source_sink => patch%source_sink_list%first
   sum_connection = 0
+  count = 0
+  qsrc_sum=0.d0
+  itest1 = 0
+  flow_sum1=0.d0
+  flow_sum2=0.d0
+  ss_flow_fluxes_sum=0.d0
+
+  write(*,*)''
   do
     if (.not.associated(source_sink)) exit

@@ -1835,6 +1845,8 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)

 #if defined(CLM_PFLOTRAN) || defined(CLM_OFFLINE)
       qsrc = qflx_pf_p(local_id)
+      count = count + 1
+      qsrc_sum=qsrc_sum+qflx_pf_p(local_id)
 #endif
       if (source_sink%flow_condition%itype(1)/=HET_VOL_RATE_SS .and. &
           source_sink%flow_condition%itype(1)/=HET_MASS_RATE_SS .and. &
@@ -1862,6 +1874,8 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)
         case(HET_MASS_RATE_SS)
                                                            ! kg/sec -> kmol/sec
           qsrc_mol = source_sink%flow_aux_real_var(ONE_INTEGER,iconn)/FMWH2O
+          itest1 = itest1 + 1
+          flow_sum1=flow_sum1+source_sink%flow_aux_real_var(ONE_INTEGER,iconn)

         case(WELL_SS) ! production well, SK 12/19/13
           ! if node pessure is lower than the given extraction pressure,
@@ -1898,6 +1912,7 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)
             endif
           endif
       end select
+      write(*,*)'In richards.F90: ',count, local_id, qsrc, qsrc_mol*FMWH2O

       if (option%compute_mass_balance_new) then
         ! need to added global auxvar for src/sink
@@ -1921,6 +1936,10 @@ subroutine RichardsResidualSourceSink(r,realization,ierr)
     enddo
     source_sink => source_sink%next
   enddo
+  write(*,*)'qsrc_sum  ',qsrc_sum
+  write(*,*)'flow_sum1 ',flow_sum1
+  write(*,*)'itest1    ',itest1
+  stop