anchor.diff:      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
  +          write(17, "(A, I0)") "NaN detected at time step ", n
             stop 1
           end if
         end do
  
positive.code:     end if
    end do
  
    write(17, "(A, I0)"), 'Warning: maximum SOR iterations exceeded at time step ', n
  
    return
  end subroutine SOR_solver
  
negative.code: 
    if (num_procs .ne. nProcX * nProcY) then
      if (myid .eq. 0) then
         write(17, "(A)") "number of processors in run command must equal nProcX * nProcY - fix this and try again"
         write(17, "(A, I0)") 'num_procs = ', num_procs
         write(17, "(A, I0)") 'nProcX = ', nProcX
         write(17, "(A, I0)") 'nProcY = ', nProcY
      end if
      stop 1
    end if
  