anchor.diff:    INTEGER                                      :: i,j, k
   
     IF (Number_of_bonds == 0 ) THEN
  -     predicted_frequencies = (/ 0.0, 0.0, 0.0 /) 
  +     ! predicted_frequencies has DIMENSION(0) so we can't assign anything to it
     
     ELSEIF (Number_of_bonds == 1) THEN
        predicted_frequencies(1) = typical_frequencies(1)
  
positive.code:   INTEGER                                      :: i,j, k
  
    IF (Number_of_bonds == 0 ) THEN
       predicted_frequencies = (/ 0.0, 0.0, 0.0 /) 
    
    ELSEIF (Number_of_bonds == 1) THEN
  
  
negative.code:   INTEGER                                      :: i,j, k
  
    IF (Number_of_bonds == 0 ) THEN
       ! predicted_frequencies has DIMENSION(0) so we can't assign anything to it
    
    ELSEIF (Number_of_bonds == 1) THEN
  
  