PROGRAM ERRBIN
C This is the error estimation program for fits to the visibilities
C from a binary star during a given night, as found by BINFIT.
C     Version 1.0     23 Jan 1991     J.T. Armstrong
C

      INCLUDE 'BINFIT.INC'

      INTEGER*4 IER, L, GETIN, VLBOPE, WHERE

      INCLUDE 'MODEL.INC'

      CALL PUTOUT('Program BINERR V1.0')
      CALL PUTOUT('Estimates uncertainties for binary star data.')
      CALL PUTOUT(' ')

C Zero some arrays:

      DO IG = 1, MXFILT
          LAMBDA0(IG) = 0.0
          DLAMBDA(IG) = 0.0
      END DO

C Get input data file name, star number; calculate UDATA, VDATA

      CALL READFILE
      FRSTRY = .TRUE.

C Reset all the VARY flags

      DO IC = 1, NCOMP
          DO IG = 1, NFILT
              DO J = 1, 8
                  VARY(IC,IG,J) = .FALSE.
              END DO
          END DO
      END DO

C Get final model

   10 CALL GETMNAME
      IF (NEWMOD) THEN
          WRITE(OUTC,*) ' Don''t use a new model; use the fitted model.'
          GO TO 10
      END IF
      CALL ROLDMOD
      CALL CADJMOD

C Save these parameters

      DO IC = 1, NCOMP
          DO IG = 1, NFILT
              PARMSTRT(IC,IG,1) = SEP    (IC,IG)
              PARMSTRT(IC,IG,2) = POSANG (IC,IG)
              PARMSTRT(IC,IG,3) = DIAM   (IC,IG)
              PARMSTRT(IC,IG,4) = DRATIO (IC,IG)
              PARMSTRT(IC,IG,5) = DELMAG (IC,IG)
              PARMSTRT(IC,IG,6) = COLOR  (IC,IG)
              PARMSTRT(IC,IG,7) = FLUX   (IC,IG)
              PARMSTRT(IC,IG,8) = FLUXTOT(IC,IG)
          END DO
      END DO

C Get initial chi squared

D     WRITE(OUTC,*) ' ERRBIN: getting initial chi squared'
      CALL BINVIS
      CALL CAGRFAC
      TCHISQ0 = TCHISQ
      DO IG = 1, NFILT
          FCHISQ0(IG) = FCHISQ(IG)
      END DO

C Get parameters to vary

      CALL VARYMOD

C Determine the uncertainties

      CALL CHIBNDRY

C Reinstate the original model parameters and write results to a file

      DO IC = 1, NCOMP
          DO IG = 1, NFILT
              SEP    (IC,IG) = PARMSTRT(IC,IG,1)
              POSANG (IC,IG) = PARMSTRT(IC,IG,2) 
              DIAM   (IC,IG) = PARMSTRT(IC,IG,3) 
              DRATIO (IC,IG) = PARMSTRT(IC,IG,4) 
              DELMAG (IC,IG) = PARMSTRT(IC,IG,5) 
              COLOR  (IC,IG) = PARMSTRT(IC,IG,6) 
              FLUX   (IC,IG) = PARMSTRT(IC,IG,7) 
              FLUXTOT(IC,IG) = PARMSTRT(IC,IG,8) 
          END DO
      END DO
      CALL CLISTERR

 1000 CONTINUE
      STOP
      END