anchor.diff:    print*,'          task. These are useful to evaluate how good of model is made. '
     print*,'          This can be done for IQA (self and interaction) and multipole'
     print*,'          models from FEREBUS kriging engine.'
  +  print*,' TASK 5 : Place nodes at high solvent density position in a system. Order the ' 
  +  print*,'          features such that they are labelled relative to a fixed node position'  
  +  print*,'          not labelled as free distingushable molecules. Reduce the features by'    
  +  print*,'          analysing the imporatnce. Cut down the feature list, and keep a backup of' 
  +  print*,'          the original feature list.' 
     print*,'------------------------------------------------------------------------------'
   
     open (status="replace", unit=500, file="Hermes.log", action="write",iostat=ierr)                                ! Unit 40 new output file BONDS
  
positive.code:      print*,'SUCCESS - Hermes.log has been created'
    end if
  
    write(500,'(A)'), 'Hermes Version 1.1, authors J. L. McDonagh, R. Coates and S. Davie (contribution from Dale Stuchfield) '
    write(500,'(A)'), 'Acknowledge the use of adapted routines : '
    write(500,'(A)'), 'strip spaces from Jauch 27 November'
    write(500,'(A)'), 'sort_asend from Rossetta code'
    write(500,'(A)')'------------------------------------ HERMES Version 1.1 -----------------------------------------------'
    write(500,'(A)')' TASK 1 : Determines the structural deviations over a DL POLY MD trajectory or normal mode distortions '
    write(500,'(A)')'          from Tyche. This is sampling statistcs and information. Pre-processing required HISTORY-XYZ.py'
    write(500,'(A)')'          or prep_Tychegjf_for_structure.py. Post-processing MDplot.R over the directories of csv files'
  
negative.code:   character (len=50)   :: TSname
    logical              :: file_ex, file_ex2
  
    print*,'------------------------------------ HERMES Version 1.2 ----------------------'
    print*,' TASK 1 : Determines the structural deviations over a DL POLY MD trajectory or' 
    print*,'          normal mode distortions from Tyche. This is sampling statistcs and '
    print*,'          information. Pre-processing required HISTORY-XYZ.py'
  