anchor.diff:    printcalib: block
       type(source), dimension(:,:), allocatable :: src
       real :: a(5)
  -    real(dp) :: err(100), erra, errb
  +    real(dp) :: err(ndata), erra, errb
   
       allocate( src(size(mapi,1), size(mapi,2)) )
       call mksources(mapi, gti, maph, gth, src)
  
positive.code:     write (34,'(A8,A7,A8,A6,A7,3A9)') 'SCTR', 'ABS', 'HSCAL', 'SIDE', 'BG', &
      & 'DEV', 'D(21.5)', 'D(18.5)'
  
      iterate_parameters: do j = 1, 16000
        call random_number(a)
  
        par % alpha =  M(a(1), 1.0, 10.0) * 1e-6
  
negative.code:       end do iterate_datpoints
        !$omp end parallel do
  
        err = dat % model - dat % sky
        call regr(dat % sky, err, erra, errb)
  
        write (34,'(F8.3,F7.3,F8.1,F6.3,F7.3,3F9.4)') &
          & par % alpha * 1e6, par % relabs, &
          & par % hscale, par % beta, (par % skybg) * 1e7, &
          & sqrt(sum(err**2) / ndata), &
          & errb + erra * 21.5, errb + erra * 18.5
        flush(34)
      end do iterate_parameters
  