anchor.diff:    ! Enable exit on error.  See README for details.
     CALL turblibSetExitOnError(1)
   
  -  do i = 1, 10, 1
  +  do i = 1, 10
       points(1, i) = 0.20 * i
       points(2, i) = 0.50 * i
       points(3, i) = 0.15 * i 
  
positive.code:   write(*,*)
    write(*,'(a)') 'Requesting velocity at 10 points...'
    rc = getvelocity(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout3)
    do i = 1, 10, 1 
      write(*,format3) i, ': ', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i)
    end do
  
    write(*,*)
    write(*,'(a)') 'Requesting forcing at 10 points...'
    rc = getforce(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout3)
    do i = 1, 10, 1 
      write(*,format3) i, ': ', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i)
    end do
  
    write(*,*)
    write(*,'(a)') 'Requesting velocity and pressure at 10 points...'
    rc = getvelocityandpressure(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout4)
    do i = 1, 10, 1 
      write(*,format4) i, ': ', dataout4(1,i), ', ', dataout4(2,i), ', ', dataout4(3,i), ', ', dataout4(4,i)
    end do
  
    write(*,*)
    write(*,'(a)') 'Requesting velocity gradient at 10 points...'
    rc = getvelocitygradient(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout9)
    do i = 1, 10, 1 
      write(*,format9) i, ': duxdx=', dataout9(1,i), ', duxdy=', dataout9(2,i), &
         ', duxdz=', dataout9(3,i), ', duydx=', dataout9(4,i),  &
         ', duydy=', dataout9(5,i), ', duydz=', dataout9(6,i),  &
  
negative.code: 
    implicit none
    
    integer, parameter :: RP=4 ! Number of bytes for reals (single precision)
  
    ! ---- Temporal Interpolation Options ----
    integer, parameter :: NoTInt = 0   ! No temporal interpolation
  