anchor.diff:      ! h2av = 0.0
   
     end if
  +
  +
   
     return
   end subroutine maybe_dump_output
  
positive.code:   write(num, '(i10.10)') n
  
    ! Output the data to a file
    open(unit=10, status='replace', file='output/'//name//num, &
        form='unformatted')
    write(10) array(1:nx+xstep, 1:ny+ystep)
    close(10)
  
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
  