diff --git a/schainf/Ffiles/fitacf.f b/schainf/Ffiles/fitacf.f
index df699eb..d6f6f39 100644
--- a/schainf/Ffiles/fitacf.f
+++ b/schainf/Ffiles/fitacf.f
@@ -485,6 +485,9 @@ c     between k and the magnetic field (radians)
         sum2=sum2+fi(i)*yi(i)
       end do
       dl=ak**2*dlf*te/densmks
+c      write(*,fmt='("Before YE")')
+c      write(*,*) imode
+c      call exit
 
       if(imode.eq.1.or.imode.eq.3) then
       
@@ -493,6 +496,8 @@ c     between k and the magnetic field (radians)
          phi=(omegae/ak)/(sqrt(2.0)*vte)
          psie=(ven/ak)/(sqrt(2.0)*vte)
          ye=y_electron(thetae,phi,psie,alpha)
+	 write(*,fmt='("AFTER YE")')
+c         call exit
   
       else if(imode.eq.2) then
 c
@@ -503,7 +508,8 @@ c
          alpha2=abs(pi/2.0-alpha)*180.0/pi
 c         write(*,*) "ye: ", ye
          call collision(densmks, te, freq, alpha2, ye)
-c         write(*,fmt='(" geobfield: time is before earliest model.")')
+c         write(*,fmt='("AFTER COLLISION")')
+c         call exit
          ye=ye*omega+cmplx(0.0,1.0)        
 
       end if
@@ -512,12 +518,12 @@ c         write(*,fmt='(" geobfield: time is before earliest model.")')
       p=(cabs(ye))**2*real(sum2)+cabs(sum1+cmplx(0.0,dl))**2*real(ye)
       p=p/(cabs(yed+sum1))**2
       spect1=p*2.0e0/(omega*pi)
-      write(*,*) "spect1 DONE "
+      write(*,*) "spect1:",spect1
       return
       end
 
-      subroutine acf2(wl, tau, te1, ti1, fi1, ven1, vin1, wi1, nion1,
-     &     alpha1, dens1, bfld1, acf)
+      subroutine acf2(wl, tau, te1, ti1, fi1, ven1, vin1, wi1,
+     &     alpha1, dens1, bfld1, acf, nion1)
 c
 c     computes autocorrelation function for given plasma parameters
 c     by integrating real spectrum
@@ -531,33 +537,44 @@ c
       real pi
       integer nion1
       integer wi1(nion1)
-      integer i,j,k
+      integer i,j,k,imode
+      common /mode/imode
 c
+c      write(*,*) "INITIAL acf:",wl,tau,te1,ti1,fi1,ven1,vin1,wi1,alpha1
+      write(*,*) "INITIAL acf:",dens1, bfld1, acf, nion1
+c      write(*,fmt='("INIT")')
       pi=4.0*atan(1.0)
 c
 c     copy arguments to common block
 c
       ak=2.0*pi/wl
       imode=2
+      write(*,*) "imode:",imode
 
       nion=nion1
       alpha=alpha1
       te=te1
       ven=ven1
+c      write(*,fmt='("INIT2")')
       do i=1,nion
+c      write(*,fmt='("INIT2.5")')
         ti(i)=ti1(i)
         fi(i)=fi1(i)
         vin(i)=vin1(i)
         wi(i)=wi1(i)
       end do
+c      write(*,fmt='("INIT3")')
       dens=dens1
       bfld=bfld1
 
 c      write(*,*) wl,alpha1,bfld1,dens
 c      call exit
-
+c      write(*,fmt='("Before Gauss")')
       call gaussq(tau,acf)
 
+      write(*,*) "FINAL acf:",acf
+
+c      write(*,fmt='("After Gauss")')
       return
       end
 
@@ -590,7 +607,7 @@ c      write(*,*) leniw,lenw
       nrmom=0
       ksave=0
       momcom=0
- 
+      write(*,*) "Before qc25f:",acf,imode
 c     much faster, more robust
 c      write(*,*) "acf_in: ",acf
       call qc25f(spect1,a,b,tau,integr,nrmom,maxp1,ksave,acf,
@@ -598,7 +615,7 @@ c      write(*,*) "acf_in: ",acf
 c      write(*,*) "acf_out: ",acf
 c      call qawf(spect1,a,tau,integr,epsabs,acf,abserr,neval,
 c     &     ier,limlst,lst,leniw,maxp1,lenw,iwork,work)
- 
-      
+      write(*,*) "After qc25f:",acf
+c      call exit
       return
       end