readngh_utils.f

Programs in the irregular grid design package described in this manual are used to carry out five ma

        subroutine ReadNghHdr(unit, error,
     +          mrec, nbtot, numrec, nbtotr,
     +		xmax, ymax, xmin, ymin)
	implicit none
	integer unit
	integer error
	integer mrec
	integer nbtot
	integer numrec
	integer nbtotr
	real	xmax, ymax, xmin, ymin
	character*256 msg

      read (unit,*,err=9999) NUMREC
      if (NUMREC .gt. MREC) then
	write(msg,'(a,i7,a,i7)') 
     +    'Too many nodes (',NUMREC,') in file: max allowed=',MREC
	call PigPutMessage(msg)
	goto 9999
      endif
c
      read (unit,*,err=9999) NBTOTR
      if (NBTOTR .gt. NBTOT) then
	write(msg,'(a,i7,a,i7)') 
     +    'Too many neighbours (',NBTOTR,') in file: max allowed='
     +    ,NBTOT
	call PigPutMessage(msg)
	goto 9999
      endif
      write(msg,'(a,i7,a,i4,a)') 
     +		'Reading file with',NUMREC,' nodes,',
     +		NBTOTR,' neighbours.'
      call PigPutMessage(msg)

      read (unit,*,err=9999) XMAX,YMAX,XMIN,YMIN

	error = 0
	return

9999	continue
	call PigPutMessage(
     +    'Error reading neighbour file header - in s/r ReadNghHdr')
	error = 9999
	end

	subroutine ReadNghNode(unit, error, x, y, code,
     +          dep, neighs, maxngh)
	implicit none
	integer unit
	integer error
	real	x
	real	y
	integer code
	real 	dep
	integer maxngh
	integer	neighs(maxngh)

	integer	irec
	integer j

      read (unit,*,end=999,err=9999) IREC, x, y, code, dep, 
     +		(neighs(j),j=1,maxngh)
	error = 0
	return

9999	continue
	call PigPutMessage('Error reading node - in s/r ReadNghNode')
	error = 9999
	return

999	continue
	call PigPutMessage(
     +          'End of file - reading node - in s/r ReadNghNode')
	error = 999
	end