anchor.diff:    !--------------------------------------------------------------------!
     !                          print final fit                           !
     !--------------------------------------------------------------------!
  -  
  +
     subroutine print_results_csplines()
   
       implicit none
  
positive.code:     write(*,'(1x,"rmsd(mean error)     : ",ES15.8)') rmsd_err
      write(*,'(1x,"(bias                : ",ES15.8,")")') bias
      write(*,*)
    
      select case(ftype)
      case('csplines')
         call print_results_csplines()
  
negative.code:     double precision,                              intent(out) :: dmean_E
      double precision,                              intent(out) :: rmsd_err
  
      double precision, parameter       :: dE_max = 0.05d0 ! Ha
  
      double precision, dimension(np) :: param
      integer                           :: ip, i
  