anchor.diff:              
               call allocate_data(output, n_obs_variables, noutput, nx, ny, n_atm_variables, tr_size, options)
   
  -            call CPU_TIME(master_timeone)
  +            call System_Clock(master_timeone)
               write(*,*), "Entering Parallel Region and Normalizing"
               !$omp parallel default(shared)                  &
               !$omp private(train_data, pred_data, i, j, v, timeone, timetwo)   &
  
positive.code:                 endif
                  
              elseif (options%analog_regression) then
                  call CPU_TIME(inner_timeone)
                  call CPU_TIME(timeone)
                  analogs = find_analogs(predictor(i,:), atm, n_analogs)
                  do a=1,n_analogs
                      obs_analogs(a) = obs(analogs(a))
                      regression_data(a,:) = atm(analogs(a),:)
                  enddo
                  call CPU_TIME(timetwo)
                  timers(1) = timers(1) + (timetwo-timeone)
  
                  call CPU_TIME(timeone)
                  output(i) = compute_regression(predictor(i,:), regression_data, obs_analogs, coefficients, errors(i))
                  output_coeff(1:nvars,i) = coefficients
                  call CPU_TIME(timetwo)
                  timers(2) = timers(2) + (timetwo-timeone)
  
                  call CPU_TIME(timeone)
                  if (logistic_threshold/=kFILL_VALUE) then
                      if (options%logistic_from_analog_exceedance) then
                          
                          ! if the user specified using fewer analogs specifically for the logistic component, find them here
                          if (options%n_log_analogs /= n_analogs) then
                              analogs(1:options%n_log_analogs) = find_analogs(predictor(i,:), atm, options%n_log_analogs)
                              call CPU_TIME(timetwo)
                              timers(9) = timers(9) + (timetwo-timeone)
                              
                              call CPU_TIME(timeone)
                              logistic(i) = compute_analog_exceedance(obs, analogs(1:options%n_log_analogs), logistic_threshold)
                          else
                              logistic(i) = compute_analog_exceedance(obs, analogs, logistic_threshold)
  
negative.code:                     output_coeff(i+nvars,:) = coefficients(i)
                  enddo
              endif
              call System_Clock(timetwo)
              timers(3) = timers(3) + (timetwo-timeone)
  
          endif
          
          do i = 1, n
              if (options%pure_analog) then
                  
  