--- trunk/OOPSE/libmdtools/calc_reaction_field.F90	2004/01/06 19:49:18	901
+++ trunk/OOPSE/libmdtools/calc_reaction_field.F90	2004/01/30 15:01:09	999
@@ -129,7 +129,7 @@ contains
        if (rij.lt.rt) then
           taper = 1.0d0
        else
-          write(*,*) 'rf in taper region'
+  !        write(*,*) 'rf in taper region'
           taper = (rrf + 2.0d0*rij - 3.0d0*rt)*(rrf-rij)**2/ ((rrf-rt)**3)
        endif
        
@@ -234,7 +234,7 @@ contains
        rfpot = rfpot - 0.5d0 * pre * mu1 * &
             (rf(1,a1)*u_l(1,a1) + rf(2,a1)*u_l(2,a1) + rf(3,a1)*u_l(3,a1))
     endif
-    
+ 
     return
   end subroutine reaction_field_final
   
@@ -275,7 +275,7 @@ contains
        if (rij.lt.rt) then
           dtdr = 0.0d0
        else
-          write(*,*) 'rf correct in taper region'
+ !         write(*,*) 'rf correct in taper region'
           dtdr = 6.0d0*(rij*rij - rij*rt - rij*rrf +rrf*rt)/((rrf-rt)**3)
        endif