anchor.diff:                              timers(6) = timers(6) + (timetwo-timeone)
                           enddo
   
  +                        !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  +                        !!
  +                        !!  Downscale all variables for this point
  +                        !!
  +                        !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                           do v=1,n_obs_variables
                               
                               if (options%debug) then
  
positive.code:         integer(8),dimension(:), intent(inout):: timers
          
          real,    dimension(:),   allocatable :: output
          real,    dimension(:),   allocatable :: obs_analogs
          real,    dimension(:,:), allocatable :: regression_data
          integer, dimension(:),   allocatable :: analogs
          real(8), dimension(:),   allocatable :: coefficients
          
          integer(8)  :: timeone, timetwo
          integer     :: i, j, n, nvars
          integer     :: a, n_analogs, selected_analog, real_analogs
          real        :: rand, analog_threshold
          
          real,    dimension(:,:),allocatable :: threshold_atm 
          real,    dimension(:),  allocatable :: threshold_obs
          logical, dimension(:),  allocatable :: threshold_packing
          integer :: n_packed
          
          n_analogs = options%n_analogs
          analog_threshold = options%analog_threshold
          nvars = size(atm,2)
          n = size(predictor,1)
          
          allocate(coefficients(nvars))
          allocate(output(n))
          if (n_analogs > 0) then
              allocate(analogs(n_analogs))
          endif
          
          if (options%analog_regression) then
              if (n_analogs > 0) then
                  ! if we are using a constant number of analogs just allocate data once
                  ! could also allocate a LARGE size and only use the bottom n...? 
                  allocate(regression_data(n_analogs, nvars))
                  allocate(obs_analogs(n_analogs))
              endif
          endif
  
          if (options%pure_regression) then
              call System_Clock(timeone)
              output(1)  = compute_regression(predictor(1,:), atm, obs, coefficients, errors(1))
  
negative.code:                         enddo
                          
                      endif
                      
                      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                      !!
                      !!  Print a status updates
                      !!
                      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                      !$omp critical (print_lock)
                      current_completed_gridcells = current_completed_gridcells + 1
                      write(*,"(A,f5.1,A$)") char(13), current_completed_gridcells/real(total_number_of_gridcells)*100.0, " %"
                      !$omp end critical (print_lock)
                  enddo
  