anchor.diff:                                    typical_frequencies(8) ) / 2.0     
        predicted_frequencies(8) = ( typical_frequencies(9) + &
                                     typical_frequencies(10) ) / 2.0
  -     predicted_frequencies(9) = ( typical_frequencies(11) + &
  -                                  typical_frequencies(12) ) / 2.0  
   
     ELSE
        call linspace(typical_frequencies(1),typical_frequencies(2), &
  
positive.code:                                                     795.0,  815.0,  &
                                                     1185.0, 1215.0   /)
  
    REAL(8), DIMENSION(6)  :: peroxy_typ_freq   = (/  470.0,  515.0,  &
                                                      800.0,  875.0,  &
                                                     1100.0, 1170.0   /)
  
    REAL(8), DIMENSION(4) :: Rings_typ_freq    = (/ 2750.0, 2950.0,   &
                                                     900.0, 1100.0    /)
  
  
  
negative.code:   REAL(8), DIMENSION(4)  :: RsCrsR2_typ_freq  = (/  360.0,  370.0,  &
                                                      300.0,  400.0   /)
  
    REAL(8), DIMENSION(4)  :: Alcohol_typ_freq  = (/ 3580.0, 3650.0,  &
                                                     1210.0, 1345.0   /)
  
    REAL(8), DIMENSION(2)  :: Ether_typ_freq    = (/  350.0,  500.0   /)
  
    REAL(8), DIMENSION(8) :: ROOH_typ_freq     = (/ 3580.0, 3650.0,  &
                                                     1300.0, 1320.0,  &
                                                      350.0,  425.0,  &
                                                      825.0,  875.0   /)
  