20,21c20,21
<      *         dir/'ue0', 'ue1', 'uu0', 'uu1', 'ua0', 'ua1', 'ub0',
<      *             'ub1', 'uu0', 'uc1', 'uu0', 'ud1', 'uw0', 'uw1',
---
>      *         dir/'uw0', 'uw1', 'uu0', 'uu1', 'ua0', 'ua1', 'ub0',
>      *             'ub1', 'uu0', 'uc1', 'uu0', 'ud1', 'ue0', 'ue1',
127a128,130
>       write (*,9991) id,n1,n2,tn1,tn2
>       pause
>  9991 format (5i10)
11a12
>  9003 FORMAT (i4)
22,28c23,39
<       DO i=1,139
<         IF ((i.ne.65).and.(i.ne.138)) THEN
<           WRITE (1,9001) i
<           WRITE (1,9002) i
<         ENDIF
<       ENDDO
<       CLOSE (1)
---
>       OPEN (
>      *      access='sequential',
>      *      carriagecontrol='list',
>      *      dispose='keep',
>      *      form='formatted',
>      *      name='x',
>      *      readonly,
>      *      shared,
>      *      status='old',
>      *      unit=2
>      *     )
>   110 READ (2,9003,end=120) i
>       WRITE (1,9001) i
>       WRITE (1,9002) i
>       GO TO 110
>   120 CLOSE (1)
>       CLOSE (2)
15,16c15
<  9003 FORMAT ('(F=savefit,I=iter,R=Rad,N=ncoef',
<      *        'S=save,T=trim,Z=zone,X=exit):' $)
---
>  9003 FORMAT ('(F=savefit,I=iter,R=rad,N=ncoef,S=save,T=trim,X=exit):'$)
55,56d53
<       ELSEIF ((c.eq.'z').or.(c.eq.'Z')) THEN
<         ok = 2
36,37d35
<       zonex0 = mx
<       zoney0 = my
19,20d18
< c...Initialization
< c
22,29d19
<       zonerad = -1
< c
< c...Zone Processing
< c
<   105 CALL zoneproc(manual)
< c
< c...Loop At This Set Of Conditions
< c
53,55d42
<         IF (zonerad.gt.0) THEN
<           IF (ABS(npair-nold).lt.(MAX(npair,nold)/10)) GO TO 105
<         ENDIF
62a50,51
>       write (*,9991) niter,npair,n1,n2
>  9991 format (' Testing ', 4i10)
66a56,59
>           IF ((niter.ge.4).and.(npair.le.(MIN(n1,n2)/10))) THEN
>             fatal = 19
>             RETURN
>           ENDIF
75,79d67
<             GO TO 110
<           ELSEIF (ok.eq.1) THEN
<             GO TO 110
<           ELSEIF (ok.eq.2) THEN
<             GO TO 105
80a69
>           GO TO 110
83,86d71
< c
< c...Try Bigger Zone
< c
<       IF (zonerad.gt.0) GO TO 105
8c8
<      *            i, f(4), l(4), j1, j2
---
>      *            i, f(4), l(4)
12,34c12,20
<   100 IF (zonerad.le.0) THEN
<         j1 = 1
<         j2 = n1
<       ELSE
<         DO j1=1,n1
<           IF (x1(j1).ge.zonexl) GO TO 110
<         ENDDO
<         j1 = 1
<   110   DO j2=n1,1,-1
<           IF (x1(j2).le.zonexh) GO TO 120
<         ENDDO
<         j2 = n1
<   120   CONTINUE
<       ENDIF
<       i = (j2+1-j1)/4
<       f(1) = j1
<       l(1) = j1+i-1
<       f(2) = l(1)+1
<       l(2) = l(1)+i
<       f(3) = l(2)+1
<       l(3) = l(2)+i
<       f(4) = l(3)+1
<       l(4) = j2
---
>   100 i = n1/4
>       f(1) = 1
>       l(1) = i
>       f(2) = i+1
>       l(2) = 2*i
>       f(3) = 2*i+1
>       l(3) = 3*i
>       f(4) = 3*i+1
>       l(4) = n1
38,39d23
<       write (*,9991) zonerad,zonexl,zonexh,j1,j2
<  9991 format (' PARPF ', 5i12)
36,38d35
<       IF (zonerad.gt.0) THEN
<         IF ((yl.lt.zoneyl).or.(yh.gt.zoneyh)) GO TO 110
<       ENDIF
