anchor.diff:          stop
        end select
     end do args_loop
  +
  +  write(frmt,*) 
   
     call units_open_unf(rst_file,'old',unt,err)
     
  
positive.code:         ! prm file
          iarg = iarg + 1
          call read_opt_arg(iarg,nargs,rst_file,err)
          if (err == 1) then
             write(0,*) 'ERROR ! missing argument: '//trim(arg)//' <filename>'
             stop
          end if
       case default
  
negative.code:      end do
    end if
    k = 0
    write(frmt,*) '(i3,1x,1000f', n_digits + 4, '.', n_digits, ')'
    do iv = 1,nvars
       !write(*,'(i3,1x,1000f8.4)') iv,prm(k+1:k+nclasses)
       !write(*,*) iv,prm(k+1:k+nclasses)
       if (n_digits >= 0) then
          write(*,frmt) iv,prm(k+1:k+nclasses)
       else
          write(*,*) iv,prm(k+1:k+nclasses)
       end if
       k = k + nclasses
    end do
    write(frmt,*) '(2(i3,1x),1000f', n_digits + 4, '.', n_digits, ')'
    do iv = 1,nvars-1
       do jv = iv+1,nvars
          ! order is inverted for printing 
          !write(*,'(2(i3,1x),1000f8.4)') jv,iv,prm(k+1:k+nclasses**2)
          ! (iv, jv) pair of features; inner loop over the states of jv
          if (n_digits >= 0) then
             write(*,frmt) iv,jv,prm(k+1:k+nclasses**2)
          else
             write(*,*) iv,jv,prm(k+1:k+nclasses**2)
          end if
          k = k + nclasses**2
       end do
    end do
  
  