diff --git a/write_stats.f90 b/write_stats.f90
index 283f4f4..dc60c77 100644
--- a/write_stats.f90
+++ b/write_stats.f90
@@ -31,6 +31,7 @@ subroutine write_tp_field
      do i=1,nx,spx
         write(tp_field,'(6e18.9)') dx*real(i),dy*real(j),dz*real(k),u_(i,j,k),v_(i,j,k),w_(i,j,k)   ! u',v',w'
      enddo
+     call flush(tp_field)
     enddo
    enddo
 
@@ -49,13 +50,12 @@ subroutine write_tp_field
                                  (1.-y1(2,i,j,ii)/y1(1,i,j,ii)),u_(i,j,ii), &
                                   v_(i,j,ii),w_(i,j,ii)
      enddo
+     call flush(tp_field)
     enddo
    enddo
 
  endif
 
- call flush(tp_field)
-
 897      format('ZONE T="time= ',f10.5,'" I=  ',i4,' J=  ',i4,' K=  ',i4,' F=POINT')
 
 end subroutine write_tp_field