PROGRAM BINPLOT
C Plot data and model on screen or into graphics file in the style of
C VISPLOT, but without the VLB stuff.  Reads .CAL files and binary-star
C model files.
C     J.T. Armstrong      12 March 1991

      INCLUDE 'BINFIT.INC'

      INTEGER*4 MXPOINT

      PARAMETER(MXPOINT=100)    ! No. points for model calc; if this 
                                ! changes, change MODVIS, MAKEPLOT too.

      REAL*8 MVIS2(MXPOINT,MXFILT),T(MXPOINT),TMIN,TMAX
      INTEGER*4 NPLOTS
      LOGICAL*1   MODPLOT

      COMMON /MOD_VIS/ MVIS2,T,TMIN,TMAX

      CALL PUTOUT('Program BINPLOT V1.0')
      CALL PUTOUT('Plots visibilities from .CAL files, and models')
      CALL PUTOUT(' ')

C Get the datafile name, read the baseline file and catalog, get the
C star number, read its data from the file(s), calculate frequency 
C band parameters and u-v points.
   10 CALL READFILE
      FRSTRY = .TRUE.
C Find the time limits to be used for the x-axis
   20 CALL TLIMIT(TMIN,TMAX)
C Get the model name, read it, and calculate the model
   30 IF(QUERY(' Plot a model too (Y/N)?')) THEN
          MODPLOT = .TRUE.
      ELSE
          MODPLOT = .FALSE.
      END IF
      IF(MODPLOT) THEN
          CALL CGETMOD
C         Calculate the model visibilities at the times of observation
          CALL BINVIS
          CALL CAGRFAC
C         Calculate the visibilities for the model curves
          CALL MODVIS
      END IF
C Plot the visibilities
      CALL MAKEPLOT(MODPLOT)
C Repeat?
   33 WRITE(OUTC,*) ' Read new file (N), change limits (L),'
      WRITE(OUTC,*) ' repeat plot (R), plot u-v coverage (U),'
      write(outc,*) ' or quit (default): '
      READ(INC,'(A1)') ANSWER
      CALL CAPS(ANSWER)
      IF(ANSWER.EQ.'N') GO TO 10
      IF(ANSWER.EQ.'L') GO TO 20
      IF(ANSWER.EQ.'R') GO TO 30
      if(answer.eq.'U')then
       call buvplot
       goto 33
      endif

      STOP
      END