anchor.diff:    ! Currently, only valid datasets are:
     !   'isotropic1024fine' and 'isotropic1024coarse'
     !
  -  character*100 :: dataset = 'isotropic1024coarse' // CHAR(0)
  +  character(*), parameter :: dataset = 'isotropic1024coarse' // CHAR(0)
   
     !
     ! Specify your access key here.
     ! If you need one, please visit http://turbulence.pha.jhu.edu/
     ! (We just want to know a bit about our users!)
     !
  -  character*100 :: authkey = 'edu.jhu.pha.turbulence.testing-201104' // CHAR(0)
  +  character(*), parameter :: authkey = 'edu.jhu.pha.turbulence.testing-201104' // CHAR(0)
   
     integer, parameter :: timestep = 182
  -  real :: time = 0.002 * timestep
  -  real :: startTime = 0.002 * timestep
  -  real :: endTime = 0.002 * (timestep + 5)
  -  integer :: nt = 3     ! fraction of database timestep to use for getposition method
  -  real points(3, 10)    ! input
  -  real dataout1(10)     ! p
  -  real dataout3(3, 10)  ! x,y,z
  -  real dataout4(4, 10)  ! x,y,z,p
  -  real dataout6(6, 10)  ! results from Pressure Hessian
  -  real dataout9(9, 10)  ! results from Velocity Gradient
  -  real dataout18(18, 10) ! results from Velocity Hessian
  +  real(rp), parameter :: dt = 0.002_rp
  +  real(rp), parameter :: time = 0.002_rp * timestep 
  +  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.
     ! This is required for custom error handling (see the README).
  -  integer getvelocity, getforce, getvelocityandpressure, getvelocitygradient
  -  integer getvelocitylaplacian, getvelocityhessian
  -  integer getpressuregradient, getpressurehessian
  -  integer getposition
  +  integer :: getvelocity, getforce, getvelocityandpressure, getvelocitygradient
  +  integer :: getvelocitylaplacian, getvelocityhessian
  +  integer :: getpressuregradient, getpressurehessian
  +  integer :: getposition
   
     ! return code
  -  integer rc
  +  integer :: rc
   
     ! loop iterator
  -  integer i
  +  integer :: i
   
  +  ! Formatting rules
  +  character(*), parameter :: format3='(i3,3(a,e13.6))'
  +  character(*), parameter :: format4='(i3,4(a,e13.6))'
  +  character(*), parameter :: format6='(i3,6(a,e13.6))'
  +  character(*), parameter :: format9='(i3,9(a,e13.6))'
  +  character(*), parameter :: format18='(i3,18(a,e13.6))'
     !
     ! Intialize the gSOAP runtime.
     ! This is required before any WebService routines are called.
  
positive.code: 
  
  program TurbTest
    implicit none
  
    ! ---- Temporal Interpolation Options ----
    integer, parameter :: NoTInt = 0   ! No temporal interpolation
  
negative.code:        ', d2uzdxdz=', dataout18(15,i), &
         ', d2uzdydy=', dataout18(16,i), &
         ', d2uzdydz=', dataout18(18,i), &
         ', d2uzdzdz=', dataout18(18,i)
    end do
  
    write(*,*)
    write(*,'(a)') 'Requesting pressure gradient at 10 points...'
    rc = getpressuregradient(authkey, dataset,  time, FD4Lag4, NoTInt, 10, points, dataout3)
    do i = 1, 10, 1 
      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 
      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)
    end do
  
    !
  