anchor.diff:          ! For every 50 reports, the optimal mean and standard deviation
           ! and the optimal error in error is calculated and printed.
           if (mod(ireport,50) ==0 .and. ireport >= bl%start_ireport .and. &
  -        bl%start_ireport>0) then
  +                    bl%start_ireport>0) then
               call change_start(bl, ireport, bl%start_point)
               call mean_std_cov(bl)
               call find_optimal_block(bl)
  
positive.code:             do j = 0, (bl%lg_max)
                  B = 2**(j)
  ! [review] - CJCS: From what you said above shouldn't this be B**3 >?
                  if (B > (2*bl%reblock_data_2(1,j,i)&
                      *real(B)*((bl%block_std(j,i)/&
                       bl%block_std(0,i))**4))**(1.0/3.0)) then
                      bl%optimal_size(i) = j
                      exit
  
                  end if
              end do
              if (bl%optimal_size(i) == 1) then
  
negative.code: 
          do i = 1, bl%n_saved_startpoints
              if (mod(bl%n_reports_blocked,(bl%save_fq * i)) == 0 .and. &
                          all(bl%reblock_save(i,:,:,:) == 0)) then
                  bl%reblock_save(i,:,:,:) = bl%reblock_data(:,:,:)
                  bl%product_save(i,:) = bl%data_product(:)
              end if
  