check4nans.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	check4nans (meqn, mbc, mx, my, q, t, ichecknan)
	Check for NANs in solution q. More...

Function/Subroutine Documentation

check4nans()

subroutine check4nans	(	integer, intent(in)	meqn,
		integer, intent(in)	mbc,
		integer, intent(in)	mx,
		integer, intent(in)	my,
		real(kind=8), dimension(meqn,1-mbc:mx+mbc,1-mbc:my+mbc), intent(in)	q,
		real(kind=8), intent(in)	t,
		integer, intent(in)	ichecknan
	)

Check for NANs in solution q.

Definition at line 4 of file check4nans.f90.

    implicit none
    
    ! Input
    integer, intent(in) ::  meqn, mbc, mx, my, ichecknan
    real(kind=8), intent(in) :: t, q(meqn,1-mbc:mx+mbc,1-mbc:my+mbc)

    ! Locals
    integer :: i, j, m

    ! print *,'Checking for NANs at ichecknan = ',ichecknan
    ! print *,'  mx = ',mx,'  my = ',my,'  meqn = ',meqn
      
    do i=1-mbc,mx+mbc
        do j=1-mbc,my+mbc
            do m=1,meqn
                if (.not. (q(m,i,j) == q(m,i,j))) then
                    ! true if q(i,j,m) = NAN
                    print *, 'SOLUTION ERROR --- ABORTING CALCULATION'
                    print *, 'At ichecknan = ',ichecknan
                    print *, '   mx,my,t:',mx,my,t
                    print *, '   m,i,j:',m,i,j
                    print *, '   q(m,i,j) = ',q(m,i,j)
                    stop
                endif
            enddo
        enddo
    enddo

    ! Uncomment the next line if desired when debugging:
    ! print *,'No NANs at ichecknan = ',ichecknan,' at t = ',t