anchor.diff:         ', d2uzdxdz=', dataout18(15,i), &
          ', d2uzdydy=', dataout18(16,i), &
          ', d2uzdydz=', dataout18(18,i), &
  -       ', d2uzdzdz=', dataout18(18,i), ')'
  +       ', d2uzdzdz=', dataout18(18,i)
     end do
   
     write(*,*)
  -  write(*,*) 'Pressure gradient at 10 particle locations'
  +  write(*,'(a)') 'Requesting pressure gradient at 10 points...'
     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), ')'
  +    write(*,format3) i, ': dpdx=', dataout3(1,i), ', dpdy=', dataout3(2,i), ', dpdz=', dataout3(3,i)
     end do
   
     write(*,*)
  -  write(*,*) 'Pressure hessian at 10 particle locations'
  +  write(*,'(a)') 'Requesting pressure hessian at 10 points...'
     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), &
  +    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), ')'
  +       ', d2pdydz=', dataout6(5,i), ', d2pdzdz', dataout6(6,i)
     end do
   
     !
  
positive.code:   ! Currently, only valid datasets are:
    !   'isotropic1024fine' and 'isotropic1024coarse'
    !
    character*100 :: 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)
  
    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
  
  
    ! 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
  
    ! return code
    integer rc
  
    ! loop iterator
    integer i
  
    !
    ! Intialize the gSOAP runtime.
    ! This is required before any WebService routines are called.
  
negative.code:     points(3, i) = 0.15 * i 
    end do
  
    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(*,format3) i, ': ', dataout3(1,i), ', ', dataout3(2,i), ', ', dataout3(3,i)
    end do
  
    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),  &
         ', duzdx=', dataout9(7,i), ', duzdy=', dataout9(8,i),  &
         ', duzdz=', dataout9(9,i)
    end do
  
    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(*,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, 1 
      write(*,format18) i, ': d2uxdxdx=', dataout18(1,i), &
         ', d2uxdxdy=', dataout18(2,i), &
         ', d2uxdxdz=', dataout18(3,i), &
         ', d2uxdydy=', dataout18(4,i), &
  