anchor.diff:    write(num, '(i10.10)') n
   
     ! Output the data to a file
  -  open(unit=10, status='replace', file='output/'//name//num, &
  +  open(unit=10, status='replace', file=name//num, &
         form='unformatted')
     write(10) array(1:nx+xstep, 1:ny+ystep)
     close(10)
  
positive.code:   call evaluate_zeta(zeta, u, v, nx, ny, layers, dx, dy)
    if (debug_level .ge. 4) then
      call write_output_3d(zeta, nx, ny, layers, 1, 1, &
        n, 'snap.zeta.')
    end if
  
    ! Calculate dhdt, dudt, dvdt at current time step
  
negative.code:   ! print *, maxval(abs(etastar))
    if (debug_level .ge. 4) then
      call write_output_2d(etastar, nx, ny, 0, 0, &
        n, 'output/snap.eta_star.')
    end if
  
    ! Prevent barotropic signals from bouncing around outside the
  