diff --git a/src/write_global_output.f90 b/src/write_global_output.f90
index 7aa5c861c1092fec469190cfd543a3334c52e34b..a7fbe112edef995f7102b329f5f18c4a51c86b63 100644
--- a/src/write_global_output.f90
+++ b/src/write_global_output.f90
@@ -76,7 +76,6 @@ call mpi_comm_size (mpi_comm_world,nproc,ierr)
 call mpi_comm_rank (mpi_comm_world,iproc,ierr)
 
 if (iproc.eq.0) then
-   print *,'Entered write_global_output'
    call int_to_char (cistep,4,istep)
    select case (outputtype)
       case ('debug')
@@ -88,7 +87,6 @@ if (iproc.eq.0) then
       case default
          call stop_run ('pb with argument in write_global_output$')
    end select
-   print *,'Case chosen'
    ! write array lengths
    write (9) osolve%octree(3), &
                osolve%nnode,     &
@@ -97,7 +95,6 @@ if (iproc.eq.0) then
                osolve%nlsf,      &
                cl%np,            &
                current_time
-   print *,'Array lengths written'
    ! write info on octree solve nodes (x,y,z,u,v,w,lsf,temp)
    write (9) (osolve%x(i),                         &
                 osolve%y(i),                         &  
@@ -115,7 +112,6 @@ if (iproc.eq.0) then
                 osolve%kfix((i-1)*3+3),              &
                 osolve%kfixt(i),                     &
                 i=1,osolve%nnode)
-   print *,'osolve nodal values written'
    ! write icon array
    write (9) ((osolve%icon(k,i),k=1,8),   &
 ! Line below uncommented by dwhipp - 12/09 
@@ -128,13 +124,10 @@ if (iproc.eq.0) then
                  osolve%is_plastic(i),      & 
                  ov%whole_leaf_in_fluid(i), &
                  i=1,osolve%nleaves)
-   print *,'osolve elemental values written'
    ! write octree information
    write (9) (osolve%octree(i),i=1,osolve%octree(3))
-   print *,'octree info written'
    ! write bad face information
    write (9) ((osolve%iface(k,i),k=1,9),i=1,osolve%nface)
-   print *,'bad faces written'
    ! write void information
    write (9) (vo%node(i),    &
                 vo%leaf(i),    &
@@ -142,17 +135,13 @@ if (iproc.eq.0) then
                 vo%rtf(i),     &
                 vo%influid(i), &
                 i=1,osolve%nnode)
-   print *,'void written'                
    ! write bad faces information
    write (9) (vo%face(i),i=1,osolve%nface)
-   print *,'bad faces written'
    ! write surface information (r,s,x,y,z,xn,yn,zn)
    do is=1,osolve%nlsf
-      print *,'In surface loop'
       write (9) surface(is)%nsurface,          &
                   surface(is)%activation_time, &
                   surface(is)%nt
-      print *,'Basic info written'
       select case (outputtype)
       case ('debug')
       write (9) (surface(is)%r(i),  &
@@ -178,9 +167,7 @@ if (iproc.eq.0) then
                    surface(is)%w(i),  &
                    i=1,surface(is)%nsurface)
       end select
-      print *,'Surface info written'
       write (9) (surface(is)%icon(1:3,i),i=1,surface(is)%nt)
-      print *,'Surface connectivity written'
    enddo
    ! write cloud information
    write (9) (cl%x(i),      &
@@ -195,21 +182,16 @@ if (iproc.eq.0) then
                 cl%press(i),  &
                 cl%tag(i),    &
                 i=1,cl%np)
-      print *,'Cloud written'
 
    if (params%isobc) then 
      ! write isostasy basal displacement array - dwhipp 11/09 
        write (9) 2**params%levelmax_oct
-       print *,'isobc size written'
        write (9) ((zi%zisodisp(i,j)+surface(surface(is)%nsurface)%sp01,&
                  j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1) 
-       print *,'zisodisp written'
    end if
 
    close (9)
 
-   print *,'And we out...'
-
 !   do i=1,osolve%nnode
 !      if (abs(osolve%x(i)-0.5d0)<1.d-6 .and. osolve%z(i)<0.08d0) then
 !         write (1234,'(5f30.15)') osolve%y(i),osolve%z(i),osolve%v(i),osolve%w(i),ov%temporary_nodal_pressure(i)