anchor.diff:      points(3, i) = 0.15 * i 
     end do
   
  -  write(*,*) 'Advancing 10 particles through time'
  -  rc = getposition(authkey, dataset, startTime, endTime, nt, Lag6, PCHIPInt, 10, points, dataout3)
  +  write(*,'(2(a,f8.6),a)') 'Requesting position at 10 points, starting at time ', &
  +       startTime, ' and ending at time ', endTime, '...'
  +  rc = getposition(authkey, dataset, startTime, endTime, lag_dt, Lag6, 10, points, dataout3)
     do i = 1, 10, 1 
  -    write(*,*) i, ': (', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i), ')'
  -  end do
  -
  -  write(*,*) 'Velocity at 10 particle locations'
  -  rc = getvelocity(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout3)
  -  do i = 1, 10, 1 
  -    write(*,*) i, ': (', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i), ')'
  -  end do
  -
  -  write(*,*) 'Forcing at 10 particle locations'
  -  rc = getforce(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout3)
  -  do i = 1, 10, 1 
  -    write(*,*) i, ': (', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i), ')'
  -  end do
  -
  -
  -  write(*,*)
  -  write(*,*) 'Velocity and pressure at 10 particle locations'
  -  rc = getvelocityandpressure(authkey, dataset,  time, Lag6, NoTInt, 10, points, dataout4)
  -  do i = 1, 10, 1 
  -    write(*,*) i, ': (', dataout4(1,i), ', ', dataout4(2,i), ', ', dataout4(3,i), ', ', dataout4(4,i), ')'
  +    write(*,format3) i, ': ', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i)
     end do
   
     write(*,*)
  -  write(*,*) 'Velocity gradient at 10 particle locations'
  +  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(*,*) i, ': (duxdx=', dataout9(1,i), ', duxdy=', dataout9(2,i), &
  +    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),  &
          ', duzdx=', dataout9(7,i), ', duzdy=', dataout9(8,i),  &
  -       ', duzdz=', dataout9(9,i), ')'
  +       ', duzdz=', dataout9(9,i)
     end do
   
  -  write(*,*) 'Velocity laplacian at 10 particle locations'
  +  write(*,*)
  +  write(*,'(a)') 'Requesting velocity laplacian at 10 points...'
     rc = getvelocitylaplacian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout3)
     do i = 1, 10, 1 
  -    write(*,*) i, ': (grad2ux=', dataout3(1,i), ', grad2uy=', dataout3(2,i), ', grad2uz=', dataout3(3,i), ')'
  +    write(*,format3) i, ': grad2ux=', dataout3(1,i), ', grad2uy=', dataout3(2,i), ', grad2uz=', dataout3(3,i)
     end do
   
  -  write(*,*) 'Velocity hessian at 10 particle locations'
  +  write(*,*)
  +  write(*,'(a)') 'Requesting velocity hessian at 10 points...'
     rc = getvelocityhessian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout18)
     do i = 1, 10, 1 
  -    write(*,*) i, ': (d2uxdxdx=', dataout18(1,i), &
  +    write(*,format18) i, ': d2uxdxdx=', dataout18(1,i), &
          ', d2uxdxdy=', dataout18(2,i), &
          ', d2uxdxdz=', dataout18(3,i), &
          ', d2uxdydy=', dataout18(4,i), &
  
positive.code:        ', d2uzdxdz=', dataout18(15,i), &
         ', d2uzdydy=', dataout18(16,i), &
         ', d2uzdydz=', dataout18(18,i), &
         ', d2uzdzdz=', dataout18(18,i), ')'
    end do
  
    write(*,*)
    write(*,*) 'Pressure gradient at 10 particle locations'
    rc = getpressuregradient(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout3)
    do i = 1, 10, 1 
      write(*,*) i, ': (dpdx=', dataout3(1,i), ', dpdy=', dataout3(2,i), ', dpdz=', dataout3(3,i), ')'
    end do
  
    write(*,*)
    write(*,*) 'Pressure hessian at 10 particle locations'
    rc = getpressurehessian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout6)
    do i = 1, 10, 1 
      write(*,*) i, ': (d2pdxdx=', dataout6(1,i), ', d2pdxdy=', dataout6(2,i), &
         ', d2pdxdz=', dataout6(3,i), ', d2pdydy=', dataout6(4,i),  &
         ', d2pdydz=', dataout6(5,i), ', d2pdzdz', dataout6(6,i), ')'
    end do
  
    !
  
negative.code: 
  
  program TurbTest
  
    implicit none
    
    integer, parameter :: rp=4 ! Number of bytes for reals (single precision)
  
    ! ---- Temporal Interpolation Options ----
    integer, parameter :: NoTInt = 0   ! No temporal interpolation
  