anchor.diff:        write(num, '(i10.10)') n
   
         ! Output the data to a file
  -      open(unit = 10, status='replace', file='output/snap.h.'//num, &
  +      open(unit=10, status='replace', file='output/snap.h.'//num, &
             form='unformatted')
         write(10) h(1:nx, 1:ny, :)
         close(10)
  -      open(unit = 10, status='replace', file='output/snap.u.'//num, &
  +      open(unit=10, status='replace', file='output/snap.u.'//num, &
             form='unformatted')
         write(10) u(1:nx+1, 1:ny, :)
         close(10)
  -      open(unit = 10, status='replace', file='output/snap.v.'//num, &
  +      open(unit=10, status='replace', file='output/snap.v.'//num, &
             form='unformatted')
         write(10) v(1:nx, 1:ny+1, :)
         close(10)
         if (.not. RedGrav) then
  -        open(unit = 10, status='replace', file='output/snap.eta.'//num, &
  +        open(unit=10, status='replace', file='output/snap.eta.'//num, &
               form='unformatted')
           write(10) eta(1:nx, 1:ny)
           close(10)
         end if
   
         if (DumpWind .eqv. .true.) then
  -        open(unit = 10, status='replace', file='output/wind_x.'//num, &
  +        open(unit=10, status='replace', file='output/wind_x.'//num, &
               form='unformatted')
           write(10) wind_x(1:nx+1, 1:ny)
           close(10)
  -        open(unit = 10, status='replace', file='output/wind_y.'//num, &
  +        open(unit=10, status='replace', file='output/wind_y.'//num, &
               form='unformatted')
           write(10) wind_y(1:nx, 1:ny+1)
           close(10)
  
positive.code:   double precision array_small(nx, ny+1, layers)
  
    if (name.ne.'') then
      open(unit = 10, form='unformatted', file=name)
      read(10) array_small
      close(10)
  
  
negative.code:   integer nx, ny, layers, n, i, j, k
    double precision data(0:nx+1, 0:ny+1, layers)
  
    do k = 1, layers
      do j = 1, ny
        do i = 1, nx
          if (data(i,j,k) .ne. data(i,j,k)) then
            ! write a file saying so
            open(unit=10, file='NaN detected.txt', action="write", &
  