diff --git a/write_stats.f90 b/write_stats.f90
index 56b651f..26b1a02 100644
--- a/write_stats.f90
+++ b/write_stats.f90
@@ -19,7 +19,13 @@ subroutine write_tp_field
  integer :: i, j, k, ii
 
  if (fdm_sw.eq.0) then
-   write(tp_field,897) time,nx/spx,ny/spy,nz/spz
+
+   if(mod(nx,spx).eq.0) then
+    write(tp_field,897) time,nx/spx,ny/spy,nz_all/spz
+   else
+    write(tp_field,897) time,nx/spx+1,ny/spy+1,nz_all/spz+1
+   endif
+
    do k=1,nz_all,spz
     do j=1,ny,spy
      do i=1,nx,spx
@@ -27,13 +33,19 @@ subroutine write_tp_field
      enddo
     enddo
    enddo
- else
-   write(tp_field,897) time,nx/svf,ny/svf,nz_all/svf
 
-   do ii=1,nz_all,svf
-    do j=1,ny,svf
-     do i=1,nx,svf
-           WRITE(tp_field,'(7e30.20)') REAL(i)*hx,REAL(j)*hx,REAL(ii)*hx, &
+ else
+
+   if(mod(nx,spx).eq.0) then
+    write(tp_field,897) time,nx/spx,ny/spy,nz_all/spz
+   else
+    write(tp_field,897) time,nx/spx+1,ny/spy+1,nz_all/spz+1
+   endif
+
+   do ii=1,nz_all,spz
+    do j=1,ny,spy
+     do i=1,nx,spx
+           write(tp_field,'(7e20.10)') real(i)*hx,real(j)*hx,real(ii)*hx, &
                                  (1.-y1(2,i,j,ii)/y1(1,i,j,ii)),u_(i,j,ii), &
                                   v_(i,j,ii),w_(i,j,ii)
      enddo