anchor.diff:    use units,       only: units_initialize,units_open
     use data,        only: data_read,data_average
     use random,      only: random_initialize,random_seq
  -  use map,         only: map_learn
  -  use scoring,     only: gauge,compute_scores,print_scores
  -  use mcmc,        only: mcmc_simulate,mcmc_compute_energy
  +  use mcmc,        only: mcmc_compute_energy
     implicit none
     integer                 :: nvars        ! total number of variables
     integer                 :: nclasses     ! total number of classes
  
positive.code: 
  
    !================================================ print a header
  
    if ( iproc == 0 ) then
       write(ulog,'(a)')         '#'
       write(ulog,'(a)')         '#==========================================='
  
negative.code:   real(kflt)                 :: wid          ! %id for weights calculation
    integer                    :: uprm         ! prm unit
    integer                    :: urst         ! rst unit
    character(len=string_size) :: mode         ! run mode: EVAL, SIM, LEARN
    integer                    :: rseed        ! random seed
    integer                         :: iproc=0,nproc=1
    integer                         :: err
    character(long_string_size)     :: err_string
    integer                         :: dim1,dim2
    integer                         :: uene,ulog
    character(len=long_string_size) :: filename
    integer                         :: j
    real(kflt)                      :: etot,efields,ecouplings
  
    !================================================ set defaults
  
  