moved file close until after zisodisp is written

src/write_global_output.f90

@@ -116,13 +116,13 @@ if (iproc.eq.0) then
 ! Line below uncommented by dwhipp - 12/09 
                  osolve%pressure(i),        & 
 ! Line below added by dwhipp - 12/09 
-                 osolve%spressure(i),  &
+                 osolve%spressure(i),       &
                  osolve%crit(i),            &
                  osolve%e2d(i),             &
                  osolve%eviscosity(i),      & 
-                 osolve%is_plastic(i),               & 
+                 osolve%is_plastic(i),      & 
                  ov%whole_leaf_in_fluid(i), &
-                  i=1,osolve%nleaves)
+                 i=1,osolve%nleaves)
    ! write octree information
    write (9) (osolve%octree(i),i=1,osolve%octree(3))
    ! write bad face information
@@ -181,13 +181,14 @@ if (iproc.eq.0) then
                 cl%press(i),  &
                 cl%tag(i),    &
                 i=1,cl%np)
-   close (9)
 
    if (params%isobc) then 
      ! write isostasy basal displacement array - dwhipp 11/09 
        write (9) ((zisodisp(i,j),j=1,2**params%levelmax_oct+1),i=1,& 
                  2**params%levelmax_oct+1) 
-   end if 
+   end if
+
+   close (9)
 
 !   do i=1,osolve%nnode
 !      if (abs(osolve%x(i)-0.5d0)<1.d-6 .and. osolve%z(i)<0.08d0) then