PROGRAM checkwat

c	Reads output of second pass through cover2 and lists
c	those solvents which need to checked

	CHARACTER*103	filnam, line

c	Get files

	WRITE (6, 1000) 'Enter input file name:'
1000	FORMAT (a)
	READ (5, 1000) filnam
	len = length(filnam)
	WRITE (6, 1000) filnam(1:len)	
	OPEN (unit=1, form='FORMATTED', status='OLD', readonly,
	1	file=filnam(1:len))

	WRITE (6, 1000) 'Enter output file name:'
	READ (5, 1000) filnam
	len = length(filnam)
	WRITE (6, 1000) filnam(1:len)
	OPEN (unit=2, form='FORMATTED', status='UNKNOWN', 
	1	file=filnam(1:len))

	WRITE (6, 1000) 'Enter levels for Fo-Fc and 2Fo-Fc maps:'
	READ (5, 1400) crit1, crit2
1400	FORMAT (2f5.1)
	WRITE (6, 1400) crit1, crit2

	WRITE (6, 1000) 'Enter level for Bs:'
	READ (5, 1400) bcrit
	WRITE (6, 1400) bcrit

c	Loop through file

	DO WHILE (.TRUE.)
10	  READ (1, 1000, end=99) line
	  IF (line(66:68) .NE. 'WAT' .AND.
	1	line(66:68) .NE. 'TIP') GOTO 10
	  READ (line, 1001) b, df, f
1001	  FORMAT (30x, f10.3, 48x, 3f8.1)
	  IF (df .LE. crit1 .OR. f .LE. crit2 .OR.
	1	b .GT. bcrit) THEN
	    WRITE (2, 1000) line
	  endIF
	endDO

99	CLOSE (unit=1)
	CLOSE (unit=2)

	end
      INTEGER FUNCTION length(line)

      CHARACTER*103   line

      DO i = 103,1,-1
        IF (line(i:i) .NE. ' ') goto 9
      endDO

9     length = i

      return
      end