anchor.diff:              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
  +                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
  
positive.code:             bl%optimal_mean(3) = 0
              bl%optimal_std(3) = 0
          else
              bl%optimal_mean(3) = bl%block_mean(size_e, 1)/bl%block_mean(size_e,2)
  
              bl%optimal_std(3) = fraction_error(bl%block_mean(size_e,1), bl%block_mean(size_e, 2), &
              bl%reblock_data_2(1, size_e,1), bl%block_std(size_e,1), bl%block_std(size_e, 2), &
              bl%block_cov(size_e))
          end if
  
      end subroutine find_optimal_block
  
negative.code:         ! 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
              call change_start(bl, ireport, bl%start_point)
              call mean_std_cov(bl)
              call find_optimal_block(bl)
  