anchor.diff:      call read_input_fileH_2D(initEtaFile, eta, 0.d0, nx, ny)
       ! Check that depth is positive - it must be greater than zero
       if (minval(depth) .lt. 0) then
  -      write(17, "(A)"), "Depths must be positive."
  +      write(17, "(A)") "Depths must be positive."
         stop 1
       end if
     end if
  
positive.code:     do j = 1, ny
        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
            stop 1
          end if
        end do
  
negative.code:           hnew(i, j, k) = hmin
            counter = counter + 1
            if (counter .eq. 1) then
              write(17, "(A, I0)") &
                  "Layer thickness dropped below hmin at time step ", n
            end if
          end if
  