anchor.diff:    write(*,*)
     write(*,'(a)') 'Requesting velocity at 10 points...'
     rc = getvelocity(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout3)
  -  do i = 1, 10, 1 
  +  do i = 1, 10
       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 
  +  do i = 1, 10 
       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 
  +  do i = 1, 10 
       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 
  +  do i = 1, 10 
       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),  &
  
positive.code:   ! Enable exit on error.  See README for details.
    CALL turblibSetExitOnError(1)
  
    do i = 1, 10, 1
      points(1, i) = 0.20 * i
      points(2, i) = 0.50 * i
      points(3, i) = 0.15 * i 
  
negative.code:   write(*,*)
    write(*,'(a)') 'Requesting velocity laplacian at 10 points...'
    rc = getvelocitylaplacian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout3)
    do i = 1, 10 
      write(*,format3) i, ': grad2ux=', dataout3(1,i), ', grad2uy=', dataout3(2,i), ', grad2uz=', dataout3(3,i)
    end do
  
    write(*,*)
    write(*,'(a)') 'Requesting velocity hessian at 10 points...'
    rc = getvelocityhessian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout18)
    do i = 1, 10 
      write(*,format18) i, ': d2uxdxdx=', dataout18(1,i), &
         ', d2uxdxdy=', dataout18(2,i), &
         ', d2uxdxdz=', dataout18(3,i), &
  