PROGRAM    SPECTRA_BROAD
      IMPLICIT   NONE 
      INTEGER*4  NS
      PARAMETER  ( NS = 10 )
      REAL*8     FREQ(NS), FLUX(NS)
      INTEGER*4  NN
!
      DATA   ( FREQ(NN), FLUX(NN), NN=1,NS ) / &
     &           1.3675D9,    1.670, &
     &           1.6555D9,    1.130, &
     &         230.0D9,       0.108, &
     &          13.621D12,    0.0284, &
     &          24.962D12,    0.0101, &
     &          65.031D12,    0.00398, &
     &         385.830D12,    0.00046, &
     &         481.980D12,    0.00030, &
     &         586.600D12,    0.00031, &
               1294.4D12,     0.00022  &
     &        /
!
      REAL*8     T1(NS), X1(NS)
      INTEGER*4  J1, J2, IUER
!
      DO 410 J1=1,NS
         T1(J1) = DLOG10 ( FREQ(J1) )
         X1(J1) = DLOG10 ( FREQ(J1)*FLUX(J1) )
 410  CONTINUE 
      CALL DIAGI_1 ( NS, T1, X1, IUER )
      END  !#!