anchor.diff:    write(*,*)
     write(*,'(a)') 'Requesting pressure gradient at 10 points...'
     rc = getpressuregradient(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout3)
  -  do i = 1, 10, 1 
  +  do i = 1, 10 
       write(*,format3) i, ': dpdx=', dataout3(1,i), ', dpdy=', dataout3(2,i), ', dpdz=', dataout3(3,i)
     end do
   
     write(*,*)
     write(*,'(a)') 'Requesting pressure hessian at 10 points...'
     rc = getpressurehessian(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout6)
  -  do i = 1, 10, 1 
  +  do i = 1, 10 
       write(*,format6) 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)
  
positive.code: 
    implicit none
    
    integer, parameter :: rp=4 ! Number of bytes for reals (single precision)
  
    ! ---- Temporal Interpolation Options ----
    integer, parameter :: NoTInt = 0   ! No temporal interpolation
  
negative.code:   character(*), parameter :: authkey = 'edu.jhu.pha.turbulence.testing-201104' // CHAR(0)
  
    integer, parameter :: timestep = 182
    real(RP), parameter :: dt = 0.002_RP
    real(RP), parameter :: time = timestep * dt
    real(RP), parameter :: startTime = 0.285656_RP;
    real(RP), parameter :: endTime = 0.295697_RP;
    real(RP), parameter :: lag_dt = 0.0004_RP     ! fraction of database timestep to use for 
                                                  ! getposition method (Lagrangian integration time step)
    real(RP) :: points(3, 10)    ! input
    real(RP) :: dataout1(10)     ! p
    real(RP) :: dataout3(3, 10)  ! x,y,z
    real(RP) :: dataout4(4, 10)  ! x,y,z,p
    real(RP) :: dataout6(6, 10)  ! results from Pressure Hessian
    real(RP) :: dataout9(9, 10)  ! results from Velocity Gradient
    real(RP) :: dataout18(18, 10) ! results from Velocity Hessian
  
  
    ! Declare the return type of the turblib functions as integer.
  