anchor.diff:        do i = 1, nx
           if (data(i,j,k) .ne. data(i,j,k)) then
             write(17, "(A, I0)") "NaN detected at time step ", n
  +          call clean_stop(n, .FALSE.)
             stop 1
           end if
         end do
  
positive.code:        write(17, "(A, I0)") 'nProcX = ', nProcX
         write(17, "(A, I0)") 'nProcY = ', nProcY
      end if
      stop 1
    end if
  
  
negative.code:     ! Check that depth is positive - it must be greater than zero
      if (minval(depth) .lt. 0) then
        write(17, "(A)") "Depths must be positive."
        call clean_stop(0, .FALSE.)
        stop 1
      end if
    end if
  