COMMON /GRUNGE/I,TD,ZO,EMA,PI,W,WC
      COMMON /FIT/X,Y,SIGMAY,NPTS,MODE,A,SIGMAA,B,SIGMAB,R,SCALEX,SCALEY
      DOUBLE PRECISION X,Y,P1,P2,SIGMAY,A,SIGMAA,B,SIGMAB,R
      DOUBLE PRECISION TD,ZO,EMA,PI,W,WC,SP2,SCALEX,SCALEY
      INTEGER NPTS,MODE,I,J,K
      DIMENSION X(50),Y(50),SIGMAY(50)
      WC=2./DSQRT(5D-12)
      PI=3.1415926535
1     TYPE 100
      ACCEPT 101,NPTS
      IF (NPTS.EQ.0) GOTO 999
      TYPE 102
      ACCEPT 103,MODE
      TYPE 500
      ACCEPT 501,J
      TYPE 502
      ACCEPT 503,SCALEX,SCALEY
      TYPE 104
      DO 10 I=1,NPTS
      ACCEPT 105,P1,P2,SP2
      W=P1*PI*2.D3
      TD=P2
      ZO=P2
      EMA=P2
      GOTO (51,52,53,54),J
51    CALL GRUNG1
      GOTO 10
52    CALL GRUNG2
      GOTO 10
53    CALL GRUNG3
      GOTO 10
54    CALL GRUNG4
10    CONTINUE
      DO 600 K=1,NPTS
600   TYPE 601,X(K),Y(K)
      CALL LINFIT
      TYPE 200,A,SIGMAA
      TYPE 201,B,SIGMAB
      TYPE 202,R
      TYPE 203
      GOTO 1
999   STOP
100   FORMAT(' NUMBER OF POINTS='$)
101   FORMAT(I)
102   FORMAT(' MODE='$)
103   FORMAT(I)
104   FORMAT('0INPUT DATA: INDEPENDENT/DEPENDENT[/SIGMA(DEPENDENT)]'/)
105   FORMAT(3D)
200   FORMAT(' A= ',1PD20.5,9X,'SIGMAA=',1PD20.5)
201   FORMAT(' B= ',1PD20.5,9X,'SIGMAB=',1PD20.5)
202   FORMAT(' CORRELATION= ',F15.7)
203   FORMAT('1')
500   FORMAT(' GRUNGE= '$)
501   FORMAT(I)
502   FORMAT(' SCALEX/SCALEY= ')
503   FORMAT(2D)
601   FORMAT(1X,1P2D15.5)
      END