diff --git a/src/forces.f90 b/src/forces.f90
index 44915cf6161cb50ded0bebc219edd4d9b96dc501..bb68e9f06000df046953e79e1205b1d6f519b755 100644
--- a/src/forces.f90
+++ b/src/forces.f90
@@ -451,6 +451,9 @@ contains
 
     if (itime.eq.1) then
       do iv=1,nvol
+         ! if (nrank.eq.0) then 
+         !    open(12,file="Body1.dat",status='unknown',form='formatted')
+         ! endif
          if ((nrank .eq. 0).and.(record_var.eq.1)) then
             write(filename,"('forces.dat',I1.1)") iv
             open(38+(iv-1),file=filename,status='unknown',form='formatted')
@@ -1039,6 +1042,8 @@ contains
          write(38+(iv-1),*) t,dra1,dra2,dra3, sum(tdiffx), sum(tdiffx2), tp1, -mom1, -sum(tunstx(:)), -sum(tconvx(:)), -sum(tconvx2(:))
          !  write(*,*) 'written to file number', 38+(iv-1), t, dra1,dra2,dra3
           call flush(38+(iv-1))
+
+         !  write(12 ,*) t, position(1), position(2), position(3), orientation(1), orientation(2), orientation(3), orientation(4), linearVelocity(1), linearVelocity(2), linearVelocity(3), angularVelocity(1), angularVelocity(2), angularVelocity(3), linearAcceleration(1), linearAcceleration(2), linearAcceleration(3)
        endif
       !  if (mod(itime, ioutput).eq.0) then
       !     if (nrank .eq. 0) then
diff --git a/src/xcompact3d.f90 b/src/xcompact3d.f90
index d757c64c02638fccc76d8f6335960ac25ff3ee42..8634f356f71cd3f768ee8cc63e101189b2b106ac 100644
--- a/src/xcompact3d.f90
+++ b/src/xcompact3d.f90
@@ -29,14 +29,17 @@ program xcompact3d
   iounit = 135
   !Print forces out on ellip
   if ((nrank==0).and.(force_csv.eq.1)) then 
-   open(unit=20, file='force_out.dat', status='replace',action='write',iostat=ierr)
-   if (ierr /= 0) then
-      print *, 'Error opening file.'
-      stop
-   end if
-   write(*,*) 'Outputting forces'
-  end if 
+   open(unit=20, file='force_out.dat', status='unknown',form='formatted')
+   ! if (ierr /= 0) then
+   !    print *, 'Error opening file.'
+   !    stop
+   ! end if
+   write(*,*) 'Outputting forces' 
+  end if
 
+  if (nrank==0) then
+   open(unit=12, file='body.dat1', status='unknown', form='formatted')
+  endif
 !   do i = 1,100
 !    x(i) = i
 !   enddo
@@ -134,10 +137,10 @@ program xcompact3d
         endif
 
         if ((nrank==0).and.(force_csv.eq.1)) then
-         open(unit=20, file='force_out.dat', action='write')
-         write(20, *) linearForce(1), linearForce(2), linearForce(3), '\n'
+         ! open(unit=20, file='force_out.dat', action='write')
+         write(20, *) linearForce(1), linearForce(2), linearForce(3)
          write(*,*) 'Writing forces', linearForce(1), linearForce(2), linearForce(3)
-         close(20)
+         flush(20)
         endif 
         
         
@@ -160,6 +163,13 @@ program xcompact3d
 
          orientation = orientation_1
          angularVelocity = angularVelocity_1
+
+
+
+        if (nrank==0) then
+         write(12 ,*) t, position(1), position(2), position(3), orientation(1), orientation(2), orientation(3), orientation(4), linearVelocity(1), linearVelocity(2), linearVelocity(3), angularVelocity(2), angularVelocity(3), angularVelocity(4), drag, lift, lat, xtorq, ytorq, ztorq
+         flush(12)
+        endif
       !   endif 
 
          if ((nrank==0).and.(mod(itime,ilist)==0)) then