anchor.diff:      end select
       
       ! Routine name
  -    if(present(fname))sl_errmsg = sl_errmsg // fname // ': '
  +    if(present(fname))sl_errmsg = trim(sl_errmsg) // ' ' // trim(fname) // ': '
       
       ! Main message
  -    if(present(msg))sl_errmsg = sl_errmsg // msg
  +    if(present(msg))sl_errmsg = trim(sl_errmsg) // trim(msg)
       
       ! Lapack
       zla_msg = ''
       if(present(la_info).and.la_info/=0 .and. &
           & present(la_fname).and.len_trim(la_fname)>0)then
  -        if(len_trim(sl_errmsg)>0)sl_errmsg = sl_errmsg // ': '
  +        if(len_trim(sl_errmsg)>0)sl_errmsg = trim(sl_errmsg) // ': '
           write(zla_msg,*)'Error running LAPACK routine '//&
               & trim(la_fname)//' (info=', la_info, ')'
  -        sl_errmsg = sl_errmsg // zla_msg
  +        sl_errmsg = trim(sl_errmsg) // trim(zla_msg)
       endif
                   
   end function sl_errmsg
  
positive.code:         ziend   = zistart
          zistart = itmp
          if(present(errmsg))&
              & errmsg = errmsg // 'istart > iend => inversion | '
      end if
      if(present(errmsg).and.len_trim(errmsg)>0)then
          if(len_trim(errmsg)>3) errmsg = errmsg(:len_trim(errmsg)-3)        
  
negative.code: 
      ! Covariances and variances
      ! -------------------------
      
      ! Covariances
      allocate(nn(nc,nc))
      allocate(cov(nc,nc))
      cov = 0d0
      nn = 0d0
      call dsyrk('U', trflag, nc, no, 1d0, zvar, nsv, 0d0, cov, nc)
      call dsyrk('U', trflag, nc, no, 1d0, dble(valid(iselects,iselectt)), nsv, 0d0, nn, nc)
      where(nn>0d0) cov = cov / nn
      deallocate(nn)
  
      ! Variances
      zevsums = 0d0
      do i = 1, nsv
          zevsums = zevsums + sum(zvar(i, :)**2)/dble(sum(valid(iselects(i), iselectt)))
      enddo
      if(zuseteof==1)then
          zevsumt = 0d0
          do i = 1, ntv
              zevsumt = zevsumt + sum(zvar(:, i)**2)/dble(sum(valid(iselects, iselectt(i))))
          enddo
      endif
      if(.not.present(pc))deallocate(valid)
  